X-Git-Url: https://git.donarmstrong.com/?p=rsem.git;a=blobdiff_plain;f=rsem-calculate-expression;h=c01c70ab2e7067f2a4b2f3683215fd87314645e1;hp=5614ecda945f8fe2b261fff8e3c9ed5bc891979d;hb=23f7edf618b2cdcd08cf6d78092350ae69817eb0;hpb=237bbdf363c9e42ee24e2fd63106dccf20d9bf2f diff --git a/rsem-calculate-expression b/rsem-calculate-expression index 5614ecd..c01c70a 100755 --- a/rsem-calculate-expression +++ b/rsem-calculate-expression @@ -173,7 +173,7 @@ if (!(-d $stat_dir) && !mkdir($stat_dir)) { print "Fail to create folder $stat_d $imdName = "$temp_dir/$sampleToken"; -if (!$is_sam && !$is_bam && $phred33 + $phred64 + $solexa == 0) { $phred33 = 1; } +if (!$is_sam && !$is_bam && !$no_qual && ($phred33 + $phred64 + $solexa == 0)) { $phred33 = 1; } my ($mate_minL, $mate_maxL) = (1, $maxL); @@ -184,13 +184,12 @@ my $command = ""; if (!$is_sam && !$is_bam) { $command = $bowtie_path."bowtie"; - if ($read_type == 0 || $read_type == 2) { $command .= " -f"; } + if ($no_qual) { $command .= " -f"; } else { $command .= " -q"; } if ($phred33) { $command .= " --phred33-quals"; } elsif ($phred64) { $command .= " --phred64-quals"; } elsif ($solexa) { $command .= " --solexa-quals"; } - else { print "Oh, no!!!"; exit(2); } $command .= " -n $C -e $E -l $L"; if ($read_type == 2 || $read_type == 3) { $command .= " -I $minL -X $maxL"; }