Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We appreciate your help!

Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

MuTect2 ->Reference name for '1283' not found in sequence dictionary

t1000t1000 Member
edited November 2017 in Ask the GATK team

Hi, I am trying to run
java -jar GenomeAnalysisTK.jar -T MuTect2 -R ucsc.hg19.fasta -I:tumor tumor_rg.bam -I:normal normal_rg.bam -o output.vcf

I am using ucsc.hg19.fasta as a reference and I used:
1. samtools faidx to create the ucsc.hg19.fasta.fai file
2. java -jar picard.jar CreateSequenceDictionary REFERENCE=ucsc.hg19.fasta OUTPUT=ucsc.hg19.fasta.dict

I added a ReadGroup to my two BAM Files:
java -jar picard.jar AddOrReplaceReadGroups I=normal.bam O=normal_rg.bam RGID=4 RGLB=lib1 RGPL=illumina RGPU=unit1 RGSM=20

I also validated the BAM files

Now I get the following error:

##### ERROR MESSAGE: SAM/BAM/CRAM file /analysis/tumnormalor.bam is malformed. Error details: Reference name for '1283' not found in sequence dictionary.

If I exchange the tumor and normal .bam files as command line arguments I get the error
##### ERROR MESSAGE: SAM/BAM/CRAM file /analysis/tumor.bam is malformed. Error details: Reference name for '1283' not found in sequence dictionary.

I tried to google it and found nothing. Can someone give me a clue where to look ?

Is the file ucsc.hg19.fasta.dict the sequence dictionary ?

Tagged:
· Share on Twitter Share on Google+

Answers

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @t1000
    Hi,

    Can you post the BAM file header and the FASTA dict file here? Specifically, I need to see if the @SQ lines in the BAM file match the ones in your dict file.

    Thanks,
    Sheila

    · Share on Twitter Share on Google+
  • michaelmcmichaelmc Member ❠

    Hi,

    I'm having a similar issue with MuTect2 (GATK version 4.0.3.0).

    I'm trying to run:

    gatk Mutect2 \
 -R human_g1k_v37_decoy.fasta \
 -I ${TUMOUR}.recal.bam \
 -tumor ${TUMOUR} \
 -I ${NORMAL}.recal.bam \
 -normal ${NORMAL} \
 --germline-resource af-only-gnomad.raw.sites.b37.vcf \
 --af-of-alleles-not-in-resource 0.00003125 \
 -O ${TUMOUR}_${NORMAL}.vcf

    And I'm getting this error, with a different number after reference name each time:

    [August 8, 2018 3:55:11 PM BST] org.broadinstitute.hellbender.tools.walkers.mutect.Mutect2 done. Elapsed time: 0.14 minutes.
