diff --git a/README.md b/README.md index 0d33c7c..8455810 100644 --- a/README.md +++ b/README.md @@ -94,7 +94,7 @@ for LSTrAP can be generated. ## Adapting LSTrAP to other cluster managers -LSTrAP is designed and tested on an Oracle Grid Engine, though with minimal effort it can be adopted to run on PBS and Torque +LSTrAP is designed and tested on a cluster running the Oracle Grid Engine, though with minimal effort it can be adopted to run on PBS and Torque based systems (and likely others). First, in the configuration file, check the qsub parameters (e.g. jobs that require multiple CPUs to run *-pe cores 4*), that differ between systems are set up properly (the nodes and cores on Torque and PBS need to be set using *-l nodes=4:ppn=2* to request 4 nodes with 2 processes per node). diff --git a/helper/htseq_count_stats.py b/helper/htseq_count_stats.py index 52eb5d4..3071102 100644 --- a/helper/htseq_count_stats.py +++ b/helper/htseq_count_stats.py @@ -23,7 +23,7 @@ mapped_reads = 0 for line in fin: - gene, value = line.strip().split() + gene, value = line.strip().rsplit(maxsplit=1) if gene not in quality_fields: mapped_reads += int(value) diff --git a/pipeline/check/quality.py b/pipeline/check/quality.py index d5422a6..3f6eff8 100644 --- a/pipeline/check/quality.py +++ b/pipeline/check/quality.py @@ -46,7 +46,7 @@ def check_htseq(filename, cutoff=0, log=None): with open(filename) as fin: mapped_reads = 0 for line in fin: - gene, value = line.strip().split() + gene, value = line.strip().rsplit(maxsplit=1) if gene not in quality_fields: mapped_reads += int(value)