Question: How to filter SNP and Indel under these criteria?
gravatar for r00628112
6 months ago by
r0062811210 wrote:

Genome-wide polymorphism analysis: after using freebayes to identify SNP and Indel, I will follow these criteria for filtering: (1) the minimum variant frequency of ≥ 90%, (2)average quality of the SNP base ≥ 30 and minimum read depth of 10.

Question1: what kind the tools can be used for these criteria filtering?

Question2: After variation filtering, I'm curious about how to compare the SNPs/Indel difference between cultivars? I analyzed re-sequencing data between four cultivars.

For example, there are A, B, C and D cultivars in the vcf file. But I would like to compare the SNP/Indel variance between A/B and A/C (such as SNP/Indel number)

Thank you very much

ADD COMMENTlink modified 4 months ago by devarora150 • written 6 months ago by r0062811210

bcftools can be used for filtering (on nearly any criteria). Your second question could use some clarification. Do you want to compare allele frequencies between the different cultivars? Or just whether variants were called/not called between each?

ADD REPLYlink modified 6 months ago • written 6 months ago by jared.andrews075.5k

Thanks to your suggestion. Because of limited data, I would like follow this paper to do the analysis of Genome-wide DNA polymorphisms by using four cultivars re-sequencing data. I would like make this kind of table and picture in this paper, thus I'm curious about how to comparing the SNP or Indel number between two cultivars in a combined VCF file.

ADD REPLYlink written 6 months ago by r0062811210
gravatar for devarora
4 months ago by
devarora150 wrote:

You may try GATK filtering and check if this help.

gatk SelectVariants -V vcffile.vcf.gz -select-type SNP -O snps.vcf.gz

gatk SelectVariants -V vcffile.vcf.gz -select-type INDEL -O indels.vcf.gz

ADD COMMENTlink written 4 months ago by devarora150
Please log in to add an answer.


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