yat
0.17.3pre
|
Filter bam reads. More...
#include <yat/omic/BamReadFilter.h>
Public Member Functions | |
BamReadFilter (unsigned char mapQ, unsigned int filtered_flags=BAM_DEF_MASK, unsigned int required_flags=0) | |
bool | operator() (const BamRead &bam) const |
Filter bam reads.
Class is specifically designed to be used with boost::filter_iterator. The code below for example
iterates over all read with mapping quality 20 or greater.
theplu::yat::omic::BamReadFilter::BamReadFilter | ( | unsigned char | mapQ, |
unsigned int | filtered_flags = BAM_DEF_MASK , |
||
unsigned int | required_flags = 0 |
||
) |
Create a functor which accepts reads with mapping quality mapQ (or greater), flag must be unset in bits described by a filter_flags, and flag must be set in bits described by required_flags.
bool theplu::yat::omic::BamReadFilter::operator() | ( | const BamRead & | bam | ) | const |
return true if 1) bam mapping quality is at least mapQ AND 2) bam flags & filter_flags == 0 (no filter bit set) AND 3) bam flags & required_flags == required_flags (all req bits set)