Loading [MathJax]/jax/output/CommonHTML/config.js
前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >专栏 >tablet以及bandage评估

tablet以及bandage评估

作者头像
生信喵实验柴
发布于 2022-05-23 03:33:26
发布于 2022-05-23 03:33:26
1.6K00
代码可运行
举报
文章被收录于专栏:生信喵实验柴生信喵实验柴
运行总次数:0
代码可运行

背景

之前介绍了quast和bosco评估软件,今天带来的是两款可视化拼接的序列的软件,均在windows下运行,得到输入文件即可操作,linux下图形化界面好像也有更新好的,自行探索。

1.tablet 可视化

利用 tablet 可视化检查基因组每个位点覆盖度情况。将测序属于与拼接好的基因组直接进行 bwa 比对。查看每个位点的覆盖情况。通过 tablet 进行可视化扫描。

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
#与拼接结果比对
bwa-mem2 index spades.fa
bwa-mem2 mem -t 12
spades.fa ../data/clean.1.fq.gz ../data/clean.2.fq.gz >spades.sam
#samtools 处理比对结果
samtools sort -@ 12 -o spades.sorted.bam spades.sam
samtools index spades.sorted.bam
#对 spades.fa 建立索引
samtools faidx spades.fa
#tablet 可视化
mkdir tablet
mv spades.fa spades.fa.fai spades.sorted.bam spades.sorted.bam.bai tablet

将4个文件(其中两个是索引)拷贝至 windows 下使用 tablet 可视化。

2.Bandage 可视化

Bandage (a Bioinformatics Application for Navigating De novo Assembly Graphs Easily),是一款用于交互式可视化基因组组装图的软件。可以直接使用例如 Velvet, SPAdes, Trinity和 MEGAHIT)的输出结果,gfa 格式。通过 bandage 可以直接查看有问题的连接区域,改善基因组拼接效果。

Bandage 可以直接进行缩放和平移图形,自定义可视化,搜索序列,提取序列等等操作。

网址:https://rrwick.github.io/Bandage/

使用说明:https://github.com/rrwick/Bandage/wiki

使用spades拼接结果中的fastg和gfa文件导入软件均可以。

写在最后:有时间我们会努力更新的。大家互动交流可以前去论坛,地址在下面,复制去浏览器即可访问,弥补下公众号没有留言功能的缺憾。

代码语言:javascript
代码运行次数:0
运行
AI代码解释
复制
bioinfoer.com

