List of Exons per gene bad covered
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

Login before adding your answer.

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

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

Powered by the version 2.3.6