Whole Exome From Whole Genome

Question

[Kasthuri]

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

Answer 1

[Steve Lianoglou]

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.

[DG]

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

[Swbarnes2]

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.

[Sean Davis]

.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.

[Kasthuri]

Thank you very much for the inputs!

Answer 2

[Malachi Griffith]

Another option is the tool: RefCov