atac - comparative mapping between genome sequences
atac -meryldir alignment_cache_dir
-dir result_dir -id1 seqid1
-seq1 seqfile1.fasta -id2 seqid2
-seq2 seqfile2.fasta
atac computes a one-to-one pairwise alignment of large DNA sequences. It first
finds the unique k-mers in each sequence, chains them to larger blocks, and
fills in spaces between blocks. It was written primarily to transfer
annotations between different assemblies of the human genome.
The output is a set of ungapped 'matches', and a set of gapped 'runs' formed
from the matches. Each match or run associates one sequence with the other
sequence. The association is 'unique', in that there is no other (sizeable)
associations for either sequence. Thus, large repeats and duplications are not
present in the output - they appear as unmapped regions.
Though the output is always pairwise, atac can cache intermediate results to
speed a comparisons of multiple sequences.
Run with no options for a list of parameters.
leaff(1),
/usr/share/doc/atac/README.atac
http://kmer.sourceforge.net/wiki/index.php/Getting_Started_with_ATAC