Question: how to calculate the coverage over a specif coordinate
2
gravatar for Lila M
16 months ago by
Lila M 410
UK
Lila M 410 wrote:

Hi everybody, I have bam files for RNAseq, and I know the coordinates for specific exons of interest. I would like to know if there is any way to overlap my bam files over those specific coordinates and get the coverage for those regions. I usually use bedtools coverage, but to compare between samples. I was wondering if I can create a bed files with my coordinates and do:

bedtools coverage -a my.bam -b mycoordintes.bed -s

will that work?

Many thanks!

rna-seq coverage exon • 931 views
ADD COMMENTlink written 16 months ago by Lila M 410
1
gravatar for Medhat
16 months ago by
Medhat8.0k
Poland
Medhat8.0k wrote:

you should use -abam

for eexample:

bedtools coverage -abam sorted.bam -b regions.bed.gz | grep ^all > mycoverage.tab
ADD COMMENTlink modified 16 months ago • written 16 months ago by Medhat8.0k
1

what about?

java -jar GenomeAnalysisTK.jar -T DepthOfCoverage -R [reference fasta] -I [bam file] -L [bed file]

which one is better?

ADD REPLYlink modified 16 months ago • written 16 months ago by Lila M 410
1

I think they both do the same job

ADD REPLYlink modified 16 months ago • written 16 months ago by Medhat8.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1007 users visited in the last hour