| | 24 | |
| | 25 | def test_ceba_bovin(self): |
| | 26 | # ../psCoils.py -f ./ferritin.fasta -l F > ./ferritin.Coils |
| | 27 | weight='un' |
| | 28 | wlambda=0.5 |
| | 29 | win=21 |
| | 30 | modeprofile=False |
| | 31 | seqs=biocomp.pscoils.readFasta("tests/ceba_bovin.fasta") |
| | 32 | modefasta=True |
| | 33 | labels='F' |
| | 34 | params=biocomp.pscoils.Params.Params(weight,win) |
| | 35 | output = StringIO.StringIO() |
| | 36 | for name in seqs.keys(): |
| | 37 | gg,gcc,prob,hept_seq,score=biocomp.pscoils.pred_coil(seqs[name],len(seqs[name]),params,biocomp.pscoils.seqScore) |
| | 38 | biocomp.pscoils.printCoil(seqs[name],len(seqs[name]),hept_seq,score,prob,gcc,gg,labels,io=output) |
| | 39 | bench_output = open("tests/ceba_bovin.Coils").read() |
| | 40 | self.assertEqual(output.getvalue(),bench_output) |