在使用bwa index chr1.fa
给参考基因组构建索引的时候,遇到问题,日志如下:
[bwa_index] Pack FASTA... 6.32 sec
[bwa_index] Construct BWT for the packed sequence...
[bwa_index] 0.00 seconds elapse.
[bwa_index] Update BWT... 0.00 sec
[bwa_index] Pack forward-only FASTA... 6.15 sec
[bwa_index] Construct SA from BWT and Occ... 0.01 sec
[main] Version: 0.7.17-r1188
[main] CMD: bwa index chr1.fa
[main] Real time: 12.711 sec; CPU: 12.488 sec
结果文件如下:并没有构建起来索引。
-rw-rw-r-- 1 liurj 56 Jul 27 16:13 chr1.fa.sa
-rw-rw-r-- 1 liurj 6 Jul 27 16:13 chr1.fa.amb
-rw-rw-r-- 1 liurj 7 Jul 27 16:13 chr1.fa.ann
-rw-rw-r-- 1 liurj 2 Jul 27 16:13 chr1.fa.pac
-rw-rw-r-- 1 liurj 72 Jul 27 16:13 chr1.fa.bwt
-rw-rw-r-- 1 liurj 2.3G Jul 25 10:36 chr1.fa
尝试使用bwa index -a bwtsw chr1.fa
构建索引,得到以下报错:
[bwa_index] Pack FASTA... 6.40 sec
[bwa_index] Construct BWT for the packed sequence...
Floating point exception (core dumped)
一共12条染色体,1-5超过2个G,6-12少于2个G,1-5没有构建起来而6-12构建起来了,猜测是不是染色体序列超过2个G的关系。但是一个文件已经是一条序列了,除了切割序列还有什么其他的方法可以解决吗,切割序列的话一条染色体又应该怎么处理呢?
染色体大小如下:
-rw-rw-r-- 1 liurj 1.4G Jul 25 10:55 chr12.fa
-rw-rw-r-- 1 liurj 1.6G Jul 25 10:53 chr11.fa
-rw-rw-r-- 1 liurj 1.7G Jul 25 10:52 chr10.fa
-rw-rw-r-- 1 liurj 2.0G Jul 25 10:51 chr9.fa
-rw-rw-r-- 1 liurj 2.0G Jul 25 10:50 chr8.fa
-rw-rw-r-- 1 liurj 2.0G Jul 25 10:49 chr7.fa
-rw-rw-r-- 1 liurj 2.0G Jul 25 10:47 chr6.fa
-rw-rw-r-- 1 liurj 2.1G Jul 25 10:46 chr5.fa
-rw-rw-r-- 1 liurj 2.1G Jul 25 10:44 chr4.fa
-rw-rw-r-- 1 liurj 2.2G Jul 25 10:43 chr3.fa
-rw-rw-r-- 1 liurj 2.2G Jul 25 10:38 chr2.fa
-rw-rw-r-- 1 liurj 2.3G Jul 25 10:36 chr1.fa
1 回答
这家伙很懒,还没有设置简介