SEQanswers (
-   Bioinformatics (
-   -   average coverage in target regions (

cmccabe 06-18-2015 06:54 AM

average coverage in target regions
I am trying to get the _X coverage in a bam file for target regions in a bed.


find *bam | parallel 'bedtools coverage -hist -abam IonXpress_009_150603.bam -b epilepsy70_medex.bed | grep ^all > output.bam.hist.txt' > all.txt
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 :).

Maybe this:

coverageBed -abam IonXpress_009_150603.bam -b epilepsy70_medex.bed -counts > count.txt

All times are GMT -8. The time now is 07:45 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.