annotate commons/pyRepetUnit/blastnForClassifierStep1/tests/Test_RepbaseBLRnForClassifierStep1.py @ 31:0ab839023fe4

Uploaded
author m-zytnicki
date Tue, 30 Apr 2013 14:33:21 -0400
parents 94ab73e8a190
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
1 import unittest
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
2 import os
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
3 import shutil
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
4 import ConfigParser
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
5 import sys
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
6 from commons.pyRepetUnit.blastnForClassifierStep1.RepbaseBLRnForClassifierStep1 import RepbaseBLRnForClassifierStep1
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
7 from commons.core.utils.FileUtils import FileUtils
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
8 from commons.core.sql.DbFactory import DbFactory
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
9 from commons.core.sql.TablePathAdaptator import TablePathAdaptator
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
10 from commons.core.coord.Path import Path
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
11 import pyRepet.launcher.programLauncher
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
12
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
13 if not os.environ.has_key( "REPET_PATH" ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
14 print "*** Error: no environment variable REPET_PATH"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
15 sys.exit(1)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
16 sys.path.append( os.environ["REPET_PATH"] )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
17
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
18 NB_EXPECTED_LINES_IN_PATH_FILE = 10
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
19 CURRENT_DIR = os.getcwd()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
20
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
21
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
22 class Test_RepbaseBLRnForClassifierStep1( unittest.TestCase ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
23
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
24
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
25 def setUp( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
26 os.chdir(CURRENT_DIR)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
27 self._inFileName = "dummyFileName"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
28 self._launch_1 = "log = os.system( \""
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
29 self._launch_2 = "\" )\n"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
30 self._launch_2 += "if log != 0:\n"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
31 self._launch_2 += "\tsys.exit(1)\n"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
32 self._cDir = "/home/user/dummy_cdir"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
33 self._tmpDir = "/home/user/dummy_tmpDir"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
34 self._configFileName = "configFile"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
35 self._repbaseName = "dummyRepbase_nt.fa"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
36 f = open(self._configFileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
37 f.write("[repet_env]\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
38 f.write("repet_host: %s\n" % os.environ["REPET_HOST"])
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
39 f.write("repet_user: %s\n" % os.environ["REPET_USER"])
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
40 f.write("repet_pw: %s\n" % os.environ["REPET_PW"])
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
41 f.write("repet_db: %s\n" % os.environ["REPET_DB"])
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
42 f.write("repet_port: %s\n" % os.environ["REPET_PORT"])
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
43 f.write("[detect_features]\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
44 f.write("TE_BLRn: yes\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
45 f.write("TE_nucl_bank: %s\n" % self._repbaseName)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
46 f.write("wublast: yes\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
47 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
48 self._verbose = 0
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
49 self._config = ConfigParser.ConfigParser()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
50 self._config.readfp( open(self._configFileName) )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
51 self._pL = pyRepet.launcher.programLauncher.programLauncher()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
52 self._project = "dummyProject"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
53 self._repbaseBLRn = RepbaseBLRnForClassifierStep1(self._inFileName, self._launch_1, self._launch_2, self._cDir, self._tmpDir, self._configFileName, self._verbose, self._pL, self._project)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
54 self._expFileName = "expFile"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
55 self._bank = self._config.get("detect_features","TE_nucl_bank")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
56 self._db = DbFactory.createInstance()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
57
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
58
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
59 def tearDown( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
60 self._db.dropTable("%s_TE_BLRn_path" % (self._project))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
61 os.chdir(CURRENT_DIR)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
62 self._repbaseBLRn = None
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
63 os.remove(self._configFileName)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
64 if os.path.isfile("%s_BLRn_%s.align.clean_match.path" % (self._project, self._bank)):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
65 os.remove("%s_BLRn_%s.align.clean_match.path" % (self._project, self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
66 if os.path.isfile("%s_BLRn_%s.align.clean_match.path.tmp" % (self._project, self._bank)):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
67 os.remove("%s_BLRn_%s.align.clean_match.path.tmp" % (self._project, self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
68 if os.path.isfile("batch_1.fa_BLRn_%s.*" % (self._bank)):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
69 os.remove("batch_1.fa_BLRn_%s.*" % (self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
70 if os.path.isfile("batch_2.fa_BLRn_%s.*" % (self._bank)):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
71 os.remove("batch_2.fa_BLRn_%s.*" % (self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
72 if os.path.isfile("batch_1.fa_BLRn_%s.align.clean_match.path" % (self._bank)):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
73 os.remove("batch_1.fa_BLRn_%s.align.clean_match.path" % (self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
74 if os.path.isfile("batch_2.fa_BLRn_%s.align.clean_match.path" % (self._bank)):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
75 os.remove("batch_2.fa_BLRn_%s.align.clean_match.path" % (self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
76 if os.path.isfile(self._repbaseName):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
77 os.remove(self._repbaseName)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
78 os.remove(self._repbaseName + "_cut")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
79 os.remove(self._repbaseName + "_cut.xnd")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
80 os.remove(self._repbaseName + "_cut.xns")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
81 os.remove(self._repbaseName + "_cut.xnt")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
82 os.remove(self._repbaseName + ".Nstretch.map")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
83 os.remove("last_time_stamp.log")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
84 if os.path.isdir( "TE_BLRn" ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
85 os.chdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
86 os.chdir(CURRENT_DIR)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
87 shutil.rmtree("TE_BLRn")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
88 if os.path.isfile(self._expFileName):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
89 os.remove(self._expFileName)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
90
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
91
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
92 def test_formatRepbase_ntIfNecessaryWithoutLog( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
93 f = open(self._repbaseName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
94 f.write(">seq1\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
95 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
96 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
97 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
98 f.write(">seq2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
99 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
100 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
101 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
102 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
103 f = open(self._expFileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
104 f.write(">1 seq1 {Cut} 1..180\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
105 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
106 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
107 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
108 f.write(">2 seq2 {Cut} 1..180\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
109 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
110 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
111 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
112 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
113 self._repbaseBLRn.formatRepbase_ntIfNecessary()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
114 obsFileName = self._repbaseName + "_cut"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
115 self.assertTrue(FileUtils.are2FilesIdentical(self._expFileName, obsFileName))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
116
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
117
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
118 def test_formatRepbase_ntIfNecessary_withLogFile( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
119 f = open(self._repbaseName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
120 f.write(">seq1\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
121 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
122 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
123 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
124 f.write(">seq2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
125 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
126 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
127 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
128 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
129 f = open(self._expFileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
130 f.write(">1 seq1 {Cut} 1..180\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
131 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
132 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
133 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
134 f.write(">2 seq2 {Cut} 1..180\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
135 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
136 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTTTTTTTGCGCGTGGCGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
137 f.write("ATGCGTGCGTAAATGCGTAATGCGTAAATGCGTAAATTTGCGCGTAAAGTATGCGTGCAT\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
138 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
139 repbaseBLRnWithLog = RepbaseBLRnForClassifierStep1(self._inFileName, self._launch_1, self._launch_2, self._cDir, self._tmpDir, self._configFileName, self._verbose, self._pL, self._project)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
140 repbaseBLRnWithLog.formatRepbase_ntIfNecessary()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
141 obsFileName = self._repbaseName + "_cut"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
142 self.assertTrue(FileUtils.are2FilesIdentical(self._expFileName, obsFileName))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
143
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
144
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
145 def test_createCmdToLaunch( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
146 bank = self._config.get("detect_features","TE_nucl_bank")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
147 obsCmd = self._repbaseBLRn.createCmdToLaunch()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
148 expCmd = self._launch_1 + os.environ["REPET_PATH"] + "/bin/blaster"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
149 expCmd += " -q %s" % ( self._inFileName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
150 expCmd += " -s %s/%s" % ( self._cDir, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
151 expCmd += " -B %s_BLRn_%s" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
152 expCmd += " -n blastn"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
153 if self._config.get("detect_features","wublast") == "yes":
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
154 expCmd += " -W"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
155 expCmd += " -r"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
156 expCmd += " -v 1"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
157 expCmd += self._launch_2
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
158 expCmd += "if not os.path.exists( \"%s/%s_BLRn_%s.param\" ):\n" % ( self._cDir, self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
159 expCmd += "\tos.system( \"mv %s_BLRn_%s.param %s\" )\n" % ( self._inFileName, bank, self._cDir )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
160 expCmd += "if os.path.exists( \"%s_cut\" ):\n" % ( self._inFileName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
161 expCmd += "\tos.system( \"rm -f %s_cut*\" )\n" % ( self._inFileName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
162 expCmd += "if os.path.exists( \"%s.Nstretch.map\" ):\n" % ( self._inFileName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
163 expCmd += "\tos.remove( \"%s.Nstretch.map\" )\n" % ( self._inFileName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
164 expCmd += "if os.path.exists( \"%s_BLRn_%s.raw\" ):\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
165 expCmd += "\tos.remove( \"%s_BLRn_%s.raw\" )\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
166 expCmd += "if os.path.exists( \"%s_BLRn_%s.seq_treated\" ):\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
167 expCmd += "\tos.remove( \"%s_BLRn_%s.seq_treated\" )\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
168 expCmd += self._launch_1
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
169 expCmd += os.environ["REPET_PATH"] + "/bin/matcher"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
170 expCmd += " -m %s_BLRn_%s.align" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
171 expCmd += " -q %s" % ( self._inFileName )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
172 expCmd += " -s %s/%s" % ( self._cDir, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
173 expCmd += " -j"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
174 expCmd += " -v 1"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
175 expCmd += self._launch_2
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
176 expCmd += "if not os.path.exists( \"%s/%s_BLRn_%s.align.clean_match.path\" ):\n" % ( self._cDir, self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
177 expCmd += "\tos.system( \"mv %s_BLRn_%s.align.clean_match.path %s\" )\n" % ( self._inFileName, bank, self._cDir )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
178 expCmd += "if not os.path.exists( \"%s/%s_BLRn_%s.align.clean_match.param\" ):\n" % ( self._cDir, self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
179 expCmd += "\tos.system( \"mv %s_BLRn_%s.align.clean_match.param %s\" )\n" % ( self._inFileName, bank, self._cDir )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
180 expCmd += "if os.path.exists( \"%s_BLRn_%s.align\" ):\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
181 expCmd += "\tos.remove( \"%s_BLRn_%s.align\" )\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
182 expCmd += "if os.path.exists( \"%s_BLRn_%s.align.clean_match.fa\" ):\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
183 expCmd += "\tos.remove( \"%s_BLRn_%s.align.clean_match.fa\" )\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
184 expCmd += "if os.path.exists( \"%s_BLRn_%s.align.clean_match.map\" ):\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
185 expCmd += "\tos.remove( \"%s_BLRn_%s.align.clean_match.map\" )\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
186 expCmd += "if os.path.exists( \"%s_BLRn_%s.align.clean_match.tab\" ):\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
187 expCmd += "\tos.remove( \"%s_BLRn_%s.align.clean_match.tab\" )\n" % ( self._inFileName, bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
188 if self._tmpDir != self._cDir:
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
189 expCmd += "if os.path.exists( \"%s\" ):\n" % ( bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
190 expCmd += "\tos.remove( \"%s\" )\n" % ( bank )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
191 self.assertEquals( expCmd, obsCmd )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
192
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
193
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
194 def test_concatPathFile( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
195 os.chdir(CURRENT_DIR)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
196 if not FileUtils.isRessourceExists("TE_BLRn"):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
197 os.mkdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
198 self._createPathFiles(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
199 os.chdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
200 f = open(self._expFileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
201 f.write("1\tQueryName1\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
202 f.write("2\tQueryName1\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
203 f.write("3\tQueryName2\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
204 f.write("4\tQueryName3\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
205 f.write("5\tQueryName1\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
206 f.write("1\tQueryName4\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
207 f.write("2\tQueryName4\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
208 f.write("3\tQueryName5\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
209 f.write("4\tQueryName6\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
210 f.write("5\tQueryName7\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
211 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
212 FileUtils.sortFileContent(self._expFileName)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
213 obsFileName = "%s_BLRn_%s.align.clean_match.path.tmp" % (self._project, self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
214 self._repbaseBLRn._concatPathFile(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
215 FileUtils.sortFileContent(obsFileName)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
216 self.assertTrue(FileUtils.are2FilesIdentical(self._expFileName, obsFileName))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
217
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
218
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
219 def test_adaptIDInPathFile( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
220 f = open("%s_BLRn_%s.align.clean_match.path.tmp" % (self._project, self._bank), "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
221 f.write("1\tQueryName1\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
222 f.write("2\tQueryName1\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
223 f.write("3\tQueryName2\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
224 f.write("4\tQueryName3\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
225 f.write("5\tQueryName1\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
226 f.write("1\tQueryName4\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
227 f.write("2\tQueryName4\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
228 f.write("3\tQueryName5\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
229 f.write("4\tQueryName6\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
230 f.write("5\tQueryName7\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
231 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
232 self._repbaseBLRn._adaptIDInPathFile(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
233 f = open(self._expFileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
234 f.write("1\tQueryName1\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
235 f.write("2\tQueryName1\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
236 f.write("3\tQueryName2\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
237 f.write("4\tQueryName3\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
238 f.write("5\tQueryName1\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
239 f.write("6\tQueryName4\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
240 f.write("7\tQueryName4\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
241 f.write("8\tQueryName5\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
242 f.write("9\tQueryName6\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
243 f.write("10\tQueryName7\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
244 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
245 obsFileName = "%s_BLRn_%s.align.clean_match.path" % (self._project, self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
246 self.assertTrue(FileUtils.are2FilesIdentical(self._expFileName, obsFileName))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
247
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
248
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
249 def test_loadPathFileInTable( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
250 os.chdir(CURRENT_DIR)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
251 if not FileUtils.isRessourceExists("TE_BLRn"):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
252 os.mkdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
253 self._createPathFiles(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
254 os.chdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
255 f = open("%s_BLRn_%s.align.clean_match.path" % (self._project, self._bank), "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
256 f.write("1\tQueryName1\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
257 f.write("2\tQueryName1\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
258 f.write("3\tQueryName2\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
259 f.write("4\tQueryName3\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
260 f.write("5\tQueryName1\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
261 f.write("6\tQueryName4\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
262 f.write("7\tQueryName4\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
263 f.write("8\tQueryName5\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
264 f.write("9\tQueryName6\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
265 f.write("10\tQueryName7\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
266 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
267 self._repbaseBLRn._loadPathFileInTable(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
268 resultFromFilePathList = self._readPathResultsFromFileAndFillList()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
269 resultFromTablePathList = self._readPathResultsFromTableAndFillList()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
270 self.assertEquals(resultFromFilePathList, resultFromTablePathList)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
271
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
272
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
273 def test_findAndRemoveUselessFiles( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
274 self._createFile("%s_BLRn_%s.align.clean_match.path.tmp" % (self._project, self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
275 self._createFile("batch_1.fa_BLRn_%s.*" % (self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
276 self._createFile("batch_2.fa_BLRn_%s.*" % (self._bank))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
277 self._repbaseBLRn._findAndRemoveUselessFiles(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
278 self.assertFalse(FileUtils.isRessourceExists("%s_BLRn_%s.align.clean_match.path.tmp" % (self._project, self._bank)))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
279 self.assertFalse(FileUtils.isRessourceExists("batch_1.fa_BLRn_%s.*" % (self._bank)))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
280 self.assertFalse(FileUtils.isRessourceExists("batch_2.fa_BLRn_%s.*" % (self._bank)))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
281
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
282
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
283 def test_collectRepbaseBLRn( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
284 os.chdir(CURRENT_DIR)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
285 if not FileUtils.isRessourceExists("TE_BLRn"):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
286 os.mkdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
287 self._createPathFiles(self._bank)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
288 os.chdir( "TE_BLRn" )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
289 self._repbaseBLRn.collectRepbaseBLRn()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
290 self.assertTrue(FileUtils.isRessourceExists("%s_BLRn_%s.align.clean_match.path" % ( self._project, self._bank )))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
291 self.assertEquals(NB_EXPECTED_LINES_IN_PATH_FILE, FileUtils.getNbLinesInSingleFile("%s_BLRn_%s.align.clean_match.path" % ( self._project, self._bank )))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
292 self.assertTrue(self._db.doesTableExist("%s_TE_BLRn_path" % (self._project)))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
293 self._db.execute('select * from %s_TE_BLRn_path' % (self._project))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
294 self.assertEquals(NB_EXPECTED_LINES_IN_PATH_FILE, len(self._db.fetchall()))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
295 self._db.dropTable("%s_TE_BLRn_path" % (self._project))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
296
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
297
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
298 def _createPathFiles( self, bankFull ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
299 bank = os.path.split(bankFull)[1]
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
300 fileName = "batch_1.fa_BLRn_" + bank + ".align.clean_match.path"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
301 f = open(fileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
302 f.write("1\tQueryName1\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
303 f.write("2\tQueryName1\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
304 f.write("3\tQueryName2\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
305 f.write("4\tQueryName3\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
306 f.write("5\tQueryName1\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
307 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
308 fileName = "batch_2.fa_BLRn_" + bank + ".align.clean_match.path"
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
309 f = open(fileName, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
310 f.write("1\tQueryName4\t2\t250\tsubjectName1\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
311 f.write("2\tQueryName4\t255\t550\tsubjectName2\t5\t255\t0.0002\t32\t78.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
312 f.write("3\tQueryName5\t1\t150\tsubjectName1\t250\t400\t5.1e-59\t132\t98\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
313 f.write("4\tQueryName6\t2\t250\tsubjectName3\t5\t255\t4.1e-39\t132\t88.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
314 f.write("5\tQueryName7\t300\t450\tsubjectName1\t300\t450\t4.1e-39\t132\t80.2\n")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
315 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
316
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
317
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
318 def _readPathResultsFromTableAndFillList( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
319 tablePathAdaptatorInstance = TablePathAdaptator (self._db, "%s_TE_BLRn_path" % (self._project))
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
320 pathList = tablePathAdaptatorInstance.getListOfAllPaths()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
321 return pathList
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
322
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
323
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
324 def _readPathResultsFromFileAndFillList( self ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
325 pathInstance = Path()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
326 pathList = []
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
327 f = open( "%s_BLRn_%s.align.clean_match.path" % (self._project, self._bank) , "r")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
328 while pathInstance.read( f ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
329 pathList.append(pathInstance)
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
330 pathInstance = Path()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
331 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
332 return pathList
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
333
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
334
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
335 def _createFile( self, nameFile ):
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
336 f = open(nameFile, "w")
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
337 f.close()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
338
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
339
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
340 test_suite = unittest.TestSuite()
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
341 test_suite.addTest( unittest.makeSuite( Test_RepbaseBLRnForClassifierStep1 ) )
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
342 if __name__ == "__main__":
94ab73e8a190 Uploaded
m-zytnicki
parents:
diff changeset
343 unittest.TextTestRunner(verbosity=2).run( test_suite )