售后

编辑本页后请点击“保存”。请参阅 syntax 了解维基语法。只有在您能改进该页面的前提下才编辑它。如果您想尝试一些东西，请先到 playground 热身。

==== snp报错 ====
=== 报错提示 ===
    htsjdk.samtools.SAMException: Unexpected number of metadata chunks 3
    at htsjdk.samtools.BAMIndexMetaData.<init>(BAMIndexMetaData.java:70)
    at htsjdk.samtools.AbstractBAMFileIndex.query(AbstractBAMFileIndex.java:316)
    at htsjdk.samtools.DiskBasedBAMFileIndex.getSpanOverlapping(DiskBasedBAMFileIndex.java:61)
    at htsjdk.samtools.BAMFileReader.getFileSpan(BAMFileReader.java:935)
    at htsjdk.samtools.BAMFileReader.createIndexIterator(BAMFileReader.java:952)
    at htsjdk.samtools.BAMFileReader.query(BAMFileReader.java:612)
    at htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter.query(SamReader.java:533)
    at htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter.query(SamReader.java:538)
    at picard.sam.ReorderSam.doWork(ReorderSam.java:158)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:295)
    at org.broadinstitute.hellbender.cmdline.PicardCommandLineProgramExecutor.instanceMain
    (PicardCommandLineProgramExecutor.java:25)
    at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:162)
    at org.broadinstitute.hellbender.Main.mainEntry(Main.java:205)
    at org.broadinstitute.hellbender.Main.main(Main.java:291)
=== 解决办法 ===
    原因：参考基因组过大（hisat建立索引时的后缀为.ht2l）
    办法：将project.txt中的参考基因组选项选为large，重新刷流程；
    将bam文件重新建立.csi结尾的索引；
    清空snp目录，重新投递snp脚本。

编辑摘要

当您选择开始编辑本页，即寓示你同意将你贡献的内容按下列许可协议发布： CC Attribution-Share Alike 4.0 International

售后

用户工具

站点工具

页面工具