Skip to content

Commit

Permalink
Update bamtobigwig.sh
Browse files Browse the repository at this point in the history
  • Loading branch information
laurenchoate authored Jan 15, 2019
1 parent b598fe9 commit f5bf5d0
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions dnase/BamToBigWig/bamtobigwig.sh
Original file line number Diff line number Diff line change
Expand Up @@ -223,7 +223,7 @@ if [[ "$SEQ" == "SE" ]] ; then
bedtools genomecov -bg -i ${TMPDIR}/$j.bed.gz -g ${CHINFO} -strand - > ${TMPDIR}/$j\_minus.noinv.bedGraph


## Invert minus strand and shift reads.
## Invert minus strand.
cat ${TMPDIR}/$j\_minus.noinv.bedGraph | awk 'BEGIN{OFS="\t"} {print $1,$2,$3,$4}' > ${TMPDIR}/$j\_minus.bedGraph

## Then to bigWig
Expand All @@ -249,7 +249,7 @@ elif [[ "$SEQ" == "PE" ]] ; then
bedtools genomecov -bg -i ${TMPDIR}/$j.bed.gz -g ${CHINFO} -strand - > ${TMPDIR}/$j\_minus.noinv.bedGraph


## Invert minus strand and shift reads.
## Invert minus strand.
cat ${TMPDIR}/$j\_minus.noinv.bedGraph | awk 'BEGIN{OFS="\t"} {print $1,$2,$3,$4}' > ${TMPDIR}/$j\_minus.bedGraph

## Then to bigWig
Expand Down

0 comments on commit f5bf5d0

Please sign in to comment.