Skip to content

Commit

Permalink
Fix dedup/${id}_peaks.narrowPeak filename
Browse files Browse the repository at this point in the history
  • Loading branch information
donald committed Sep 13, 2022
1 parent aab07b4 commit e1cc50f
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions Adapting_new_script_for_computing_cluster.sh
Original file line number Diff line number Diff line change
Expand Up @@ -149,9 +149,9 @@ execute() {
mapFrac=$(printf "%.4f" $(echo $map / 100 | bc -l))
dedup=$(printf "%.4f" $(grep -Po 'Unknown Library.*' dedup/${id}_dedupMetric.txt | awk '{print $10}'))
filtered=$(expr $(samtools view -c dedup/${id}_dedup.bam) / 2)
ripReads=$(expr $(bedtools intersect -a dedup/${id}_dedup.bam -b dedup_peaks/${id}_dedup_peaks.narrowPeak -u -ubam | samtools view -c) / 2)
ripReads=$(expr $(bedtools intersect -a dedup/${id}_dedup.bam -b dedup_peaks/${id}_peaks.narrowPeak -u -ubam | samtools view -c) / 2)
FRiP=$(printf "%.4f" $(echo $fripReads / $filtered | bc -l))
peaks=$(expr $(wc -l <"dedup_peaks/${id}_dedup_peaks.narrowPeak") - 1)
peaks=$(expr $(wc -l <"dedup_peaks/${id}_peaks.narrowPeak") - 1)
echo -e "$id\t$fragments\t$filtered\t$mapFrac\t$dedup\t$FRiP\t$peaks" >> cnt_gata_qc_metrics.txt

# Here I have no idea how to transfer this
Expand Down

0 comments on commit e1cc50f

Please sign in to comment.