Skip to content

HTTPS clone URL

Subversion checkout URL

You can clone with
or
.
Download ZIP
new codebase of the PaPaRa algorithm
C++ C Other

This branch is 1 commit ahead, 106 commits behind master

Fetching latest commit…

Cannot retrieve the latest commit at this time

Failed to load latest commit information.
ivy_mike @ 41473cc
old
ublasJama-1.0.2.3
.gitmodules
CMakeLists.txt
README
align_pvec_vec.h
align_utils.h
align_vec.h
aligned_buffer.h
build_papara2.sh
dtw.h
epa_extract_qs_covered.cpp
fasta.h
fasta_random_sample.cpp
math_approx.h
pairwise_seq_distance.cpp
pairwise_seq_distance.h
papara.cpp
pars_align_gapp_seq.cpp
pars_align_gapp_seq.h
pars_align_seq.cpp
pars_align_seq.h
parsimony.cpp
parsimony.h
pvec.h
pw_dist.cpp
raw_to_phy.rb
raxml_interface.cpp
raxml_interface.h
release_papara2.sh
sequence_model.cpp
sequence_model.h
small.tree
smith_waterman.cpp
stepwise_align.h
test_bitset.cpp
tree_similarity.cpp
tree_similarity.h
tree_utils.h
vec_unit.h

README

Implementation of the PaPaRa 2.0 algorithm.
Build with 'sh build_papara2.sh' which will produce the executable 'papara'.
Alternatively you can use the supplied cmake files. To be able to compile the sources you need a resonably recent version
of the boost (www.boost.org) libraries. 

Invoke PaPaRa using "./papara -t <ref tree> -s <phylip RA> -q <fasta QS>".

The phylip file (option -s) must contain the reference alignment, consistent with the reference tree (option -t).
The FASTA file (option -q) contains the unaligned QS. Optionally, all sequences which are in <phylip RA> but do not occur in the <ref tree> are also interpreted as QS. 
The alignment parameters can be modified using the (optional) option -p <user_options>. <user options> is a string and must have the following form:
"<gap_open>:<gap_extend>:<mismatch>:<match_cgap>", so the default parameters used given in the paper correspond to the user option "-p -3:-1:2:-3".  

The output alignment will be written to papara_alignment.default (you can change the file suffix (i.e., "default") by supplying a run-name with parameter '-n'.
You can invoke the multi threaded version by adding the option '-j <num threads>'. 

The latest source code is available at https://github.com/sim82/papara_nt
Something went wrong with that request. Please try again.