Hi, I am new to Python and am trying to write a simple script to loop over sites in a multiple fasta file to identify variant sites. I would like to loop over each residue and then as an output, (1) have a list of variant positions, (2) concatenate residues at SNP sites for each sample (construct a list of haplotypes). Could you point me towards code that allows me to loop over residues in a multiple alignment?
Thank you!
Thank you!
Comment