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

index 6221ee6..1d00f38 100644 (file)
@@ -315,6 +315,9 @@ class PipelineHeuristic:
       dna   = currentDNASeq
 
       alt_don,alt_acc = self.findHighestScoringSpliceSites(currentAcc,currentDon)
+      print alt_don
+      print alt_acc
+      
       alternativeScores = []
       
       # inlined
@@ -417,9 +420,10 @@ class PipelineHeuristic:
          print original_est,original_est_cut
          
          AcceptorScore = calculatePlif(d, [acc_score])[0]
-         print AcceptorScore,acc_score
+         print AcceptorScore,acc_score,acc_pos
          if acc_score<0.1:
              print currentAcc[0:50]
+             print currentDon[0:50]
          
          score = computeSpliceAlignScoreWithQuality(original_est_cut, quality, qualityPlifs, run, self.currentPhi)
          score += AcceptorScore + IntronScore + dummyDonorScore