MafFilter: a highly flexible and extensible multiple genome alignment files processor

Julien Y. Dutheil (1, 2) , Sylvain Gaillard (3) , Eva H. Stukenbrock (4)

1 Department of Organismic Interactions [Marburg]
2 UMR ISEM - Institut des Sciences de l'Evolution de Montpellier
3 IRHS - Institut de Recherche en Horticulture et Semences
4 Fungal Biodiversity Group [Marburg]

Background: Sequence alignments are the starting point for most evolutionary and comparative analyses. Full genome sequences can be compared to study patterns of within and between species variation. Genome sequence alignments are complex structures containing information such as coordinates, quality scores and synteny structure, which are stored in Multiple Alignment Format (MAF) files. Processing these alignments therefore involves parsing and manipulating typically large MAF files in an efficient way. Results: MafFilter is a command-line driven program written in C++ that enables the processing of genome alignments stored in the Multiple Alignment Format in an efficient and extensible manner. It provides an extensive set of tools which can be parametrized and combined by the user via option files. We demonstrate the software's functionality and performance on several biological examples covering Primate genomics and fungal population genomics. Example analyses involve window-based alignment filtering, feature extractions and various statistics, phylogenetics and population genomics calculations. Conclusions: MafFilter is a highly efficient and flexible tool to analyse multiple genome alignments. By allowing the user to combine a large set of available methods, as well as designing his/her own, it enables the design of custom data filtering and analysis pipelines for genomic studies.

Domaines: Sciences du Vivant [q-bio]

Julien Y. Dutheil, Sylvain Gaillard, Eva H. Stukenbrock. MafFilter: a highly flexible and extensible multiple genome alignment files processor. BMC Genomics, 2014, 15, ⟨10.1186/1471-2164-15-53⟩. ⟨hal-01210010⟩