Whole Exome From Whole Genome
2
6
Entering edit mode
12.8 years ago
Kasthuri ▴ 300

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 • 5.3k views
ADD COMMENT
4
Entering edit mode
12.8 years ago

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.

ADD COMMENT
1
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

Thank you very much for the inputs!

ADD REPLY
0
Entering edit mode
12.8 years ago

Another option is the tool: RefCov

ADD COMMENT

Login before adding your answer.

Traffic: 1452 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