I still haven't found an easy solution to this problem and I am afraid
I'm going to have to write one my own - which makes little sense as I
bet this has been solved thousands of times!
Can anybody point me to a script/software to convert a samtools pileup
file into a fasta consensus file? It would be nice to set coverage
thresholds, etc... but I'll take anything I can work with.
The best google could do for me was this:
Not that helpful,
P.S. If there is a better way of doing this (something other than
samtools) I'm all ears.
One solution, if you want fasta sequence based on the reference genome
(could be a native Galaxy genome, a custom genome in your history, or
really any fasta file in your history as long as the mapped
"chromosomes" names are identical), is to use the tool "NGS: SAM Tools
-> Pileup-to-Interval". Then, to extract fasta sequence based on these
coordinates use the tool "Fetch Sequences -> Extract Genomic DNA".
This utilizes SAMTools, but is in the Galaxy public server and perhaps
this makes it an acceptable option.
If you are interested in examining the variation in your data vs the
reference, please see the tools under "NGS: Indel Analysis". Combined
with the tool "Genome Diversity -> Extract DNA flanking chosen SNPs"
this can incorporate your SNPs into the background reference to
novel fasta sequences.
If still needed, moving from FASTQ to FASTA in Galaxy is very simple
using the tool "NGS: QC and manipulation -> FASTQ to FASTA converter
If command line is your preference, all of Galaxy's tools can be run
there, too, using the source. http://getgalaxy.org
I will post these options at BioStar at the question you quoted, for
that user and others who may have a similar analysis project.
Apologies for the delay in reply. Please let us know if we can help