Hi all, I'm trying to deal with the restriction module of BioPyhton because I need to create a report of the fragments generated by a double in silico digestion. Briefly, if I digest a .fasta sequence with let's say EcoRI and PstI, I would like to know how many Eco-Pst, Pst-Eco, Eco-Eco and Pst-Pst fragments are generated. Does anybody know how to do it with BioPython?
>Sequence CTAGCGTAATAATAGGTACACTGAATAGTAGTACAGTACGTACAGCTTTTCCTGGGGATC CTATCGCAATCGCGAATGCGACTTCACGTGAATAGATCTCATTCTGAGCTCCCTTATACG TTATAGTTCGACTGTGCTTGATACAAAACGTTTTACTGACTATAACGTGGGGGCACGGGA ATTCAACAGAACTCTCCAAGCTGTCGATTTCTGTATGTTTGAGATTAGATCAGACCTCAC AAGACTCCCTAAACCATCCAGCCCACTTTATATCCCCTCTTCTCCGCCGGAGGTGAATTC AATCCGGCACCAAGGGACTGACAATTTAGCGCAGATACGAGGCAGAACACCGGAAAGACC AGCGGCACTCGCGGGGATCTGGCCCGGTGGGCCCCGGTCCGTGAGCCCGAAGACCCCCTC CCCGAAGATTGGAGGTGCCAGGCAACTGAGGGAGGTGGCTGTCGACGCGCGCCCGGTGCC CGGCCGAGATGTGGGGCCTCCCGGACGGGTCGACCAGCAGCCGGCCGGTGCCCCCTCCGT
In this sequence there are 2 EcoRI sites and 4 TaqI sites. I want to know how many fragments Eco-Taq, Eco-Eco, Taq-Eco or Taq-Taq are made by the digestion.
The answer from ALchEmiXt is something reasonable but I can't undestand the subtraction step: It could be that there are many sites of one enzyme, contigous... how can I deal with that?
Tank you for your help.