from pydna.dseqrecord import Dseqrecord
mydna = Dseqrecord("gataaa", circular=True)
mydna
Dseqrecord(o6)
"aga" in mydna
False
concat = mydna._multiply_circular(2)
"aga" in concat
True
"aga" in str(concat.seq).lower()
True
start = str(concat.seq).lower().find("aga")
start
5
concat.add_feature(start, start + len("aga"), label="myfeature")
concat.format()
LOCUS name? 12 bp DNA linear UNK 09-FEB-2017 DEFINITION description? ACCESSION id? VERSION id? KEYWORDS . SOURCE . ORGANISM . . FEATURES Location/Qualifiers misc 6..8 /label="myfeature" ORIGIN 1 gataaagata aa //
from Bio.SeqFeature import SeqFeature, CompoundLocation, FeatureLocation