RSEM v1.1.8

diff --git a/README.md b/README.md
index 2bb60152c0492f6b8260169a72edfe877b5938c0..f8c4e2709436f97e446ac3dbbda2306818df70f4 100644 (file)
--- a/README.md
+++ b/README.md
@@ -148,7 +148,8 @@ The plots generated depends on read type and user configuration. It
 may include fragment length distribution, mate length distribution,
 read start position distribution (RSPD), quality score vs observed
 quality given a reference base, position vs percentage of sequencing
-error given a reference base and histogram of read alignments.
+error given a reference base and histogram of reads with different
+number of alignments.
 
 fragment length distribution and mate length distribution: x-axis is fragment/mate length, y axis is the probability of generating a fragment/mate with the associated length
 
@@ -158,7 +159,7 @@ Quality score vs. observed quality given a reference base: x-axis is Phred quali
 
 Position vs. percentage sequencing error given a reference base: x-axis is position and y-axis is percentage sequencing error
 
-Histogram of read alignments: x-axis is the number of alignments a read has and y-axis is the number of such reads
+Histogram of reads with different number of alignments: x-axis is the number of alignments a read has and y-axis is the number of such reads. The inf in x-axis means number of reads filtered due to too many alignments
  
 ## <a name="example"></a> Example
 

diff --git a/rsem-plot-model b/rsem-plot-model
index 8ff626a530db1808ac4067a91ac04297df423d76..7466eea2398be524e3f073b1ac7d914b3c6a8d83 100755 (executable)
--- a/rsem-plot-model
+++ b/rsem-plot-model
@@ -124,6 +124,6 @@ if (model_type == 1 || model_type == 3) {
 close(con)
 
 pair <- read.table(file = cntF, skip = 3, sep = "\t")
-barplot(pair[,2], names.arg = pair[,1], xlab = "Number of Alignments", ylab = "Number of Reads", main = "Histogram of Read Alignments")
+barplot(pair[,2], names.arg = pair[,1], xlab = "Number of Alignments", ylab = "Number of Reads", main = "Histogram of Reads with Different Number of  Alignments")
 
 dev.off()