unified genotyper targetfile
1
0
Entering edit mode
7.4 years ago
Mamta ▴ 500

HI all,

I have a quick question if someone can help.

I am trying to use UnifiedGenotypper for variant calling, however no matter what file I use for -L target interval (I downloaded from Illumina websites, UCSC) but it gives me error

File associated with name gr37Exons.bed is malformed: Problem reading the interval file caused by Badly formed genome loc: Contig chr1 given as location, but this contig isn't present in the Fasta sequence dictionary

Can someone help me with this??

Thanks in advance

Mamta

target_interval gatk unifiedgenotyper • 1.6k views
ADD COMMENT
0
Entering edit mode

Input files reads and reference have incompatible contigs: Found contigs with the same name but different lengths:

##### ERROR   contig reads = 10 / 135524748
##### ERROR   contig reference = 10 / 135534747

This is what I get now. I had used b37 for alignment and same is the ref file.

ADD REPLY
0
Entering edit mode
7.4 years ago
Ram 36k

Maybe your reference sequence has different naming conventions for the chromosomes (1 instead of the expected chr1) - might wanna check that.

If so, change your intervals file to match the naming convention from the ref file

ADD COMMENT
0
Entering edit mode

Thanks, I tried that. But now I get a different error

Input files reads and reference have incompatible contigs: Found contigs with the same name but different lengths:

Any help?

ADD REPLY

Login before adding your answer.

Traffic: 2013 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