2668,
3186,
3189,
3916,
3947,
4089,
4092,
]
count = 0
of = open(ofn, "w")
for entry in open(mtgs_fn, "U"):
count += 1
if count not in expect_to_fail:
entry = entry.strip()
entry = entry.replace("-", "")
try:
sites = seq2sites(entry)
seq = sites2seq(sites, what="all")
seq = seq.replace("-", "")
if entry != seq:
print count, "FAILED"
of.write("%d\n" % count)
else:
print count, "PASSED"
except:
print count, "FAILED"
of.write("%d\n" % count)
of.flush()
of.close()
here = os.path.abspath(os.path.dirname(__file__))
module_dir = os.path.join(here, '..')
genographic_sites_filepath = os.path.join(here, 'genographic_sites.txt')
import sys
sys.path = [module_dir] + sys.path
outfn = os.path.join(here, 'fail_genographic.txt')
outf = open(outfn, 'w')
from oldowan.mitomotifs import sites2seq
from oldowan.mitomotifs import seq2sites
from oldowan.mitomotifs import str2sites
count = 0
for line in open(genographic_sites_filepath, 'rU'):
count += 1
line = line.upper()
sites = str2sites(line)
sites.sort()
seq = sites2seq(sites, region=range(16000,16570))
roundtrip_sites = seq2sites(seq)
roundtrip_sites.sort()
if sites != roundtrip_sites:
outf.write("%d,%s,%s\n" % (count, sites, roundtrip_sites))
else:
print "%d,good" % count
outf.flush()
outf.close()
DLOOPplus = range(15800,16570) + range(1,1500)
dloops_fn = os.path.join(here, 'dloops.fasta')
ofn = os.path.join(here, 'fail_dloop.txt')
fail_fn = os.path.join(here, 'dloop_expect_to_fail.txt')
expect_to_fail = []
for line in open(fail_fn, 'U'):
expect_to_fail.append(int(line.strip()[:-1]))
of = open(ofn, 'w')
count = 0
for entry in fasta(dloops_fn):
count += 1
if count not in expect_to_fail:
try:
sites = seq2sites(entry["sequence"])
seq = sites2seq(sites, region=DLOOPplus)
seq = seq.replace('-', '')
if entry["sequence"] in seq:
print count, entry["name"], sites2str(sites)
else:
print 'FAILED', count, entry["name"]
of.write('%d, %s\n' % (count, entry["name"]))
except Exception, e:
print 'FAILED',e, count, entry["name"]
of.write('%d, %s\n' % (count, entry["name"]))
of.flush()
of.close()
