==== 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脚本。