Skip to content
Snippets Groups Projects
Commit 1902055c authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Pass endness to stranded param

parent 471380f7
Branches
Tags
2 merge requests!37v0.0.1,!16Resolve "process_getRef"
Pipeline #5836 failed with stages
in 8 minutes
......@@ -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
......
......@@ -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 = ""
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment