from Bio import AlignIO
from Bio.Seq import Seq
from Bio.Align import MultipleSeqAlignment
from Bio.SeqRecord import SeqRecord
counter = 0
Sitesofvariation = []
for i in range(0,len(alig[0])):
s = alig[:,i]
if s == len(s) * s[0]:
counter += 1
Sitesofvariation.append(i+1)
Will this work to help me find the sites of variation in my genome alignment?
source https://stackoverflow.com/questions/74965637/variable-sites-in-alignment-using-python
Comments
Post a Comment