What it does
A single genomic region can be covered by multiple alignment blocks. In many cases it is desirable to stitch these alignment blocks together. This tool accepts a list of genomic intervals. For every interval it performs the following:
- finds all MAF blocks that overlap the interval;
- sorts MAF blocks by alignment score;
- stitches blocks together and resolves overlaps based on alignment score;
- outputs alignments in FASTA format.
Example
Here three MAF blocks overlapping a single interval are stitched together. Space between blocks 2 and 3 is filled with gaps: