SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
CASIM: Variants, Filtering SNP Calls casim UK - Cambridge 0 04-18-2013 01:29 PM
dubious SNP calls from GATK caswater Bioinformatics 0 04-24-2012 06:21 PM
Hi all! [SNP calls from Illumina w/o reference] gdilla Introductions 0 08-13-2010 01:29 PM
maq 2nd and 3rd likely SNP Calls AnamikaDarwin Bioinformatics 0 07-24-2009 09:00 AM
Maq SNP calls from a large pool mimi_lupton Bioinformatics 0 10-31-2008 09:48 AM

Reply
 
Thread Tools
Old 06-23-2015, 12:12 PM   #1
skbrimer
Member
 
Location: OP Kansas

Join Date: Mar 2014
Posts: 53
Question SNP calls with IGV

Good day all!

I have a quick question about SNP calls, more of a best practice I suppose. For my project we are using single end Ion Torrent PGM data for WGS of infectious bronchitis virus. It is a +ssRNA virus so we used the RNAseq pipeline for library prep. The sequenced was aligned to the NCBI RefSeq for IBV using TMAP. In the attached file at 19,095 and 19097 there are two SNPs that are linked. i.e. if its A at 19095 its G at 19097, similarly if its C its T in the other.

The reads are sorted by map quality so the reads at the top are the best mapped to the reference, and when I make the call in the draft I will make the corresponding call as well, but how do you make this call?
Attached Images
File Type: png SNPs.png (30.7 KB, 22 views)
skbrimer is offline   Reply With Quote
Old 06-24-2015, 06:55 AM   #2
Jessica_L
Senior Member
 
Location: Washington, D.C. metro area

Join Date: Feb 2010
Posts: 118
Default

I'm slowly picking up more data analysis and bioinformatics at work, so I've been doing some of the same kind of work lately. Here's how I would look at that:

I see 11 of the A and G reads, 9 of the C and T reads (along with some sequencing errors that are too probably too hard to call). Based on experience with Sanger data, it looks like a heterozygote, but IGV doesn't (as far as I know) really make it easy to see total read depth. To make that call, I'd look at something like the vcf file and check out the allele depth at the indicated SNPs to see how many reads you have of each genotype.

Since I'm still new at this, I'd be interested to see how others would analyze similar data.
Jessica_L is offline   Reply With Quote
Old 06-24-2015, 08:57 AM   #3
skbrimer
Member
 
Location: OP Kansas

Join Date: Mar 2014
Posts: 53
Default

Thank you Jessica, that is a good suggestion. I will give it a try and see if it does shake out for me.
skbrimer is offline   Reply With Quote
Reply

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 01:23 AM.


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