This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Whole Exome From Whole Genome

Is there a way I can extract whole exome portions from the whole genome data using any standard tools? I need the whole exome in the bam format.

Thanks -Kasthuri

exome genome

2 answers

If you define your "exome regions" in BED format, you can use something like intersectBed from bedtools to extract the reads from whole genome BAM file.

That was my thought too, but do you know of a place where such a file can be found for a whole genome, like mouse? We do exon capture only, and I've been using the probe .beds, but they are bigger than exons.

One of the example tutorials on the galaxy server is extracting exonic regions only from the human genome.

.bed files can typically be created from UCSC using the Table Browser. For mouse, you could use the Table Browser to get the regions for any number of transcript models in .bed format.

Thank you very much for the inputs!

Another option is the tool: RefCov

Log in to answer this question.