diff --git a/build/lib/liqa_src/PreProcess_gtf.pl b/build/lib/liqa_src/PreProcess_gtf.pl
index f9431de..6b26b80 100644
--- a/build/lib/liqa_src/PreProcess_gtf.pl
+++ b/build/lib/liqa_src/PreProcess_gtf.pl
@@ -83,6 +83,9 @@
     }
 }
 
+my $number_genes = scalar keys %knownGene;
+print "$number_genes genes are being processed...\n";
+
 #open(REF, $ref
 open (RRR, ">$output");
 
diff --git a/build/lib/liqa_src/quantify.py b/build/lib/liqa_src/quantify.py
index 63401e5..507a16c 100644
--- a/build/lib/liqa_src/quantify.py
+++ b/build/lib/liqa_src/quantify.py
@@ -141,7 +141,7 @@ def auto_dict():
                 tmpisolength[iii] = tmpisolength[iii]/2566
             #print(gene+"\t"+str(tmpisolength[iii]))
     else:
-        genelength = 1
+        genelength = 100
         tmpisolength = tmpisoinf[0].split(",")
         for iii in range(len(tmpisolength)-1):
             tmpisolength[iii] = 0
@@ -525,11 +525,13 @@ def auto_dict():
 
     print(gene+"\t"+str(iterCount)+" iterations\tDone!")
 
-    rpg_lengthcorrected = readCount/genelength*100
+    #rpg_lengthcorrected = readCount/genelength*100
     for i in range(len(Alpha)):
         #isoformRelativeAbundances[i] = (Alpha[i]/isoformLength[isoformNames[i]]+tmpisolength[i]*weightF) /(sumTheta)
         isoformRelativeAbundances[i] = (Alpha[i] + tmpisolength[i]*weightF) /(sumTheta)
         #print(gene+"\t"+str(geneCount)+"\t"+str(iterCount)+"\t"+isoformNames[i]+"\t"+str(isoformRelativeAbundances[i])+"\t"+str(tmpTime))
+        
+        rpg_lengthcorrected = readCount/genelength*100*isoformRelativeAbundances[i]
 
         #OUT.write(gene+"\t"+isoformNames[i]+"\t"+str(readCount)+"\t"+str(isoformRelativeAbundances[i])+"\t"+str(rpg_lengthcorrected)+"\n") ## write results into specified file
         OUT.write(gene+"\t"+isoformNames[i]+"\t"+str(isoformRelativeAbundances[i])+"\t"+str(rpg_lengthcorrected)+"\n")
diff --git a/build/scripts-3.7/PreProcess_gtf.pl b/build/scripts-3.7/PreProcess_gtf.pl
index f9431de..6b26b80 100755
--- a/build/scripts-3.7/PreProcess_gtf.pl
+++ b/build/scripts-3.7/PreProcess_gtf.pl
@@ -83,6 +83,9 @@
     }
 }
 
+my $number_genes = scalar keys %knownGene;
+print "$number_genes genes are being processed...\n";
+
 #open(REF, $ref
 open (RRR, ">$output");
 
diff --git a/dist/liqa-1.1.15-py3-none-any.whl b/dist/liqa-1.1.15-py3-none-any.whl
new file mode 100644
index 0000000..4c76b0c
Binary files /dev/null and b/dist/liqa-1.1.15-py3-none-any.whl differ
diff --git a/dist/liqa-1.1.15.tar.gz b/dist/liqa-1.1.15.tar.gz
new file mode 100644
index 0000000..152755b
Binary files /dev/null and b/dist/liqa-1.1.15.tar.gz differ
diff --git a/liqa.egg-info/PKG-INFO b/liqa.egg-info/PKG-INFO
index 73eacda..2d00aac 100644
--- a/liqa.egg-info/PKG-INFO
+++ b/liqa.egg-info/PKG-INFO
@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: liqa
-Version: 1.1.13
+Version: 1.1.15
 Summary: A statistical tool to quantify isoform-specific expression using long-read RNA-seq
 Home-page: https://github.com/WGLab/LIQA
 Author: Yu Hu
diff --git a/liqa_src/quantify.py b/liqa_src/quantify.py
index 507a16c..a641725 100644
--- a/liqa_src/quantify.py
+++ b/liqa_src/quantify.py
@@ -108,7 +108,7 @@ def auto_dict():
 startTime = time.time()
 
 #OUT.write("GeneName\tIsoformName\tNumberOfReads\tRelativeAbundance\n") ## Header of Results
-OUT.write("GeneName\tIsoformName\tRelativeAbundance\tReadPerGene_corrected\n")
+OUT.write("GeneName\tIsoformName\tReadPerGene_corrected\tRelativeAbundance\n")
 
 for gene in geneStructureInformation:
 
@@ -534,7 +534,7 @@ def auto_dict():
         rpg_lengthcorrected = readCount/genelength*100*isoformRelativeAbundances[i]
 
         #OUT.write(gene+"\t"+isoformNames[i]+"\t"+str(readCount)+"\t"+str(isoformRelativeAbundances[i])+"\t"+str(rpg_lengthcorrected)+"\n") ## write results into specified file
-        OUT.write(gene+"\t"+isoformNames[i]+"\t"+str(isoformRelativeAbundances[i])+"\t"+str(rpg_lengthcorrected)+"\n")
+        OUT.write(gene+"\t"+isoformNames[i]+"\t"+str(rpg_lengthcorrected)+"\t"+str(isoformRelativeAbundances[i])+"\n")
 
 OUT.close()
             
diff --git a/setup.py b/setup.py
index 06f2f28..8c97cac 100644
--- a/setup.py
+++ b/setup.py
@@ -5,7 +5,7 @@
 
 setuptools.setup(
     name="liqa",
-    version="1.1.13",
+    version="1.1.15",
     author="Yu Hu",
     author_email="huyu999999@gmail.com",
     description="A statistical tool to quantify isoform-specific expression using long-read RNA-seq",