Regarding Gatk output

I want to identify snps in my sequencing data.I am new in this field so i am just following guideline that you provided. i have paired end sequencing data (one sample) and according to documentation. i ran the programs and in the end i ran
java -Xmx32g -jar GenomeAnalysisTK.jar -T HaplotypeCaller -R genome.fasta -I realigned.bam -nct 14 -stand_call_conf 30 -o DK_new.vcf

After getting vcf file, i want to filter vcf file on the basis of allelic depth (DP) that describe in format column. Can you tell me how can i do this.

One more question about the DP value.. why DP value is different in info and format column???
i read the discussion about this and find that because of informative reads and uninformative reads. So it will always possible that allelic depth < total read depth at particular site. so it will better that to filter snps on basis of allelic depth than the total read depth.
I am totally confused that how can i filter snps ????

Please give suggestions about that.

Thanks & Regards


Sign In or Register to comment.