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()