More posts like this
-
sambamba-depth: All files must be coordinate-sorted error?
written by MobiusT •
-
Samtools index and sambamba-depth error?
written by MobiusT •
-
What is read depth exactly?
written by MobiusT •
-
Find mean read length of a bam file dataset of reads?
written by MobiusT •
-
differences in covered regions in exome and whole genome sequencing
written by miaowzai
-
Statistical analysis of a whole exome VCF file
written by adam.waring •
-
Input file for MDR(Multifactoral dimension reduction) for exome data
written by vvividhya93 •
-
predicting Alternative splicing
written by rajesh
-
Merging multiplexed bam files
written by Nikleotide
-
Analysing Multiple Exome Datasets
written by User 9126 •
Do you mean that you have a bam file generated from WGS data and that you want to extract the exonic alignments? If so take a look into intersectBed.
While possible the solution from Pierre Lindenbaum will be much faster as (if the bam is sorted and indexed) samtools will use random access to pull the entries directly whereas intersect first needs to find them manually.
This would be a theoretically equivalent dataset not a real exome dataset (since that generally uses capture techology).