I'm trying to pull out protein sequences from BAM files containing genome alignments. What is the best way to do this? Right now, I have a BAM file containing the neanderthal alignment to the human genome and a fasta file containing the human reference sequence. I wasn't sure whether the next step would be to pull out the neanderthal consensus sequence using SAMtools. If so, where do I go from there, and if that's not the case, what should I do next?