Question: how to calculate the coverage over a specif coordinate
1
gravatar for Lila M
3 months ago by
Lila M 370
UK
Lila M 370 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 • 341 views
ADD COMMENTlink written 3 months ago by Lila M 370
1
gravatar for Medhat
3 months ago by
Medhat7.1k
Poland
Medhat7.1k wrote:

you should use -abam

for eexample:

bedtools coverage -abam sorted.bam -b regions.bed.gz | grep ^all > mycoverage.tab
ADD COMMENTlink modified 3 months ago • written 3 months ago by Medhat7.1k
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 3 months ago • written 3 months ago by Lila M 370
1

I think they both do the same job

ADD REPLYlink modified 3 months ago • written 3 months ago by Medhat7.1k
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: 1365 users visited in the last hour