Filter maf blocks to remove ambiguously aligned or non-informative regions. More...

#include <Bpp/Seq/Io/Maf/AlignmentFilterMafIterator.h>

Inheritance diagram for bpp::AlignmentFilter2MafIterator:

Collaboration diagram for bpp::AlignmentFilter2MafIterator:

Public Member Functions
	AlignmentFilter2MafIterator (MafIterator *iterator, const std::vector< std::string > &species, unsigned int windowSize, unsigned int step, unsigned int maxGap, unsigned int maxPos, bool keepTrashedBlocks, bool missingAsGap)
	Create a new AlignmentFilter2MafIterator with absolute thresholds. More...

	AlignmentFilter2MafIterator (MafIterator *iterator, const std::vector< std::string > &species, unsigned int windowSize, unsigned int step, double maxPropGap, unsigned int maxPos, bool keepTrashedBlocks, bool missingAsGap)
	Create a new AlignmentFilterMafIterator with relative thresholds. More...

MafBlock *	nextRemovedBlock ()
	Get the next available removed alignment block. More...

void	setLogStream (std::shared_ptr< OutputStream > logstream)

void	addIterationListener (IterationListener *listener)

MafBlock *	nextBlock ()
	Get the next available alignment block. More...

bool	isVerbose () const

void	setVerbose (bool yn)

Protected Member Functions
virtual void	fireIterationStartSignal_ ()

virtual void	fireIterationMoveSignal_ (const MafBlock &currentBlock)

virtual void	fireIterationStopSignal_ ()

Protected Attributes
MafIterator *	iterator_

MafBlock *	currentBlock_

std::shared_ptr< OutputStream >	logstream_

std::vector< IterationListener * >	iterationListeners_

bool	started_

bool	verbose_

Private Member Functions
MafBlock *	analyseCurrentBlock_ ()

Private Attributes
std::vector< std::string >	species_

unsigned int	windowSize_

unsigned int	step_

unsigned int	maxGap_

double	maxPropGap_

unsigned int	maxPos_

std::deque< MafBlock * >	blockBuffer_

std::deque< MafBlock * >	trashBuffer_

std::deque< std::vector< bool > >	window_

bool	keepTrashedBlocks_

bool	missingAsGap_

bool	relative_

Detailed Description

Filter maf blocks to remove ambiguously aligned or non-informative regions.

This iterators offers a different algorithm than AlignmentFilterMafIterator. It takes two parameters: g=maxGap and n=maxPos. Windows with more than n positions containing each of them more than g=maxPos gaps will be discarded. In addition, consecutives patterns are only counted once. In case a sequence from the list is missing, it can be either ignored or counted as a full sequence of gaps.

Definition at line 176 of file AlignmentFilterMafIterator.h.

Constructor & Destructor Documentation

◆ AlignmentFilter2MafIterator() [1/2]

bpp::AlignmentFilter2MafIterator::AlignmentFilter2MafIterator	(	MafIterator *	iterator,
		const std::vector< std::string > &	species,
		unsigned int	windowSize,
		unsigned int	step,
		unsigned int	maxGap,
		unsigned int	maxPos,
		bool	keepTrashedBlocks,
		bool	missingAsGap
	)

inline

Create a new AlignmentFilter2MafIterator with absolute thresholds.

Parameters

iterator	Input iterator
species	Selection of species on which filtering criteria are applied. Results of filtering will be applied to all species.
windowSize	Size of the sliding window (nt).
step	Step by which windows are moved (nt).
maxGap	Maximum number of gaps allowed in the window.
maxPos	Maximum number of gaps "events" allowed.
keepTrashedBlocks	Removed windows are kept as separate blocks.
missingAsGap	Add missing species as gap sequences where needed.

Definition at line 208 of file AlignmentFilterMafIterator.h.

◆ AlignmentFilter2MafIterator() [2/2]

bpp::AlignmentFilter2MafIterator::AlignmentFilter2MafIterator	(	MafIterator *	iterator,
		const std::vector< std::string > &	species,
		unsigned int	windowSize,
		unsigned int	step,
		double	maxPropGap,
		unsigned int	maxPos,
		bool	keepTrashedBlocks,
		bool	missingAsGap
	)

inline

Create a new AlignmentFilterMafIterator with relative thresholds.

Parameters

iterator	Input iterator
species	Selection of species on which filtering criteria are applied. Results of filtering will be applied to all species.
windowSize	Size of the sliding window (nt).
step	Step by which windows are moved (nt).
maxPropGap	Maximum proportion of gaps allowed in the window.
maxPos	Maximum number of gaps "events" allowed.
keepTrashedBlocks	Removed windows are kept as separate blocks.
missingAsGap	Add missing species as gap sequences where needed.