parasail
parasail copied to clipboard
Question only - pairwise local alignment without gaps
What would be the command line parameters for parasail_aligner
for doing local pairwise alignment (without gaps within the sequences) of two DNA sequences
For example, 1.fa
is like this
>1
GCCCCAGTAGGCAACACGCTGGCTGGCTGGCTGTCCGGGC
and 2.fa
is like this:
>2
GCCCCAGTAGGCAACACGCTGGCTGGCTGGCTGTCCGGGC
and I want to get them aligned like this:
ACGGTGACCAGGGTTCCCTGGCCCCAGTAGGCAACACGCTGG------------------
--------------------GCCCCAGTAGGCAACACGCTGGCTGGCTGGCTGTCCGGGC