Dpuiu alignment

Filtering

Programs:

• filter-anc.pl -max -percent -w -W

 max:      max distance between the ref_intercept (default 1,000,000)
 percent:  min percent agreement for keeping an alignment (default 25)
 w:        min window size (default 4)   => 1+w*2
 W:        max window size (default 100) => 1+W*2

 Bacterial genomes (few Mbp) use -max 100000  -W 10
 Mammalian genomes (few Gbp) use -max 1000000 -W 100 
 

• merge-anc.pl -max

 max:      max distance between the ref_intercept (default 1,000,000)

• Example:

 cat PA14-PA7.delta     | ~/bin/shrinkIds.pl | ~/bin/DELTA/delta2anc.pl > PA14-PA7.anc 
 cat PA14-PA7.anc       | ~/bin/DELTA/filter-anc.pl -max 100000 -W 10   > PA14-PA7.filter.anc
 PA14-PA7.filter.anc    | ~/bin/DELTA/merge-anc.pl -max 100000          > PA14-PA7.merge.anc
 cat PA14-PA7.merge.anc | ~/bin/DELTA/anc2delta.pl                      > PA14-PA7.merge.delta

 wc -l *anc
  691 PA14-PA7.anc
  528 PA14-PA7.filter.anc
    9 PA14-PA7.merge.anc

---

Pairwise Alignment Processing (poster)

Daniela Puiu*, Arthur Delcher, Steven Salzberg University of Maryland, College Park

Abstract

We present a software tool for context based sequence alignment processing. The tool was initially developed for finding syntenic regions between large eukariotic genomes but could be also used for analyzing rearrangements in microbial genomes. The sofware takes as input alignment output generated by di

Motivation

Sequence alignment is one of the oldest bioinformatics problems. Due to its complexity there is no one general solution and multiple algorithms have been implemented for aligning different data set types. There is always a tradeoff between speed and accuracy.