SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   average coverage in target regions (http://seqanswers.com/forums/showthread.php?t=60547)

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.

Code:

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 :).

EDIT:
Maybe this:
Code:

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.