SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
pileup output format charlie_sequencing Bioinformatics 5 09-27-2011 08:25 AM
pileup file format Hena Bioinformatics 0 08-03-2011 03:30 AM
samtools pileup format foxyg Bioinformatics 3 11-16-2010 07:25 AM
pileup format- what is ^F EHC Genomic Resequencing 0 07-29-2010 12:01 PM
Samtools pileup SNPs questions Melissa Bioinformatics 2 03-24-2010 10:15 AM

Reply
 
Thread Tools
Old 11-03-2010, 09:53 AM   #1
krobison
Senior Member
 
Location: Boston area

Join Date: Nov 2007
Posts: 747
Default 2 pileup format questions

In the samtools documentation on pileup format, it states
Quote:
A pattern `\+[0-9]+[ACGTNacgtn]+' indicates there is an insertion between this reference position and the next reference position. The length of the insertion is given by the integer in the pattern, followed by the inserted sequence. Here is an example of 2bp insertions on three reads:
seq2 156 A 11 .$......+2AG.+2AG.+2AGGG <975;:<<<<<
Similarly, a pattern `-[0-9]+[ACGTNacgtn]+' represents a deletion from the reference. Here is an exmaple of a 4bp deletions from the reference, supported by two reads:
seq3 200 A 20 ,,,,,..,.-4CACC.-4CACC....,.,,.^~. ==<<<<<<<<<<<::<;2<<
I'm trying to implement a filter on pileup lines & I'm unclear whether insert or deletions have an effect on the quality string. I'm guessing that insert operations have the quality value of the first base of the insert in the quality string, but wanted to double check that.

Also, what do > and < mean in a pileup sequence line (line below was generated by samtools)?

Code:
chr10   53746655        A       15      ><>>>..+1CC...GG^L.^~.  (&:/"P`"R`T````
thanks in advance
krobison is offline   Reply With Quote
Old 11-03-2010, 10:52 AM   #2
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

in the latest online document. >< denote "reference skip", depending on the strand. they are added in samtools-0.1.9. quality corresponding to indels actually means nothing. ignore them.
lh3 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 03:47 AM.


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