diff --git a/workflow/rna-seq.nf b/workflow/rna-seq.nf index f36c6e36131912d924e0febb31bb2cb1fe77d84c..46abbf1eee0aadaeb764305ef7a39a0de3ca1bb6 100755 --- a/workflow/rna-seq.nf +++ b/workflow/rna-seq.nf @@ -144,7 +144,7 @@ process parseMetadata { echo "LOG: strandedness metadata parsed: \${stranded}" >>${repRID}.parseMetadata.err # Get spike-in metadata - spike=\$(python3 ${script_parseMeta} -r ${repRID} -m "${experimentSettingsMeta}" -p spike) + spike=\$(python3 ${script_parseMeta} -r ${repRID} -m "${experimentSettingsMeta}" -p spike -e \${endsManual}) echo "LOG: spike-in metadata parsed: \${spike}" >>${repRID}.parseMetadata.err # Get species metadata diff --git a/workflow/scripts/parseMeta.py b/workflow/scripts/parseMeta.py index 3667387fbfe89ad88425d2dcbc163b499786f2c7..eaa2656f9266a5c3de07fb0fc80ad5aeb2ea8422 100644 --- a/workflow/scripts/parseMeta.py +++ b/workflow/scripts/parseMeta.py @@ -10,6 +10,7 @@ def get_args(): parser.add_argument('-r', '--repRID',help="The replicate RID.",required=True) parser.add_argument('-m', '--metaFile',help="The metadata file to extract.",required=True) parser.add_argument('-p', '--parameter',help="The parameter to extract.",required=True) + parser.add_argument('-e', '--endsManual',help="The endness.",required=False) args = parser.parse_args() return args @@ -55,9 +56,9 @@ def main(): # Get strandedness metadata from 'Experiment Settings.csv' if (args.parameter == "stranded"): if (metaFile.Has_Strand_Specific_Information.unique() == "yes"): - if (len(metaFile[metaFile["File_Type"] == "FastQ"]) == 1): + if (args.endsManual=="se"): stranded = "--rna-strandness F" - elif (len(metaFile[metaFile["File_Type"] == "FastQ"]) == 2): + elif (args.endsManual=="pe"): stranded = "--rna-strandness FR" elif (metaFile.Has_Strand_Specific_Information.unique() == "no"): stranded = ""