aineniamh / snipit

snipit: summarise snps relative to your reference sequence
GNU General Public License v3.0
148 stars 28 forks source link

Snipit on a.acids #42

Closed woguta closed 2 months ago

woguta commented 3 months ago

I was investigating snps in my sequence data, but got a challenge with amino acids/protein sequences. Seems like snipit only works well with nucleotides (4 bases), can snipit be extended to work with protein sequences as well? I noticed that only amino acids with the same letters as nucleotide bases were the only ones being picked. Merci.

woguta commented 3 months ago

Or has it been implemented?

aineniamh commented 2 months ago

If you specify --sequence-type aa it should run in amino acid mode, but otherwise will default assume your alignment is nucleotides!

woguta commented 2 months ago

Okay, thank you for this. Merci

aineniamh commented 2 months ago

No worries!