I am trying to get the _X coverage in a bam file for target regions in a bed.
The files are too large to attach but after the -abam is the bam file to use, after the -b is the terget bed to use with the baits, after the grep^all is the bedTools coverage output and the > all.txt are the results. Is there something missing? Thank you .
EDIT:
Maybe this:
Code:
find *bam | parallel 'bedtools coverage -hist -abam IonXpress_009_150603.bam -b epilepsy70_medex.bed | grep ^all > output.bam.hist.txt' > all.txt
EDIT:
Maybe this:
Code:
coverageBed -abam IonXpress_009_150603.bam -b epilepsy70_medex.bed -counts > count.txt