We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

Genotyping 1000 samples 1000x

grafalgrafal MPI TuebingenMember

Dear GATK specialists,

I would like to ask you for advice; I am genotyping combined gvcf for about 1000 samples. For each of them I have chloroplast genome of about 150kb. Problem is that coverage for those chloroplasts reach 500-1000x. When I try to genotype them altogether - I wait 16h and not even a single position is genotyped for the dataset:

My command:

java -Xmx200G -jar /ebio/abt6/rgutaker/bin/GenomeAnalysisTK.jar -T GenotypeGVCFs -R ref.fa -V combined.gvcf.gz -o genotyped.gvcf -nt 5

Would you have any idea how should I genotype that dataset whilst still being able to include all the information for genotyping?

Many thanks.


Best Answer


Sign In or Register to comment.