Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

skip chromosomes GATK HaplotypeCaller

$
0
0

I am trying to use GATK HaplotypeCaller to identify variants on a set of bam files.
however, my bam files contain only data on chrX, but the headers have the @SQ tag for all chromosomes, including the various _gl000etcetc.
so when I run the tool with

java -jar /GATK/v3.6/GenomeAnalysisTK.jar \
-T HaplotypeCaller \
-R broad_genome/hg19_broad.fasta \
-I myBam.bam \
--genotyping_mode DISCOVERY \
-stand_emit_conf 10 \
-stand_call_conf 30 \
-o ../calls/myBam.vcf

I get the following error:

ERROR MESSAGE: Badly formed genome location: Contig chr1_gl000191_random given as location, but this contig isn't present in the Fasta sequence dictionary

...so, can I tell GATK to ignore all chr BUT chrX?

thank you


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>