from chemsh import *

# this is a demonstrative random sequence
#seg1 = 'dAdAdAdA'
#seg2 = 'dTdTdT'
#seg1 = 'dAdAdAdAdAdAdAdAdA'
seg1 = 'dTdAdTdGdTdAdCdCdAdTdAdGdTdGdAdGdCdAdGdGdCdCdCdGdCdTdTdGdGdGdCdGdGdCdCdCdAdCdAdTdGdCdGdGdGdAdGdTdTdCdCdGdGdAdTdAdGdGdCdTdCdCdAdCdGdAdAdAdTdTdTdAdAdA'
seg1 = 'dAdGdCdTdAdCdCdAdAdAdTdCdG'

#strand1_a = Polynucleotide(seg1, strand2='auto', periodic=10.4, conformation='A')

strand1_b = Polynucleotide(seg1, strand2='auto', periodic=10.4, conformation='B', backbone='drew1981',
#                           alpha=-46.0,
#                           beta=-147.0,
#                           gamma=36.0,
#                           delta=157.0,
#                           epsilon=155.0,
#                           zeta=-96.0,
#                           chi=-102.0
                           )

