Mutect2 error - Cannot merge sequence dictionaries
0
0
Entering edit mode
6.8 years ago
asd ▴ 20

I downloaded two BAM files from the TCGA Legacy repository, a normal BAM and a tumor BAM for the same patient. I tried to process the files with Mutect2, but I got an error:

<h5>ERROR AN INPUT FILE ERROR has occurred (version 3.7-0-gcfedb67):</h5> <h5>ERROR MESSAGE: Cannot merge sequence dictionaries because sequence X and 22 are in different orders in two input sequence dictionaries.</h5>

I regenerated the BAM indexes with samtools index and checked the BAMs with picard's validation tool, I got this error.

HISTOGRAM java.lang.String

Error Type Count

ERROR:CIGAR_MAPS_OFF_REFERENCE 16

ERROR:INVALID_INDEXING_BIN 6961

I tried to process the files with hg19 and GRCh38.d1.vd1 reference genomes.

What is the problem and how can I solve it?

mutect2 tcga bam wxs • 1.8k views
ADD COMMENT

Login before adding your answer.

Traffic: 2537 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6