SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
N in pileup file Firebird Bioinformatics 9 08-04-2013 10:29 PM
samtool sort: truncated file error? shawpa Bioinformatics 0 12-29-2011 11:21 AM
output format of SAMtool pileup SF_mallish Bioinformatics 3 08-26-2011 03:47 PM
interpretation insertion/deletion file SAMTOOL labmtl Bioinformatics 0 03-28-2011 07:39 AM
SAMtool's pileup format - reference base question mfischer Bioinformatics 11 10-19-2010 04:12 PM

Reply
 
Thread Tools
Old 08-25-2010, 09:35 AM   #1
wuhoucdc
Member
 
Location: Nashville

Join Date: Oct 2009
Posts: 14
Default Samtool Pileup file

Hello All,

Have you seen a pileup file like this?
chr1 862 N G 4 0 3 1 ^$G 8
chr1 863 N C 4 0 3 1 C 8
chr1 864 N A 4 0 3 1 A 8
chr1 865 N A 4 0 3 1 A 9
chr1 866 N G 4 0 3 1 G 8
chr1 867 N C 4 0 3 1 C /
chr1 868 N C 4 0 3 1 C 5
chr1 869 N T 4 0 3 1 T *
chr1 870 N A 8 0 4 2 A^&A 58
chr1 871 N C 8 0 4 2 CC 53

This is ten rows of my file, actually all the data in chr1 look like this.
I used $ samtools pileup -f chr1.fa -c sample.chr1.bam > sample.chr1.raw.txt.

Thank you very much!

Wu
wuhoucdc is offline   Reply With Quote
Old 08-25-2010, 11:36 AM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

That looks about what you'd expect if that area of your reference has a sequence of n's (or, Samtools didn't figure out your reference file properly) which is covered by a single read.

The only surprising part is the last column. That should be one single letter Ascii letter per read crossing that point, and the ascii letter represents the quality score of that exact read letter. But the Ascii ranges usually used today, at least on Illumina instruments, start at @, which is after the numbers. So it's possible that some piece of software might not be scaling your fastq quality scores properly.
swbarnes2 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 11:34 PM.


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