Question: Per chromosome gVCF
0
gravatar for win
4.3 years ago by
win810
India
win810 wrote:

I am using the GATK UnifiedGenotyper to generate a gVCF file based using the Emit_All_Stites option. 

The problem I have is that if I have a single chromosome BAM such as BAM for only say chromosome 3 and I want gVCF data for only chromosome 3, is it possible to do that? So I want a gVCF file with variant and non variant sites only for chromosome 3

Right now when I run GATK then it starts processing all the chromosomes. Is there a way to only generate gVCF data for one chromosome even if genome truly means all the chromosomes?

gvcf gatk • 2.0k views
ADD COMMENTlink modified 4.3 years ago by David W4.7k • written 4.3 years ago by win810
2
gravatar for David W
4.3 years ago by
David W4.7k
New Zealand
David W4.7k wrote:

You want the "intervals" argument 'L', as described in the docs

 

ADD COMMENTlink modified 4.3 years ago • written 4.3 years ago by David W4.7k

The intervals argument as shown in docs it text file and so in that case how do i define a single chromosome?

The interval format is specified here

https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_engine_CommandLineGATK.html#--intervals

ADD REPLYlink modified 4.3 years ago • written 4.3 years ago by win810

From the detailed description, linked from the UnifiedGenotyper page 

"You can use samtools-style intervals either explicitly on the command line (e.g. -L chr1 or -L chr1:100-200) or by loading in a file containing a list of intervals (e.g. -L myFile.intervals)."

 

 

ADD REPLYlink written 4.3 years ago by David W4.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1118 users visited in the last hour