I am very new to dealing with NGS data. Have followed the instructions on how to get paired end data mapped to a reference genome (albeit one that is a related species). Combined all the resulting bam files, cleaned it up but am stumped on how to get the resulting bam file into a fasta consensus. I assume you use samtools and mpile but can't seem to get there. Any suggestions.
Sorry if this is rather basic.