git-svn-id: http://svn.tuebingen.mpg.de/ag-raetsch/projects/QPalma@8626 e1793c9e...
authorraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 20:56:40 +0000 (20:56 +0000)
committerraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 20:56:40 +0000 (20:56 +0000)
scripts/PipelineHeuristic.py

index e9985e7..1af7974 100644 (file)
@@ -353,16 +353,16 @@ class PipelineHeuristic:
          exons[1,1]     = acc_pos+1+(self.read_size-don_pos)
 
          _dna = dna[:int(exons[1,1])]
-         _dna = _dna[:exons[1,0]] + orig_seq[don_pos:]
+         _dna = _dna[:exons[1,0]] + est[don_pos:]# only correct if there are no indels!!!
 
-         _currentAcc = currentAcc[:int(exons[1,1])]
-         _currentAcc = [0.25]*len(_currentAcc) 
+         #_currentAcc = currentAcc[:int(exons[1,1])]
+         _currentAcc = [0.25]*len(currentAcc) 
 
-         _currentDon = currentDon[:int(exons[1,1])]
-         _currentDon = [0.25]*len(_currentDon)
+         #_currentDon = currentDon[:int(exons[1,1])]
+         #_currentDon = [0.25]*len(_currentDon)
 
          currentVMatchAlignment = _dna, exons, est, original_est, quality,\
-         _currentAcc, _currentDon
+         _currentAcc, currentDon
 
          stop = cpu()
          self.array_stuff += stop - start