有些板块也可以预设为大家日常趣事的分享等,欢迎大家来提建议。

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2022-04-13,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 生信喵实验柴 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
暂无评论
推荐阅读
编辑精选文章
换一批
数据分析-cuttag分析流程分享1-linux代码流程分析
老板最近比较痴迷于各种seq,由于俩师姐外加一师妹的chip-seq建库老不成功,于是改成了CUTTAG建库,其实读了文献,发现都是相似的原理,只不过是CUTTAG用的细胞起始量远低于chip-seq,同时用的酶不同,建库的时间相对少很多。具体可以参考一下这篇文献CUT&Tag for efficient epigenomic profiling of small samples and single cell(doi:10.1038/s41467-019-09982-5)。
小胡子刺猬的生信学习123
2022/04/03
7.9K1
数据分析-cuttag分析流程分享1-linux代码流程分析
HLA-VBSeq:对全基因组数据进行HLA分型
HLA-VBseq 利用全基因组测序的数据,可以提供8位的HLA分型结果,其文献链接如下
生信修炼手册
2020/05/11
1.7K0
短序列比对练习
针对人基因组RNA-Seq,read跨越外显子,此时要使用支持切除read的软件,例如tophat2,升级的hisat2.而bowtie2则不行。下面我们可以来对比下。
生信喵实验柴
2022/10/25
1.8K0
分享 | ATAC-Seq 分析流程
ATAC-Seq 是“Assay for Transposase-Accessible Chromatin with high-throughput Sequencing”的缩写。 ATAC-Seq 方法依赖于使用高活性转座酶 Tn5 的下一代测序(NGS)文库的构建。将 NGS 接头连接到转座酶上,该转座酶可以使染色质断裂并同时将这些接头整合到开放的染色质区域中。构建的文库可通过 NGS 测序,并使用生物信息学分析具有可及或可访问染色质的基因组区域。
用户7010445
2024/06/18
1.1K0
分享 | ATAC-Seq 分析流程
ATAC-seq经典分析流程(上)
--split-3 把双端测序数据拆分成两个文件,对单端测序数据不起作用.fastq-dump默认会把双端测序结果保存到一个文件里
生信菜鸟团
2023/12/28
2.6K0
ATAC-seq经典分析流程(上)
【直播】我的基因组 43:简单粗糙的WGS数据分析流程
前面我们扯到bam文件的各种操作,vcf文件的各种操作,基础知识不牢固的同学可能已经云里雾里了。这次我们来讲一个简单的。就是拿到了fastq的测序数据,如何把全基因组分析给跑一遍。(不谈细节!) 首先就是fastq文件比对到参考基因组变成sam文件: head -40 read1.fq >tmp/read1.fq head -40 read2.fq >tmp/read2.fq ~/biosoft/bwa/bwa-0.7.15/bwa mem -t 20 -M ~/reference/index/bwa/
生信技能树
2018/03/08
1.9K0
【直播】我的基因组 43:简单粗糙的WGS数据分析流程
snakemake杂记:多个转录组比对到多个基因组得到多个bam文件然后合并
我有10个基因组,然后又12个转录组数据,然后将这个12个基因组数据分别比对到这个10个基因组,每个基因组得到12个bam文件,然后将每个基因组的12个bam文件合并 ,最终得到10个合并的bam文件
用户7010445
2023/08/23
3620
snakemake杂记:多个转录组比对到多个基因组得到多个bam文件然后合并
看优秀本科生如何一周内学会Linux进而搞定RNA-seq上游分析
我是武汉大学基础医学专业第一届的学生,2016年9月刚进大学的时候就选了导师进入实验室接受科研训练。虽然我们实验室不是专门做生物信息学的,但第一次和导师正式交流的时候,她就建议我要学点生信。(巧合的是2016年9月也是生信菜鸟团转型生信技能树的时间点,如果所有的导师都如此明智就好了)
生信技能树
2020/04/14
8.9K1
看优秀本科生如何一周内学会Linux进而搞定RNA-seq上游分析
PSMC分析
A1_R1.fa.gz,A1_R2.fa.gz: 二代重测序clean_reads(个体B)【注:必须两个个体A与B】
花落花相惜
2021/12/06
1.1K0
新冠参考基因组构建
目前新冠病毒的基因组拼接主要采用与参考序列比对,生成一致性序列的方法。所以,参考序列就非常重要,那么参考序列从何而来,参考序列是否准备,遇到新物种如何构建参考序列?
生信喵实验柴
2022/04/07
8140
新冠参考基因组构建
bioconvert计算测序覆盖度并使用R语言ggplot2画折线图进行可视化展示
https://www.ncbi.nlm.nih.gov/nuccore/FN433596
用户7010445
2021/08/31
1.8K0
bioconvert计算测序覆盖度并使用R语言ggplot2画折线图进行可视化展示
拼接病毒基因组
得到参考序列之后,由于病毒基因组较小,不容易从头拼接,因此可以采用与参考序列比对生成一致性序列的方法。无论是 PCR 扩增产物还是宏基因组测序都可以使用该方法得到新冠基因组。
生信喵实验柴
2022/04/07
9710
拼接病毒基因组
HLAreporter : HLA分型软件简介
从示意图可以看出,一个HLA Allel 可以分成四个字段,在加上最后的修饰后缀,共5个字段;在定义HLA 分型结果的分辨率时,会根据分型结果的最大位数来判断,如果只给出了字段一,即血清学分类的信息,代表是2位的分型结果;如果最多给出了字段二,即对应的蛋白信息,代表是4位的分型结果;如果最多给出了字段三,即CDS区信息,代表是8位的分型结果;如果分型结果给出了最后的后缀,代表是9位的分型结果。
生信修炼手册
2020/05/11
2.1K0
生物结构变异分析软件meerkat 0.189使用笔记(一)
meerkat 0.189版本和以前的版本相比,支持bwa mem 输出的bam文件,还支持全外显子数据count SV。
用户1680321
2022/03/22
6800
生物结构变异分析软件meerkat 0.189使用笔记(一)
【直播】我的基因组52:X和Y染色体的同源区域探索
很久以前,我其实就遇到过通过NGS测序数据来判定性别的难题(搜索我博客即可查看详情),本次探究自己的基因组得到的统计结果与常识不符,所以我可以肯定是我们的常识太浅显了。 【直播】我的基因组48:我可能测了一个假的全基因组 【直播】我的基因组49:Y染色体的SNV不能用常规流程来找? 【直播】我的基因组50:从测序深度和位点间距来看SNV分布情况 通过自己的测序数据的详细分析,我才知道PAR(pseudoautosomal region)。这样的X,Y染色体大量同源,说到底是测序片段压根无法准确定位,所以说所
生信技能树
2018/03/08
2K0
【直播】我的基因组52:X和Y染色体的同源区域探索
组装结果纠错
由于三代 nanopore 测序质量比较低,原始数据中存在大量测序错误,即使拼接前进行了纠错,组装结果中仍会存在错误,用长读长或短读长的数据对组装结果进行矫正可以,提高准确率,减少 Miscalls,Indels,改善由错装(mis-assemblies)导致的低比对区域。因此,序列拼接完需要对拼接结果进行优化,根据文献报道,经过 polish 之后,拼接结果与真实基因组(其他测序数据拼接结果)的一致性可以达到 99.99%以上。即使组装工具带有纠错功能,仍建议再次进行一轮或多轮的矫正。
生信喵实验柴
2022/10/25
2K0
组装结果纠错
生信分析过程中这些常见文件的格式以及查看方式你都知道吗?
生信分析过程中,会与很多不同格式的文件打交道,除了原始测序数据fastq之外,还需要准备基因组文件fasta格式和基因注释文件gtf格式。在分析的过程中还会有众多中间文件的生成,如bed、bed12、sam、bam、wig、bigwig、bedgraph等,生成后我们一般会查看下内容了解文件每一列的含义,以此来决定需要提取哪些有用信息列来进行下一步分析。
生信宝典
2019/10/14
2.6K0
生信分析过程中这些常见文件的格式以及查看方式你都知道吗?
序列比对之BWA
BWA是一个用于将DNA序列(特别是低差异性序列)映射到大型参考基因组(例如人类基因组)上的工具。它在基因组学和生物信息学研究中尤为重要。因为它能有效处理高通量测序数据,常常集成于WES分析流程,被广泛应用于基因组学研究,如在寻找与疾病相关的基因变异、理解种系发育关系等领域。
生信菜鸟团
2023/11/22
1.6K0
序列比对之BWA
文献笔记四十三:不同形态的南瓜重测序探索与形态和有价值的农艺性状有关的基因组变异
Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits
用户7010445
2020/03/03
1K0
Sentieon BWA-Meth流程:甲基化WGBS数据分析速度与精确性的完美结合
全基因组甲基化测序(WGBS)是一种研究DNA甲基化的方法,以全面了解在基因组水平上的表观遗传变化。在进行WGBS数据分析时,通常需要使用专门的比对工具,因为这些工具需要能够处理亚硫酸盐转化后的数据。
毅硕科技
2023/04/15
1.1K0
Sentieon BWA-Meth流程:甲基化WGBS数据分析速度与精确性的完美结合
推荐阅读
相关推荐
数据分析-cuttag分析流程分享1-linux代码流程分析
更多 >
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档
本文部分代码块支持一键运行,欢迎体验
本文部分代码块支持一键运行,欢迎体验