Runtime.totalMemory()=2076049408
java.lang.IllegalArgumentException: Reference name for '3223094' not found in sequence dictionary.
        at htsjdk.samtools.SAMRecord.resolveNameFromIndex(SAMRecord.java:569)
        at htsjdk.samtools.SAMRecord.setReferenceIndex(SAMRecord.java:422)
        at htsjdk.samtools.BAMRecord.<init>(BAMRecord.java:87)
        at htsjdk.samtools.DefaultSAMRecordFactory.createBAMRecord(DefaultSAMRecordFactory.java:42)
        at htsjdk.samtools.BAMRecordCodec.decode(BAMRecordCodec.java:210)
        at htsjdk.samtools.BAMFileReader$BAMFileIterator.getNextRecord(BAMFileReader.java:829)
        at htsjdk.samtools.BAMFileReader$BAMFileIndexIterator.getNextRecord(BAMFileReader.java:981)
        at htsjdk.samtools.BAMFileReader$BAMFileIterator.advance(BAMFileReader.java:803)
        at htsjdk.samtools.BAMFileReader$BAMFileIndexIterator.<init>(BAMFileReader.java:963)
        at htsjdk.samtools.BAMFileReader.createIndexIterator(BAMFileReader.java:918)
        at htsjdk.samtools.BAMFileReader.query(BAMFileReader.java:575)
        at htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter.query(SamReader.java:528)
        at htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter.queryOverlapping(SamReader.java:400)
        at org.broadinstitute.hellbender.utils.iterators.SamReaderQueryingIterator.loadNextIterator(SamReaderQueryingIterator.java:125)
        at org.broadinstitute.hellbender.utils.iterators.SamReaderQueryingIterator.<init>(SamReaderQueryingIterator.java:66)
        at org.broadinstitute.hellbender.engine.ReadsDataSource.prepareIteratorsForTraversal(ReadsDataSource.java:404)
        at org.broadinstitute.hellbender.engine.ReadsDataSource.iterator(ReadsDataSource.java:330)
        at org.broadinstitute.hellbender.engine.MultiIntervalLocalReadShard.iterator(MultiIntervalLocalReadShard.java:134)
        at org.broadinstitute.hellbender.engine.AssemblyRegionIterator.<init>(AssemblyRegionIterator.java:109)
        at org.broadinstitute.hellbender.engine.AssemblyRegionWalker.processReadShard(AssemblyRegionWalker.java:286)
        at org.broadinstitute.hellbender.engine.AssemblyRegionWalker.traverse(AssemblyRegionWalker.java:271)
        at org.broadinstitute.hellbender.engine.GATKTool.doWork(GATKTool.java:893)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.runTool(CommandLineProgram.java:134)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMainPostParseArgs(CommandLineProgram.java:179)
        at org.broadinstitute.hellbender.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:198)
        at org.broadinstitute.hellbender.Main.runCommandLineProgram(Main.java:160)
        at org.broadinstitute.hellbender.Main.mainEntry(Main.java:203)
        at org.broadinstitute.hellbender.Main.main(Main.java:289)

    I've checked the @SQ lines of both the .dict file and the .bam files, but can't see any obvious differences. Here is my .dict file:

    @HD     VN:1.0  SO:unsorted
