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
conversion of pileup file to arlequin format catb General 0 06-27-2011 08:04 AM
SAMtools pileup format dong01 Bioinformatics 0 06-24-2011 11:45 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

Reply
 
Thread Tools
Old 08-03-2011, 03:30 AM   #1
Hena
Member
 
Location: Finland

Join Date: Nov 2009
Posts: 19
Default pileup file format

Hi all,

I'm trying to interpret the (m)?pileup file format with regarding to the indels. The documentation says that the indels have pattern [+-][0-9]+[ACGTNacgtn]+. After looking at the format it would seem that at least . or , before the pattern is part of it.

<chr> <pos> C 7 ,,-3aga,,-3aga,-3aga.-3AGA.

Seems to follow that theory pretty well. However then I also have

<chr> <pos> T 15 ,+1g,+1g,+1g,+1g,+1g,+1g,+1g,+1g,+1g,+1g,+1gc+1g,+1g,+1g,+1g

Which does not as the one pair of +1g has been separated with 'c'. How should I interpret that?

I'm using samtools mpileup version 0.1.16, command line:
mpileup -f $GENOME -l position_list.txt -ABQ0 -d 100000 $d/alignment-sorted.bam
Hena 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 07:42 PM.


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