X-Git-Url: https://git.donarmstrong.com/?p=rsem.git;a=blobdiff_plain;f=rsem-prepare-reference;h=56d6972d12ad99f3cecceda21c6dfbc6133feebd;hp=78743e946f2b2f1b6b02ecdb6cbeeaaebdf3c70e;hb=683863b75f8d8bef2461039a6911b0e9619cc113;hpb=fc69cf6af24c0550e55447fc82f01cb6f90c1c42 diff --git a/rsem-prepare-reference b/rsem-prepare-reference index 78743e9..56d6972 100755 --- a/rsem-prepare-reference +++ b/rsem-prepare-reference @@ -97,12 +97,12 @@ if (!$no_bowtie) { sub runCommand { print $_[0]."\n"; my $status = system($_[0]); - if ($status != 0) { - my $errmsg; - if (scalar(@_) > 1) { $errmsg = $_[1]; } - else { $errmsg = "\"$command\" failed! Plase check if you provide correct parameters/options for the pipeline!"; } - print $errmsg."\n"; - exit(-1); + if ($status != 0) { + my $errmsg = ""; + if (scalar(@_) > 1) { $errmsg .= $_[1]."\n"; } + $errmsg .= "\"$_[0]\" failed! Plase check if you provide correct parameters/options for the pipeline!\n"; + print $errmsg; + exit(-1); } print "\n"; } @@ -155,7 +155,7 @@ Each line of should be of the form: gene_id transcript_id with the two fields separated by a tab character. - + If you are using a GTF file for the "UCSC Genes" gene set from the UCSC Genome Browser, then the "knownIsoforms.txt" file (obtained from the "Downloads" section of the UCSC Genome Browser site) is of this format. If this option is off, then the mapping of isoforms to genes depends on whether the --gtf option is specified. If --gtf is specified, then RSEM uses the "gene_id" and "transcript_id" attributes in the GTF file. Otherwise, RSEM assumes that each sequence in the reference sequence files is a separate gene.