Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
No output from Samtools pileup pravee1216 Bioinformatics 1 12-05-2012 04:59 AM
output info from samtools pileup? shuang Bioinformatics 1 09-26-2011 05:35 PM
SAMtools pileup output zippered_ohio Bioinformatics 1 09-08-2011 02:27 PM
pileup output different using Maq and Samtools commands mr_boourns Bioinformatics 1 03-31-2011 04:05 AM
samtools pileup output for SNP detection Mansequencer Bioinformatics 1 09-08-2010 07:36 PM

Thread Tools
Old 02-10-2011, 12:54 PM   #1
Location: US

Join Date: Jan 2011
Posts: 13
Default samtools pileup output file

Hey guys, I am using the command line "samtools pileup -cv -f ref_genome snp.sorted.bam" trying to get SNP information between RNA-Seq data and reference genome. The output file has ten columns, which I've no idea for most of them. Here are some of the output:
Group1.1 225532 T G 0 3 60 1 g B
Group1.2 1381414 C Y 0 3 60 1 T$ B
Group1.3 589260 A R 0 3 60 1 ^~g B

From what I can tell so far, the 1st column is the scaffold id and the 2nd is the scaffold coordinate where SNP occurs between the scaffold and the RNA-Seq read. Anyone has previous experiences? Thanks so much!
charlie_sequencing is offline   Reply With Quote
Old 02-10-2011, 01:51 PM   #2
Location: Connecticut

Join Date: Jun 2009
Posts: 74

"pileup up is deprecated, use mpileup." (just trying to be Heng Li, heehee)
However, I remember galaxy has good docs on those fields if the manual wasn't clear for some reason.
csoong is offline   Reply With Quote

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 12:46 AM.

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