Mercurial > repos > rlegendre > ribo_tools

diff metagene_readthrough.py @ 14:344bacf6acb8 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author rlegendre
date Fri, 10 Apr 2015 03:18:56 -0400
parents 7c944fd9907e
children 702e60e819c2
line wrap: on
line diff
--- a/metagene_readthrough.py	Thu Apr 09 11:35:48 2015 -0400
+++ b/metagene_readthrough.py	Fri Apr 10 03:18:56 2015 -0400
@@ -121,7 +121,7 @@
 def check_homo_coverage(gene,GFF,start,stop, aln) :
     
     chrom = GFF[gene]['chrom']
-    ## compute nunber of nucleotides in CDS with a coverage equal to zero
+    ## compute number of nucleotides in CDS with a coverage equal to zero
     nt_cds_cov = 0
     nt_cds_num = 0
     for i in range(1,GFF[gene]['exon_number']+1) :
@@ -131,7 +131,10 @@
                 nt_cds_cov += 1
                     
     ## compute percent of CDS no covering
-    percent = nt_cds_cov*100/nt_cds_num
+    if nt_cds_cov == 0 or nt_cds_num == 0:
+        percent = 0
+    else:
+        percent = nt_cds_cov*100/nt_cds_num
     
     ## compute number of nucleotides with no coverage in extension
     nt_no_cover = 0