In [2]:
# Reads a design, prints oligo info, then some strand info.
import cadnano
from cadnano.document import Document
In [3]:
# Read design
app = cadnano.app()
doc = app.document = Document()
doc.readFile('design01.json')
Out[3]:
In [4]:
part = doc.activePart()
oligos = part.oligos()
In [5]:
# No sequence gets applied, so oligo.sequence() should return None.
for oligo in oligos:
print("{0}\t{1}\t\'{2}\'\t{3}".format(oligo, oligo.length(), oligo.getColor(), oligo.sequence()))
In [6]:
# Print Strand info
vhs = list(part.getIdNums()) # convert set to list
for vh_id in vhs[:3]:
fwd_ss, rev_ss = part.getStrandSets(vh_id)
print('VH{0}'.format(vh_id))
print('\t', fwd_ss, '\t', [s.idxs() for s in fwd_ss.strands()], '\n\t\t\t\t',
[s.getColor() for s in fwd_ss.strands()])
print('\t', rev_ss, '\t', [s.idxs() for s in rev_ss.strands()], '\n\t\t\t\t',
[s.getColor() for s in rev_ss.strands()])
In [ ]:
In [ ]: