List of Exons per gene bad covered
0
0
Entering edit mode
7.0 years ago
windsur ▴ 20

Dear all, this is my first post. I 'm new on NGS and I have some questions that I hope you can guide me.

I want to create a python script to calculate the coverage. But I would like an output with headers as: sampleID, gene, chr, exon, start, end, and avg cov.

For that I tried with bedtools

bedtools coverage [options] -a<BAM> -b<BED> > output.txt

With my bam file I converted in bed file using bamtobed, but I think I should do something before or while I create the bedfile.

Any idea that can help me? thanks a lot!

exon coverage python • 1.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 2944 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6