Entering edit mode
9.8 years ago
Jautis
▴
580
Hi, does anybody know of a way to call SNP data from RRBS or another form of bisulfite sequencing? Typical programs to call SNPs have difficulty because RRBS sequences are non-complementary and have mutated base pairs.
The system does not have to be perfect or detect all SNPs; I'm hoping to use the data primarily for mapping local ancestry.
Thank you!
So our lab tried Bis-SNP in the past but it didn't work well. Are there any alternatives out there?
Got here a bit late, but just in case someone is still looking for (faster) alternatives,
https://github.com/zwdzwd/biscuit
https://github.com/hellbelly/BS-Snper
Not that I'm aware of. Calling variants from BSseq data is problematic for a number of reasons, so it's rarely done.
Makes sense; we're just looking for an option to get population-level SNPs in non-model organisms to go along with RRBS data w/o having to pay for 2 rounds of sequencing
The best bang for your buck is going to be to just do RRDNAseq, which is an abbreviation I probably just made up for RRBS sans the bisulfite treatment. Given the quality of the many BSseq datasets I've seen, I wouldn't put that much trust in calling variants from them.