Can Some One Explain These Gatk Unifiedgenotyper Commands And Difference In Terms Of Output Also
1
0
Entering edit mode
10.6 years ago
always_learning ★ 1.1k

Can some one explain these GATK Unifiedgenotyper Commands and difference in terms of output also

This one is basically used for calling Dindel target

1) java -Xmx54g -jar GenomeAnalysisTK.jar -T UnifiedGenotyper -nt 16 -l INFO -R refGenomeFile -I fixed.bam -L exome_target_list -o indel.target.vcf -stand_call_conf 50.0 -S SILENT -D ${DBSNP} --max_alternate_alleles 1 -metrics ${resultDir}/${woId}-${sgId}.target.INDELmetrics -glm INDEL -A FisherStrand -A DepthOfCoverage -A AlleleBalance'

This one is used for calling all Dindel:

2) java -Xmx54g -jar GenomeAnalysisTK.jar \ -T UnifiedGenotyper -nt 16 -l INFO -R refGenomeFile \ -I fixed.bam -stand_call_conf 50.0 \ -S SILENT -o indel.vcf \ -D DBSN --max_alternate_alleles 1 \ -metrics INDELmetrics -glm INDEL -A FisherStrand -A DepthOfCoverage -A AlleleBalance'

Please help me in understanding these two !!

gatk variant ngs • 3.2k views
ADD COMMENT
0
Entering edit mode
10.6 years ago

I know this is not what you are looking for but GATK has well written documentation available at http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_genotyper_UnifiedGenotyper.html

You should try to read it first. Also, commands you have mentioned have more than one parameters where they differ. No one would go through your commands and read all the parameters for you and then explain it to you. If you have a specific parameter you want to talk about please mention it below.

ADD COMMENT

Login before adding your answer.

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