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

Fix echo greps

parent 821253a0
Branches
Tags
2 merge requests!37v0.0.1,!36Metadata output update
......@@ -107,7 +107,7 @@ countData:
script:
- singularity run 'docker://bicf/subread2:2.0.0' featureCounts -T 20 -a /project/BICF/BICF_Core/shared/gudmap/references/GRCh38.p12.v31/genome.gtf -G /project/BICF/BICF_Core/shared/gudmap/references/GRCh38.p12.v31/genome.fna -g 'gene_name' -o Q-Y5F6_1M.se.featureCounts -s 1 -R SAM --primary --ignoreDup ./test_data/bam/small/Q-Y5F6_1M.se.sorted.deduped.bam
- singularity run 'docker://bicf/subread2:2.0.0' Rscript ./workflow/scripts/calculateTPM.R --count ./test_data/counts/small/Q-Y5F6_1M.se.featureCounts
- assignedReads=grep -m 1 'Assigned' *.countData.summary | grep -oe '\([0-9.]*\)'
- assignedReads=$(grep -m 1 'Assigned' *.countData.summary | grep -oe '\([0-9.]*\)')
- pytest -m makeFeatureCounts
makeBigWig:
......
......@@ -877,7 +877,7 @@ process countData {
echo -e "LOG: counted" >> ${repRID}.countData.log
# extract assigned reads
assignedReads=grep -m 1 'Assigned' *.countData.summary | grep -oe '\\([0-9.]*\\)'
assignedReads=\$(grep -m 1 'Assigned' *.countData.summary | grep -oe '\\([0-9.]*\\)')
echo -e \${assignedReads} > assignedReads.csv
echo -e "LOG: assigned reads: "\${assignedReads} >> ${repRID}.countData.log
......
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