@SQ     SN:1    LN:249250621    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1b22b98cdeb4a9304cb5d48026a85128
@SQ     SN:2    LN:243199373    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:a0d9851da00400dec1098a9255ac712e
@SQ     SN:3    LN:198022430    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:fdfd811849cc2fadebc929bb925902e5
@SQ     SN:4    LN:191154276    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:23dccd106897542ad87d2765d28a19a1
@SQ     SN:5    LN:180915260    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:0740173db9ffd264d728f32784845cd7
@SQ     SN:6    LN:171115067    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1d3a93a248d92a729ee764823acbbc6b
@SQ     SN:7    LN:159138663    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:618366e953d6aaad97dbe4777c29375e
@SQ     SN:8    LN:146364022    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:96f514a9929e410c6651697bded59aec
@SQ     SN:9    LN:141213431    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:3e273117f15e0a400f01055d9f393768
@SQ     SN:10   LN:135534747    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:988c28e000e84c26d552359af1ea2e1d
@SQ     SN:11   LN:135006516    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:98c59049a2df285c76ffb1c6db8f8b96
@SQ     SN:12   LN:133851895    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:51851ac0e1a115847ad36449b0015864
@SQ     SN:13   LN:115169878    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:283f8d7892baa81b510a015719ca7b0b
@SQ     SN:14   LN:107349540    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:98f3cae32b2a2e9524bc19813927542e
@SQ     SN:15   LN:102531392    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:e5645a794a8238215b2cd77acb95a078
@SQ     SN:16   LN:90354753     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:fc9b1a7b42b97a864f56b348b06095e6
@SQ     SN:17   LN:81195210     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:351f64d4f4f9ddd45b35336ad97aa6de
@SQ     SN:18   LN:78077248     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:b15d4b2d29dde9d3e4f93d1d0f2cbc9c
@SQ     SN:19   LN:59128983     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1aacd71f30db8e561810913e0b72636d
@SQ     SN:20   LN:63025520     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:0dec9660ec1efaaf33281c0d5ea2560f
@SQ     SN:21   LN:48129895     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:2979a6085bfe28e3ad6f552f361ed74d
@SQ     SN:22   LN:51304566     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:a718acaa6135fdca8357d5bfe94211dd
@SQ     SN:X    LN:155270560    UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:7e0e2e580297b7764e31dbc80c2540dd
@SQ     SN:Y    LN:59373566     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1fa3474750af0948bdf97d5a0ee52e51
@SQ     SN:MT   LN:16569        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:c68f52674c9fb33aef52dcf399755519
@SQ     SN:GL000207.1   LN:4262 UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:f3814841f1939d3ca19072d9e89f3fd7
@SQ     SN:GL000226.1   LN:15008        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1c1b2cd1fccbc0a99b6a447fa24d1504
@SQ     SN:GL000229.1   LN:19913        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:d0f40ec87de311d8e715b52e4c7062e1
@SQ     SN:GL000231.1   LN:27386        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:ba8882ce3a1efa2080e5d29b956568a4
@SQ     SN:GL000210.1   LN:27682        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:851106a74238044126131ce2a8e5847c
@SQ     SN:GL000239.1   LN:33824        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:99795f15702caec4fa1c4e15f8a29c07
@SQ     SN:GL000235.1   LN:34474        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:118a25ca210cfbcdfb6c2ebb249f9680
@SQ     SN:GL000201.1   LN:36148        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:dfb7e7ec60ffdcb85cb359ea28454ee9
@SQ     SN:GL000247.1   LN:36422        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:7de00226bb7df1c57276ca6baabafd15
@SQ     SN:GL000245.1   LN:36651        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:89bc61960f37d94abf0df2d481ada0ec
@SQ     SN:GL000197.1   LN:37175        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:6f5efdd36643a9b8c8ccad6f2f1edc7b
@SQ     SN:GL000203.1   LN:37498        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:96358c325fe0e70bee73436e8bb14dbd
@SQ     SN:GL000246.1   LN:38154        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:e4afcd31912af9d9c2546acf1cb23af2
@SQ     SN:GL000249.1   LN:38502        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1d78abec37c15fe29a275eb08d5af236
@SQ     SN:GL000196.1   LN:38914        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:d92206d1bb4c3b4019c43c0875c06dc0
@SQ     SN:GL000248.1   LN:39786        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:5a8e43bec9be36c7b49c84d585107776
@SQ     SN:GL000244.1   LN:39929        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:0996b4475f353ca98bacb756ac479140
@SQ     SN:GL000238.1   LN:39939        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:131b1efc3270cc838686b54e7c34b17b
@SQ     SN:GL000202.1   LN:40103        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:06cbf126247d89664a4faebad130fe9c
@SQ     SN:GL000234.1   LN:40531        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:93f998536b61a56fd0ff47322a911d4b
@SQ     SN:GL000232.1   LN:40652        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:3e06b6741061ad93a8587531307057d8
@SQ     SN:GL000206.1   LN:41001        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:43f69e423533e948bfae5ce1d45bd3f1
@SQ     SN:GL000240.1   LN:41933        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:445a86173da9f237d7bcf41c6cb8cc62
@SQ     SN:GL000236.1   LN:41934        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:fdcd739913efa1fdc64b6c0cd7016779
@SQ     SN:GL000241.1   LN:42152        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:ef4258cdc5a45c206cea8fc3e1d858cf
@SQ     SN:GL000243.1   LN:43341        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:cc34279a7e353136741c9fce79bc4396
@SQ     SN:GL000242.1   LN:43523        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:2f8694fc47576bc81b5fe9e7de0ba49e
@SQ     SN:GL000230.1   LN:43691        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:b4eb71ee878d3706246b7c1dbef69299
@SQ     SN:GL000237.1   LN:45867        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:e0c82e7751df73f4f6d0ed30cdc853c0
@SQ     SN:GL000233.1   LN:45941        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:7fed60298a8d62ff808b74b6ce820001
@SQ     SN:GL000204.1   LN:81310        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:efc49c871536fa8d79cb0a06fa739722
@SQ     SN:GL000198.1   LN:90085        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:868e7784040da90d900d2d1b667a1383
@SQ     SN:GL000208.1   LN:92689        UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:aa81be49bf3fe63a79bdc6a6f279abf6
@SQ     SN:GL000191.1   LN:106433       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:d75b436f50a8214ee9c2a51d30b2c2cc
@SQ     SN:GL000227.1   LN:128374       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:a4aead23f8053f2655e468bcc6ecdceb
@SQ     SN:GL000228.1   LN:129120       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:c5a17c97e2c1a0b6a9cc5a6b064b714f
@SQ     SN:GL000214.1   LN:137718       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:46c2032c37f2ed899eb41c0473319a69
@SQ     SN:GL000221.1   LN:155397       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:3238fb74ea87ae857f9c7508d315babb
@SQ     SN:GL000209.1   LN:159169       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:f40598e2a5a6b26e84a3775e0d1e2c81
@SQ     SN:GL000218.1   LN:161147       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:1d708b54644c26c7e01c2dad5426d38c
@SQ     SN:GL000220.1   LN:161802       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:fc35de963c57bf7648429e6454f1c9db
@SQ     SN:GL000213.1   LN:164239       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:9d424fdcc98866650b58f004080a992a
@SQ     SN:GL000211.1   LN:166566       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:7daaa45c66b288847b9b32b964e623d3
@SQ     SN:GL000199.1   LN:169874       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:569af3b73522fab4b40995ae4944e78e
@SQ     SN:GL000217.1   LN:172149       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:6d243e18dea1945fb7f2517615b8f52e
@SQ     SN:GL000216.1   LN:172294       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:642a232d91c486ac339263820aef7fe0
@SQ     SN:GL000215.1   LN:172545       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:5eb3b418480ae67a997957c909375a73
@SQ     SN:GL000205.1   LN:174588       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:d22441398d99caf673e9afb9a1908ec5
@SQ     SN:GL000219.1   LN:179198       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:f977edd13bac459cb2ed4a5457dba1b3
@SQ     SN:GL000224.1   LN:179693       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:d5b2fc04f6b41b212a4198a07f450e20
@SQ     SN:GL000223.1   LN:180455       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:399dfa03bf32022ab52a846f7ca35b30
@SQ     SN:GL000195.1   LN:182896       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:5d9ec007868d517e73543b005ba48535
@SQ     SN:GL000212.1   LN:186858       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:563531689f3dbd691331fd6c5730a88b
@SQ     SN:GL000222.1   LN:186861       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:6fe9abac455169f50470f5a6b01d0f59
@SQ     SN:GL000200.1   LN:187035       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:75e4c8d17cd4addf3917d1703cacaf25
@SQ     SN:GL000193.1   LN:189789       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:dbb6e8ece0b5de29da56601613007c2a
@SQ     SN:GL000194.1   LN:191469       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:6ac8f815bf8e845bb3031b73f812c012
@SQ     SN:GL000225.1   LN:211173       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:63945c3e6962f28ffd469719a747e73c
@SQ     SN:GL000192.1   LN:547496       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:325ba9e808f669dfeee210fdd7b470ac
@SQ     SN:NC_007605    LN:171823       UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:6743bd63b3ff2b5b8985d8933c53290a
@SQ     SN:hs37d5       LN:35477943     UR:file:/humgen/1kg/reference/human_g1k_v37_decoy.fasta M5:5b6a4b3a81a2d3c134b7d14bf6ad39f1

    and here are the @SQ lines from the ${TUMOUR}.recal.bam file:

    @SQ     SN:1    LN:249250621
