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

index 3a19c03..028c248 100644 (file)
@@ -336,7 +336,7 @@ class PipelineHeuristic:
       qualityPlifs = self.qualityPlifs
       # inlined
 
-      IntronScore = calculatePlif(h, [90])[0]
+      IntronScore = calculatePlif(h, [self.intron_size])[0]
       dummyAcceptorScore = calculatePlif(a, [0.25])[0] 
       dummyDonorScore = calculatePlif(d, [0.25])[0]
       print IntronScore,dummyAcceptorScore,dummyDonorScore
@@ -345,6 +345,8 @@ class PipelineHeuristic:
       for (don_pos,don_score) in alt_don:
          start = cpu()
 
+         acc_pos = don_pos + self.intron_size
+         
          exons[0,0]     = 0
          exons[0,1]     = don_pos
          exons[1,0]     = acc_pos+1