SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Frequency of multiple SNPs? Meier1000 Bioinformatics 0 01-11-2018 01:08 AM
Detecting low frequency mutations - GATK spreeth84 Bioinformatics 3 02-13-2015 03:26 AM
How to filter out SNPs with Minor allele frequency less than 5% in VCF file jdpr_100 Bioinformatics 1 10-01-2014 02:14 PM
Detecting SNPs with helioscope lbragg Bioinformatics 3 08-25-2010 01:57 PM

Reply
 
Thread Tools
Old 01-28-2019, 11:44 AM   #1
jleehan
Junior Member
 
Location: Florida

Join Date: Jan 2019
Posts: 2
Default Detecting Known, Low Frequency SNPs in evolving bacterial communities

Hi all,

I'm using NGS data to observe how a community of congenic bacteria evolve over time. Basically, I'm starting by inputting several different alleles of one gene into the same medium and seeing how that population changes over time by sequencing the DNA of that gene.

As far as I can tell, all the SNP calling pipelines appear to be more for genome wide discovery than something like what I am doing.

Currently, I am using VarScan and LoFreq to determine the allele frequencies in this community. The agreement between the two pipelines is decent, but as of now VarScan doesn't seem to be able to detect samples that are present in less than 20% of the population (despite what was reported in this paper and all the other programs are reported to be even worse) and I would like to be able to use at least two different programs to verify the presence of these alleles.

Does anyone know of a better low frequency SNP calling pipeline or perhaps a program that allows you to call SNPs based on known alleles (this would actually be the ideal, but I haven't seen anything like this anywhere that I've read)?

Thanks!

P.S. I'm relatively new to NGS work so please bear with me if I've made some egregious oversight here.
jleehan is offline   Reply With Quote
Old 01-28-2019, 07:52 PM   #2
Meyana
Member
 
Location: Japan

Join Date: Sep 2017
Posts: 37
Default

I have been really pleased with using the variant caller that is a part of the BBTools package for rare variants in mitochondrial DNA.
https://jgi.doe.gov/data-and-tools/b...ls-user-guide/
It is not so well documented at that site, but look here for how to use:
https://github.com/abiswas-odu/Disco...allvariants.sh
or throughout this forum.
Meyana is offline   Reply With Quote
Old 01-29-2019, 07:06 AM   #3
jleehan
Junior Member
 
Location: Florida

Join Date: Jan 2019
Posts: 2
Default

Thank you very much for your recommendation, I'll be sure to check it out!
jleehan is offline   Reply With Quote
Reply

Tags
lofreq, ngs, snp, varscan

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:48 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO