Using UnifiedGenotyper on single chromosome without coordinates
0
0
Entering edit mode
3.2 years ago
shpak.max ▴ 50

I have a large number of bam files, containing subsets that were aligned against (slightly) different reference genomes that were created to account for differences in indel distributions etc across different populations.

I am specifically interested in genetic variation in the X chromosome (corresponding to contig 4 in the assembly), so I attempted to run GATK's UnifiedGenotyper using the argument

-L 4:

without coordinates (because, as described above, the coordinates for this chromosome differ slightly from sample to sample. Running UnifiedGenotyper with this argument (both with and without the ":") returns an error and terminates.

Is there some way to restrict UnifiedGenotyper to a single chromosome/contig without specifying the coordinates?

GATK • 502 views
ADD COMMENT

Login before adding your answer.

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