SNP calling from pool sequencing with freebayes
1
0
Entering edit mode
7.6 years ago

Hi all

I have NGS data from 5 pools of 5 different breeds. My idea is to find a set of SNPs that identify each breed, then a SNP that is present only in one breed and not in the other four. I used several software and filters but when I visualize the results with IGV or Genome Browser the same SNP is present in more than one breed although with a lower number of reads supporting that SNP. Please could you help me to set the filters with Freebayes? any suggestions will be very appreciated.

greetings

Marco

NGS Freebayes SNP • 2.5k views
ADD COMMENT
1
Entering edit mode
7.6 years ago
h.mon 34k

Have you tried the suggestions on How To Efficiently Remove Snps That Are Present In All Samples? thread?

ADD COMMENT
0
Entering edit mode

Thank you, I saw the suggestions But I think that my problem is a little different. In that suggestion he wants to remove the SNP that are common in all 32 samples, I'm interested to remove the SNPs that are present in at least 2 of my 5 samples, for example

A  B  C
1  2  1
2  3  3
5  5  5
7  8  9
99 99 99

where A, B and C are three different VCF files and I'd like to obtain 7, 8 e 9.

Any suggestion?

ADD REPLY

Login before adding your answer.

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