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

Update filenames to new convention

parent b3a65d91
2 merge requests!58Develop,!53Resolve "process_derivaUpload"
...@@ -988,7 +988,7 @@ process makeBigWig { ...@@ -988,7 +988,7 @@ process makeBigWig {
tuple path (bam), path (bai) from dedupBam_makeBigWig tuple path (bam), path (bai) from dedupBam_makeBigWig
output: output:
path ("${repRID}.bw") into bigwig path ("${repRID}_sorted.deduped.bw") into bigwig
script: script:
""" """
...@@ -1019,7 +1019,7 @@ process countData { ...@@ -1019,7 +1019,7 @@ process countData {
output: output:
path ("*_tpmTable.csv") into counts path ("*_tpmTable.csv") into counts
path ("*.countData.summary") into countsQC path ("*_countData.summary") into countsQC
path ("assignedReads.csv") into assignedReadsInfer_fl path ("assignedReads.csv") into assignedReadsInfer_fl
script: script:
...@@ -1055,11 +1055,11 @@ process countData { ...@@ -1055,11 +1055,11 @@ process countData {
echo -e "LOG: counted" >> ${repRID}.countData.log echo -e "LOG: counted" >> ${repRID}.countData.log
# extract assigned reads # extract assigned reads
grep -m 1 'Assigned' *.countData.summary | grep -oe '\\([0-9.]*\\)' > assignedReads.csv grep -m 1 'Assigned' *_countData.summary | grep -oe '\\([0-9.]*\\)' > assignedReads.csv
# calculate TPM from the resulting countData table # calculate TPM from the resulting countData table
echo -e "LOG: calculating TPM with R" >> ${repRID}.countData.log echo -e "LOG: calculating TPM with R" >> ${repRID}.countData.log
Rscript calculateTPM.R --count "${repRID}.countData" Rscript calculateTPM.R --count "${repRID}_countData"
# convert gene symbols to Entrez id's # convert gene symbols to Entrez id's
echo -e "LOG: convert gene symbols to Entrez id's" >> ${repRID}.countData.log echo -e "LOG: convert gene symbols to Entrez id's" >> ${repRID}.countData.log
...@@ -1463,7 +1463,7 @@ process uploadQC { ...@@ -1463,7 +1463,7 @@ process uploadQC {
exist=\$(curl -s https://${source}/ermrest/catalog/2/entity/RNASeq:mRNA_QC/Replicate=${repRID}) exist=\$(curl -s https://${source}/ermrest/catalog/2/entity/RNASeq:mRNA_QC/Replicate=${repRID})
if [ "\${exist}" != "[]" ] if [ "\${exist}" != "[]" ]
then then
rids=\$(echo $exist | grep -o '\\"RID\\":\\".\\{7\\}' | sed 's/^.\\{7\\}//') rids=\$(echo \${exist} | grep -o '\\"RID\\":\\".\\{7\\}' | sed 's/^.\\{7\\}//')
for rid in \${rids} for rid in \${rids}
do do
python3 deleteEntry.py -r \${rid} -t mRNA_QC -o ${source} -c \${cookie} python3 deleteEntry.py -r \${rid} -t mRNA_QC -o ${source} -c \${cookie}
......
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