Question: How to filter SNP and Indel data
1
gravatar for ppurkayastha2010
8 months ago by
ppurkayastha201010 wrote:

I have mapped my sequence with a wild type sequence using Bowtie2 and obtained around 8000 SNPs and 700 Indels using Varscan. Please let me know how do I further shortlist the SNPs and Indels (based on which criteria)?

The output file looks like this

POS #CHROM  ID  REF ALT QUAL    FILTER  INFO    FORMAT  Sample1
1   CP009494.1  .   A   C   .   PASS    ADP=16;WT=0;HET=1;HOM=0;NC=0    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    0/1:0:16:16:2:11:68.75%:9.8E-1:32:32:2:0:11:0
2   CP009494.1  .   A   G   .   PASS    ADP=20;WT=0;HET=1;HOM=0;NC=0    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    0/1:0:21:20:6:7:35%:9.8E-1:30:32:6:0:7:0
3   CP009494.1  .   A   G   .   PASS    ADP=21;WT=0;HET=1;HOM=0;NC=0    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    0/1:0:21:21:10:8:36.36%:9.8E-1:30:32:10:0:8:0

Now based on which criteria shall I shortlist the SNPs to a smaller number?? based on reads?

indels snp variant analysis • 243 views
ADD COMMENTlink modified 8 months ago by Guy Reeves1.0k • written 8 months ago by ppurkayastha201010
1
gravatar for Guy Reeves
8 months ago by
Guy Reeves1.0k
Germany
Guy Reeves1.0k wrote:

Bowtie outputs a list of potential variants ( a mixture of sites with low quality evidence for a true variant and those with high quality evidence). You now need to use a program to call variants freebayes and mpileup are common and avalible on usegalaxy.org. They can asses information from multiple individauls simultaniously to call the probability that you have a 'true' varient. The imput file is a BAM file or multiple BAM files. This is a nice trutorial which goes through mpileup and concepts tutorial The trutorial also mentions GATK but would sugest you ignore this as it is a lot of hassle in Galaxy. Personally I use freebayes. Cheers

Guy

ADD COMMENTlink written 8 months ago by Guy Reeves1.0k
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 16.09
Traffic: 79 users visited in the last hour