my $genGenomeBamF = 0;
my $sampling = 0;
my $calcCI = 0;
+my $var_opt = 0; # temporarily, only for internal use
my $quiet = 0;
my $help = 0;
"no-bam-output" => sub { $genBamF = 0; },
"output-genome-bam" => \$genGenomeBamF,
"sampling-for-bam" => \$sampling,
+ "var" => \$var_opt,
"calc-ci" => \$calcCI,
"ci-memory=i" => \$NMB,
"time" => \$mTime,
$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);
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"; }
if ($mTime) { $time_start = time(); }
-if ($calcCI) {
+if ($calcCI || $var_opt) {
$command = $dir."rsem-run-gibbs $refName $sampleName $sampleToken $BURNIN $NCV $SAMPLEGAP";
$command .= " -p $nThreads";
+ if ($var_opt) { $command .= " --var"; }
if ($quiet) { $command .= " -q"; }
&runCommand($command);
+}
+if ($calcCI) {
system("mv $sampleName.isoforms.results $imdName.isoforms.results.bak1");
system("mv $sampleName.genes.results $imdName.genes.results.bak1");
&collectResults("$imdName.iso_res", "$sampleName.isoforms.results"); # isoform level