美文网首页
[Call SNP]bam index坑中坑

[Call SNP]bam index坑中坑

作者: Silver_42ac | 来源:发表于2020-06-02 10:07 被阅读0次

简直就是一个历史遗留问题
GATK 2015年不支持 CSI bam index
GATK 2018 年不支持CSI bam index
GATK 4.0 依然不支持CSI bam index

sort bam 步骤报错:

[E::hts_idx_push] Region 538655237..538655267 cannot be stored in a bai index. Try using a csi index with min_shift = 14, n_lvls >= 6
samtools index: failed to create index for "xxxxxx.bam": Numerical result out of range

GATK 步骤 报错信息:

#因为我用sambamba做了排序,所以没有到上一步sort bam 步骤报错
htsjdk.samtools.SAMException: Unexpected number of metadata chunks 3

问题的开始

htsjdk.samtools.SAMException: Unexpected number of metadata chunks 3
这里是其他软件调用samtools 报错

bai 对于contig长度限制

Question: samtools 1.0 'index' : CSI index vs BAI index ?

image.png

最终方案

Make bam index for long chromosomes
Samtools+bcftools Call SNP

2015年 不支持,之后2018年还是不支持

how to let GATK support Coordinate Sorted Index (CSI) format of bam file

相关文章

网友评论

      本文标题:[Call SNP]bam index坑中坑

      本文链接:https://www.haomeiwen.com/subject/iykqahtx.html