git-svn-id: http://svn.tuebingen.mpg.de/ag-raetsch/projects/QPalma@8598 e1793c9e...
authorraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 19:27:34 +0000 (19:27 +0000)
committerraetsch <raetsch@e1793c9e-67f9-0310-80fc-b846ff1f7b36>
Thu, 17 Apr 2008 19:27:34 +0000 (19:27 +0000)
qpalma/computeSpliceAlignWithQuality.py

index 598cbc5..cd137c0 100644 (file)
@@ -235,7 +235,8 @@ def  computeSpliceAlignScoreWithQuality(dna, exons, est, original_est, quality,
                    score += numqP[(estletter-1)*6*numQualSuppPoints + dnaletter*numQualSuppPoints+0] ;
                    
                elif Lower == len(cur_plf.limits):
-                   score += numqP[(estletter-1)*6*numQualSuppPoints + dnaletter*numQualSuppPoints+numQualSuppPoints-1] ;
+                   print len(numqP),estletter,numQualSuppPoints,dnaletter,(estletter-1)*6*numQualSuppPoints + dnaletter*numQualSuppPoints + numQualSuppPoints-1
+                   score += numqP[(estletter-1)*6*numQualSuppPoints + dnaletter*numQualSuppPoints + numQualSuppPoints-1] ;
 
                else:
                    # because we count from 0 in python