二代重测序运行过程中涉及到的基因组索引文件

本文概述了在基因组比对中,BWA索引文件(如At.fna.*)的重要性,以及GATKcallsSNP/Indel所需的At.dict索引,以及delly软件对fai索引(At.fna.fai)的需求。Picard和samtools用于创建相应的序列字典和索引文件。

1 基因组比对bwa 索引

bwa index At.fna

生成以下5个文件

At.fna.amb
At.fna.ann
At.fna.bwt
At.fna.pac
At.fna.sa

2 GATK call snp indel 需要的索引

GATK call 变异时不仅需要基因组bwa的索引文件,还需要一个.dict的索引文件(注意dict前没有fna 或者 fa 的后缀)

picard CreateSequenceDictionary R=At.fna O=At.dict

3 call SV delly 软件需要的索引文件

delly使用时官方文件说需要基因组的index文件,没说哪种index(我觉得可能只有我有这样的困惑),直到我看见了example的数据文件……
delly运行时需要的是fai的索引文件

samtools faidx At.fna

生成At.fna.fai索引文件

暂时总结这么多,后续遇到再增加。

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值