Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
finding SNPs on PlasmoDB Jocelyn SP Bioinformatics 0 05-20-2014 07:19 PM
Extract reads with defined nucleotide at defined aligned position charlescoldroom Bioinformatics 2 03-22-2014 09:11 AM
Finding SNPs in a particular chromosome shasan Bioinformatics 0 11-01-2013 02:38 PM
finding SNPs using aligments! biouser General 2 08-17-2012 11:05 AM
Finding SNPs Big_Cheese Bioinformatics 9 07-11-2012 08:01 AM

Thread Tools
Old 06-13-2016, 09:03 PM   #1
Junior Member
Location: london

Join Date: Jun 2016
Posts: 2
Default Finding SNPs that have been pre-defined

Hi All,

This is a bit of a beginners question but here goes....

I have a dataset of Illumina paired-end reads mapped to a reference genome. I want to use genome locations of existing SNP's in an input file to find them in my dataset. I've used Bedtools before to find unique SNPs but I'm not sure how to ask it to use an input file with a-priori defined SNPs and output these calls (or even whether bedtools is the best tool for this job). I'd also like to set a min and max read depth if possible.

If anyone has the expertise to point me in the right direction I'd be ever so grateful.

Magpie101 is offline   Reply With Quote
Old 06-15-2016, 06:02 AM   #2
Senior Member
Location: Washington, D.C. metro area

Join Date: Feb 2010
Posts: 116

bedtools would work for that. Check out the intersect command-- you can use a bed or a vcf for locating variants in your dataset. There are several other options, too-- vcfintersect from the vcflib toolkit, etc.
Jessica_L is offline   Reply With Quote


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 07:08 AM.

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