Entering edit mode
8.8 years ago
rodrigo.nicodemos82
▴
10
Hello,
I am new to creating bash scripts to run commands. I've tried to run this simple bash script to calculate the genome coverage of a bam file I have using the bedTools genomeCoverageBed command. I am not sure what I am doing wrong, but when I run this script I can a file that contains the printed command.
Here is my script:
#!/bin/bash
echo "This is a Test:"
BAM_FILE=/RNA-seq/TopHat_Output/WT_Rep_1_accepted_hits_sorted.bam
echo genomeCoverageBed -ibam $BAM_FILE -bg > WT_Rep_1_genomecov.bedgraph
Thank you
Rodrigo