Biostar Beta. Not for public use.
Calculating coverage on a list of individual genes in WES data?
0
Entering edit mode
2.5 years ago
gvelz • 0

I have whole exome sequencing VCF and CRAM files from a few human samples and would like to calculate the sequence coverage on a list of individual genes. Is there a good way to go about doing this?

ADD COMMENTlink
0
Entering edit mode

Thanks! I installed Alfred, but am unsure how to make a list of genes in .bed format. Also, does this program work for .cram files?

ADD REPLYlink
0
Entering edit mode

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized.

This comment belongs under @trausch's answer.

ADD REPLYlink
0
Entering edit mode

You should be able to get a BED format gene list following this: How To Get Bed File Containing Exons Of Canonical Transcripts And Their Corresponding Gene Symbols

Choose the genome build you are interested. Rest of the directions should work.

ADD REPLYlink
0
Entering edit mode

I recently added cram support. If you clone the github repository there is an example Rscript that downloads ccds genes from UCSC and stores all exons in a BED file.

ADD REPLYlink
1
Entering edit mode
14 months ago
trausch ♦ 1.3k
Germany

If you want a fragment count table for these genes you can use Alfred:

alfred count_dna -o cov.tsv.gz -i genes.bed <input.bam>

For estimating the on-target rate and the avg. coverage per target region:

alfred qc -b genes.bed -o qc.tsv.gz -r <genome.fa> <input.bam>

Static binaries are available here or via Bioconda.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3