|
|
@@ -47,7 +47,7 @@ task ReadStats { |
|
|
|
cat ${in_sam_align_RNA} | grep -v '@' | awk '($2!=4)' | cut -f 3 | sed 's/.*;//g' | awk '{a[$1]++}END{for(i in a){printf "%s\t%d\n",i,a[i]}}' > $groupedReadCount |
|
|
|
|
|
|
|
|
|
|
|
if [ $(cat $groupedReadCount | grep '^mRNA' | wc -l ) -gt 0] |
|
|
|
if [ $(cat $groupedReadCount | grep '^mRNA' | wc -l ) -gt 0 ] |
|
|
|
then |
|
|
|
n_mRNA=$(cat $groupedReadCount | grep '^mRNA' | cut -f 2 ) |
|
|
|
else |
|
|
@@ -75,7 +75,7 @@ task ReadStats { |
|
|
|
n_YRNA=0 |
|
|
|
fi |
|
|
|
|
|
|
|
if [$(cat $groupedReadCount | grep -E '^misc_RNA|small|guide_RNA|vault_RNA' | wc -l ) -gt 0 ] |
|
|
|
if [ $(cat $groupedReadCount | grep -E '^misc_RNA|small|guide_RNA|vault_RNA' | wc -l ) -gt 0 ] |
|
|
|
then |
|
|
|
n_otsmall=$(cat $groupedReadCount | grep -E '^misc_RNA|small|guide_RNA|vault_RNA' | cut -f 2 | awk '{sum+=$1}END{print sum}') |
|
|
|
else |