This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Generate A Consensus Sequence On The Galaxy Platform

Greetings,

I am attempting to generate a consensus sequence using the galaxy platform. I have been able to generate a pileup using Mpileup however I have been unable to generate a consensus sequence. I have tried using the instructions found in How to generate a consensus fasta sequence from SAM tools pileup? however Galaxy is not recognizing the results of the NGS: SAM Tools -> Pileup-to-Interval as a format that can be utilized for the final suggested step of Fetch Sequences -> Extract Genomic DNA. I have attempted to convert the data again to no avail. I'm at quite a loss at the moment. I would appreciate any insight or tips anyone could provide.

Thank you

galaxy pileup
• 6,627 views • link
updated by Biostar Biostar • written by thelittledoctor • thelittledoctor

I think you are confused about the "extract genomic DNA" tool. You supply a list of intervals and this tool will fetch the reference sequence between these coordinates. For generating a consensus sequence you'll need to take a different approach - one that I'm not sure you'll find in Galaxy.

• 90 views • link
written by Matt Shirley Matt Shirley

1 answer

Generate the pileup from your bam file using "Generate pileup" and select the ten column pileup format. In this pileup, the fourth column corresponds to the consensus base. Using the cut columns 1 and 4 using "Text Manipulation>Cut" and then "Tabular-to-Fasta" you should get to what you want. Post your solution if you get it to work.

• 197 views • link
written by Joseph Hughes Joseph Hughes

This looks like a good answer.

• 0 views • link
written by Matt Shirley Matt Shirley

"Generate pileup" uses "samtools pileup" command which is no longer available on recent versions of Samtools. I wonder if 'samtools mpileup" could be used in a similar way.

• 0 views • link
written by Carlos Borroto Carlos Borroto

Log in to answer this question.

More posts like this
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • How Is Heterozygosity Determined In A Consensus Fasta?
    written by Justin Justin
    <p>I have a consensus fasta file generated by <a href='http://samtools.sourceforge.net/'>samtools</a>/bcftools/vcfutils as discussed in a previous post here: <a href='/p/1388/'>How to generate a consensus fasta sequence …
  • Is It Possible To Generate A Sam Header Without Any External Header Information?
    written by Anima Mundi Anima Mundi
    <p>Hello, I have some <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> files which lack the header. In order to convert them to <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> I need first to add the …
  • Generating Both Pileup And Vcf/Bcf Using Samtools Mpileup
    written by claresloggett • claresloggett
    <p>I'm using <a href='http://samtools.sourceforge.net/'>samtools</a> mpileup and would like to generate both a pileup file and a vcf file as output. I can see how to …
  • How To Pull Out Qual From Multi Line Fastq?
    written by Zhshqzyc • Zhshqzyc
    <p>From <a href="http://biostar.stackexchange.com/questions/1389/how-to-generate-a-consensus-fasta-sequence-from-sam-tools-pileup">http://biostar.stackexchange.com/questions/1389/how-to-generate-a-consensus-fasta-sequence-from-sam-tools-pileup</a> We use awk and perl to get fasta by</p> <pre><code>awk '/^@chr1$/,/^+$/' consensus.fastq | perl -pe "s/@/&gt;/ ; s/\\+//" &gt; chr1.fasta </code></pre> <p>The …
  • Converting Consensus Fastq To Fasta
    written by Love • Love
    <blockquote> <p><strong>Possible Duplicate:</strong> <br/> <a href="http://biostar.stackexchange.com/questions/1389/how-to-generate-a-consensus-fasta-sequence-from-sam-tools-pileup">How to generate a consensus fasta sequence from SAM tools pileup?</a> </p> <p><a href="http://biostar.stackexchange.com/questions/11957/how-common-are-multi-line-fastq-files">How common are multi-line fastq files?</a></p> </blockquote> …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …
  • Fasta File To Sam Conversion
    written by Rallapalli Sukeerthi Teja • Rallapalli Sukeerthi Teja
    <p>Hi,</p> <p>Is it possible to convert a fasta file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a>/bam format using galaxy OR <a href='http://samtools.sourceforge.net/SAM1.pdf'>Sam</a> tools in MAC ? </p> <p>And how …
  • Samtools + Picard Markduplicates
    written by Nery • Nery
    <p>Hello,</p> <p>I'm using <a href='http://samtools.sourceforge.net/'>Samtools</a> to call variants and I am using Picard MarkDuplicates to mark duplicates in my <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. My question is: …
  • How To Generate A Consensus Fasta Sequence From Sam Tools Pileup?
    written by Casey Bergman Casey Bergman
    <p>SAM tools pileup produces a consensus sequence in vertical format, which is not terribly useful for traditional sequence analysis. Does anyone know of a command …