SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools pileup BioSlayer Bioinformatics 9 01-16-2012 03:10 AM
samtools pileup raghu1990 Bioinformatics 0 06-28-2011 02:26 AM
pileup mapping quality cap and SNPs? MatthewVC Bioinformatics 2 06-06-2011 07:29 AM
Pileup shows variants - Not called as SNPs BertieWooster Bioinformatics 10 05-12-2011 10:57 AM
2 pileup format questions krobison Bioinformatics 1 11-03-2010 10:52 AM

Reply
 
Thread Tools
Old 03-23-2010, 11:14 PM   #1
Melissa
Senior Member
 
Location: Switzerland

Join Date: Aug 2008
Posts: 116
Talking Samtools pileup SNPs questions

Hi, I'm using pileup to look for SNPs from bowtie alignment. There are several questions I wanna ask. Your help is greatly appreciated.

1. I observed the maximum SNP quality is 228. Is it correct?

2. $7 is always 60. What's wrong?

3. What does the symbol "$" or "^F" means in $8?

4. Will considering a SNPs site that are covered by reads from plus & minus strand reduce false positive?

5. Also, I noticed BWA maps less reads compared to Bowtie. It gives a lower mapping score, lower SNP quality score and thus, less SNPs. Can I trust these SNPs that are found in bowtie but not BWA alignment?

Cheers,
Melissa
Melissa is offline   Reply With Quote
Old 03-24-2010, 05:12 AM   #2
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

1. it depends on your data

2. in repetitive regions, you will see smaller numbers

3. see FAQ

4. see maq paper. it considers this effect, but probably not as much as it should

5. bowtie does not have a mapping score. in my evaluation, bwa always gives better alignment and better snp calls, but you know my answer is likely to be biased.
lh3 is offline   Reply With Quote
Old 03-24-2010, 10:15 AM   #3
sowmyai
Member
 
Location: America

Join Date: Jan 2010
Posts: 27
Default

I too noticed that all most of my SNPs had quality score of 228. But I had a couple of 255s. Did you ever find out what's magical about 228 ?
sowmyai 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 08:12 AM.


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