1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100
|
# Copyright 2001 by Iddo Friedberg. All rights reserved.
# This code is part of the Biopython distribution and governed by its
# license. Please see the LICENSE file that should have been included
# as part of this package.
import unittest
from Bio import CAPS
from Bio.Restriction import EcoRI, AluI
from Bio import Alphabet
from Bio.Seq import Seq
from Bio.SeqRecord import SeqRecord
from Bio.Align import MultipleSeqAlignment
def createAlignment(sequences, alphabet):
"""Create an Alignment object from a list of sequences"""
return MultipleSeqAlignment((SeqRecord(Seq(s, alphabet), id="sequence%i"%(i+1))
for (i, s) in enumerate(sequences)),
alphabet)
class TestCAPS(unittest.TestCase):
def test_trivial(self):
enzymes = [EcoRI]
alignment = ["gaattc",
"gaactc",
]
align = createAlignment(alignment, Alphabet.generic_dna)
map = CAPS.CAPSMap(align, enzymes)
self.assertEqual(len(map.dcuts), 1)
self.assertEqual(map.dcuts[0].enzyme, EcoRI)
self.assertEqual(map.dcuts[0].start, 1)
self.assertEqual(map.dcuts[0].cuts_in, [0])
self.assertEqual(map.dcuts[0].blocked_in, [1])
def test(self):
alignment = [
"""\
AAAagaattcTAGATATACCAAACCAGAGAAAACAAATACATAATCGGAGAAATACAGAT
AGAGAGCGAGAGAGATCGACGGCGAAGCTCTTTACCCGGAAACCATTGAAATCGGACGGT
TTAGTGAAAATGGAGGATCAAGTagAtTTTGGGTTCCGTCCGAACGACGAGGAGCTCGTT
GGTCACTATCTCCGTAACAAAATCGAAGGAAACACTAGCCGCGACGTTGAAGTAGCCATC
AGCGAGGTCAACATCTGTAGCTACGATCCTTGGAACTTGCGCTGTAAGTTCCGAATTTTC
""",
"""\
AAAagaTttcTAGATATACCAAACCAGAGAAAACAAATACATAATCGGAGAAATACAGAT
AGAGAGCGAGAGAGATCGACGGCGAAGCTCTTTACCCGGAAACCATTGAAATCGGACGGT
TTAGTGAAAATGGAGGATCAAGTagctTTTGGGTTCCGTCCGAACGACGAGGAGCTCGTT
GGTCACTATCTCCGTAACAAAATCGAAGGAAACACTAGCCGCGACGTTGAAGTAGCCATC
AGCGAGGTCAACATCTGTAGCTACGATCCTTGGAACTTGCGCTGTAAGTTCCGAATTTTC
""",
"""\
AAAagaTttcTAGATATACCAAACCAGAGAAAACAAATACATAATCGGAGAAATACAGAT
AGAGAGCGAGAGAGATCGACGGCGAAGCTCTTTACCCGGAAACCATTGAAATCGGACGGT
TTAGTGAAAATGGAGGATCAAGTagctTTTGGGTTCCGTCCGAACGACGAGGAGCTCGTT
GGTCACTATCTCCGTAACAAAATCGAAGGAAACACTAGCCGCGACGTTGAAGTAGCCATC
AGCGAGGTCAACATCTGTAGCTACGATCCTTGGAACTTGCGCTGTAAGTTCCGAATTTTC
""",
]
enzymes = [EcoRI, AluI]
align = createAlignment(alignment, Alphabet.generic_dna)
map = CAPS.CAPSMap(align, enzymes)
self.assertEqual(len(map.dcuts), 2)
self.assertEqual(map.dcuts[0].enzyme, EcoRI)
self.assertEqual(map.dcuts[0].start, 5)
self.assertEqual(map.dcuts[0].cuts_in, [0])
self.assertEqual(map.dcuts[0].blocked_in, [1, 2])
self.assertEqual(map.dcuts[1].enzyme, AluI)
self.assertEqual(map.dcuts[1].start, 144)
self.assertEqual(map.dcuts[1].cuts_in, [1, 2])
self.assertEqual(map.dcuts[1].blocked_in, [0])
def testNoCAPS(self):
alignment = ["aaaaaaaaaaaaaaaaaaaa",
"aaaaaaaaaaaaaaaaaaaa",
]
enzymes = []
align = createAlignment(alignment, Alphabet.generic_nucleotide)
map = CAPS.CAPSMap(align, enzymes)
self.assertEqual(map.dcuts, [])
def test_uneven(self):
alignment = ["aaaaaaaaaaaaaa",
"aaaaaaaaaaaaaa", # we'll change this below
"aaaaaaaaaaaaaa",
]
align = createAlignment(alignment, Alphabet.generic_nucleotide)
align[1].seq = align[1].seq[:8] # evil
self.assertRaises(CAPS.AlignmentHasDifferentLengthsError,
CAPS.CAPSMap,
align)
if __name__ == "__main__":
runner = unittest.TextTestRunner(verbosity = 2)
unittest.main(testRunner=runner)
|