Entering edit mode
3.6 years ago
minifoog
▴
10
seqList = []
#basic_data_set_aligned contain one letter amino acid sequences
inFile = open("basic_data_set_aligned.fasta",'r')
for record in SeqIO.parse(inFile,'fasta'):
seqList.append(str(record.seq))
m = motifs.create(seqList, alphabet='ACDEFGHIKLMNPQRSTVWYBXZJUO-')
print(m.counts)
m.weblogo("mymotif.png")
Dashes in seqList are my aligners to the other sequences. Is it ideal, for an motif, to include dashes for analysis? Also, the .png file made from using .weblogo is for some reason not recognizing on my computer though I have open many .png files. Any help would be appreciated