Esempio n. 1
0
q = Seq('kjhas')

d = DNAseq('''atgatgaagatcacgatcgccgccctcgctctgctggccgtcgccgttgctgcatacgaa
catgaggactatcattcgcatcccagctacaagttcgagtatggagtaaaggatcctcac
accggagaccacaagagccagtgggaacatcgggacggagatgtcgtcaagggagcatat
acccttcacgaggctgacggaactgagcgagtagttgagtactcgtccgacaagcacaac
ggcttccaggctcacgtcaagcgagtgggtcatgctcaccacccggaagtgtatggacac
catgaggctggacactcgtacggccatggacacggtcatggtcatggcagcagctacgcc
aacggcaacctgtaccagcaccatcactaa
''','testSeq1')

print 'repr:\n%s' % (d)
print 'str(d):\n%s' % (str(d))
print 'len(d):\n%s' % (len(d))
print 'toString:\n%s' % (d.toString())
print 'toFasta:\n%s'  % (d.toFasta())

# ResFreqs and stuff
print 'Freq of "A":\n%s' % (d.getResFreq('A'))
print 'FreqTable:\n'
frqTab = d.getResFreqTable()
for res in frqTab:
    print res


#slice
sd = d.slice(3,9)
print 'slice:\n%s' % (sd.toString())

# rvCmp
rd = d.revCmp()