Entering edit mode
9.8 years ago
darshan.wash
•
0
I was trying get a length of each fragment produce by restriction enzyme for each contig present in my contigs.fasta. In my following example enzyme is AflII. So basically, I am digesting each contig with AflII enzyme. So my code for that is as follows
for record in Bio.SeqIO.parse(open("contigs.fasta"), "fasta"):
printrecord.id, [len(fragment) for fragment in Bio.Restriction.AflII.catalyze(record.seq)])
Now, I have more than one enzyme and I have to generate digestion for that. How can I do it using BioPython?
Please let me know, If you are unclear about the question. PS: I am new to the topic and BioPython as well.