From 0889126d07bc92e2eeec30ee218d8bfa2660dad6 Mon Sep 17 00:00:00 2001 From: Bo Li Date: Tue, 26 Jul 2011 01:51:38 -0500 Subject: [PATCH] MAPQ changed --- rsem-calculate-expression | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/rsem-calculate-expression b/rsem-calculate-expression index 93fabce..94764f4 100755 --- a/rsem-calculate-expression +++ b/rsem-calculate-expression @@ -624,7 +624,7 @@ Only generated when --out-bam is specified. genomic coordinates. Alignments of reads that have identical genomic coordinates (i.e., alignments to different isoforms that share the same genomic region) are collapsed into one alignment. The MAPQ field -of each alignment is set to max(100, floor(-10 * log10(1.0 - w) + +of each alignment is set to min(100, floor(-10 * log10(1.0 - w) + 0.5)), where w is the posterior probability of that alignment being the true mapping of a read. In addition, RSEM pads a new tag ZW:f:value, where value is a single precision floating number -- 2.39.2