@SQ     SN:2    LN:243199373
@SQ     SN:3    LN:198022430
@SQ     SN:4    LN:191154276
@SQ     SN:5    LN:180915260
@SQ     SN:6    LN:171115067
@SQ     SN:7    LN:159138663
@SQ     SN:8    LN:146364022
@SQ     SN:9    LN:141213431
@SQ     SN:10   LN:135534747
@SQ     SN:11   LN:135006516
@SQ     SN:12   LN:133851895
@SQ     SN:13   LN:115169878
@SQ     SN:14   LN:107349540
@SQ     SN:15   LN:102531392
@SQ     SN:16   LN:90354753
@SQ     SN:17   LN:81195210
@SQ     SN:18   LN:78077248
@SQ     SN:19   LN:59128983
@SQ     SN:20   LN:63025520
@SQ     SN:21   LN:48129895
@SQ     SN:22   LN:51304566
@SQ     SN:X    LN:155270560
@SQ     SN:Y    LN:59373566
@SQ     SN:MT   LN:16569
@SQ     SN:GL000207.1   LN:4262
@SQ     SN:GL000226.1   LN:15008
@SQ     SN:GL000229.1   LN:19913
@SQ     SN:GL000231.1   LN:27386
@SQ     SN:GL000210.1   LN:27682
@SQ     SN:GL000239.1   LN:33824
@SQ     SN:GL000235.1   LN:34474
@SQ     SN:GL000201.1   LN:36148
@SQ     SN:GL000247.1   LN:36422
@SQ     SN:GL000245.1   LN:36651
@SQ     SN:GL000197.1   LN:37175
@SQ     SN:GL000203.1   LN:37498
@SQ     SN:GL000246.1   LN:38154
@SQ     SN:GL000249.1   LN:38502
@SQ     SN:GL000196.1   LN:38914
@SQ     SN:GL000248.1   LN:39786
@SQ     SN:GL000244.1   LN:39929
@SQ     SN:GL000238.1   LN:39939
@SQ     SN:GL000202.1   LN:40103
@SQ     SN:GL000234.1   LN:40531
@SQ     SN:GL000232.1   LN:40652
@SQ     SN:GL000206.1   LN:41001
@SQ     SN:GL000240.1   LN:41933
@SQ     SN:GL000236.1   LN:41934
@SQ     SN:GL000241.1   LN:42152
@SQ     SN:GL000243.1   LN:43341
@SQ     SN:GL000242.1   LN:43523
@SQ     SN:GL000230.1   LN:43691
@SQ     SN:GL000237.1   LN:45867
@SQ     SN:GL000233.1   LN:45941
@SQ     SN:GL000204.1   LN:81310
@SQ     SN:GL000198.1   LN:90085
@SQ     SN:GL000208.1   LN:92689
@SQ     SN:GL000191.1   LN:106433
@SQ     SN:GL000227.1   LN:128374
@SQ     SN:GL000228.1   LN:129120
@SQ     SN:GL000214.1   LN:137718
@SQ     SN:GL000221.1   LN:155397
@SQ     SN:GL000209.1   LN:159169
@SQ     SN:GL000218.1   LN:161147
@SQ     SN:GL000220.1   LN:161802
@SQ     SN:GL000213.1   LN:164239
@SQ     SN:GL000211.1   LN:166566
@SQ     SN:GL000199.1   LN:169874
@SQ     SN:GL000217.1   LN:172149
@SQ     SN:GL000216.1   LN:172294
@SQ     SN:GL000215.1   LN:172545
@SQ     SN:GL000205.1   LN:174588
@SQ     SN:GL000219.1   LN:179198
@SQ     SN:GL000224.1   LN:179693
@SQ     SN:GL000223.1   LN:180455
@SQ     SN:GL000195.1   LN:182896
@SQ     SN:GL000212.1   LN:186858
@SQ     SN:GL000222.1   LN:186861
@SQ     SN:GL000200.1   LN:187035
@SQ     SN:GL000193.1   LN:189789
@SQ     SN:GL000194.1   LN:191469
@SQ     SN:GL000225.1   LN:211173
@SQ     SN:GL000192.1   LN:547496
@SQ     SN:NC_007605    LN:171823
@SQ     SN:hs37d5       LN:35477943

    I did use samtools reheader to change something in the header of my ${NORMAL}.recal.bam file (not the @SQ lines), and I notice that in this topic a ReadGroup was added to the bam file. Could this be the issue?

    Thanks,

    Michael

    · Share on Twitter Share on Google+
  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin

    @michaelmc
    Hi Michael,

    Can you confirm this happens with the latest version?

    -Sheila

    · Share on Twitter Share on Google+
  • michaelmcmichaelmc Member ❠

    Yes, I've tried it with 4.0.8.0 and I get the same error

    · Share on Twitter Share on Google+
  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin
    edited September 2018

    @michaelmc
    Hi Michael,

    Thanks. I need to ask the developer and get back to you.

    -Sheila

    EDIT: This has been resolved in this thread.

    Post edited by Sheila on
    · Share on Twitter Share on Google+
Sign In or Register to comment.