标签：toplevel dna fa TAIR10 基因组 bwa Arabidopsis 索引 thaliana

1、当参考基因组

大于2G时（约20亿个核苷酸，20,0000,0000）：bwa index -a bwtsw  xxx.fa

小于2G时（约20亿个核苷酸，20,0000,0000）：bwa index xxx.fa

2、统计参考基因组碱基数目

root@DESKTOP-1N42TVH:/home/test# ls
Arabidopsis_thaliana.TAIR10.dna.toplevel.fa
root@DESKTOP-1N42TVH:/home/test# grep -v "^>" Arabidopsis_thaliana.TAIR10.dna.toplevel.fa | awk 'BEGIN{sum = 0} {sum += length($0)} END {print sum}'
119667750

3、bwa对参考基因组构建索引，生成五个文件

root@DESKTOP-1N42TVH:/home/test# ls
Arabidopsis_thaliana.TAIR10.dna.toplevel.fa
root@DESKTOP-1N42TVH:/home/test# bwa index Arabidopsis_thaliana.TAIR10.dna.toplevel.fa
root@DESKTOP-1N42TVH:/home/test# ls
Arabidopsis_thaliana.TAIR10.dna.toplevel.fa      Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.bwt
Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.amb  Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.pac
Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.ann  Arabidopsis_thaliana.TAIR10.dna.toplevel.fa.sa

标签：toplevel,dna,fa,TAIR10,基因组,bwa,Arabidopsis,索引,thaliana
来源： https://www.cnblogs.com/liujiaxin2018/p/16201230.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。