search

TimoLassmann / kalign

A fast multiple sequence alignment program.
GNU General Public License v3.0
124 stars 29 forks source link

how to adjust the penalty parameter #29

Closed huangnengCSU closed 1 year ago

huangnengCSU commented 2 years ago

Hi, I use kalign to perform MSA for DNA sequence. I used the parameter kalign -i reads.fa -f clu --tgpe 0.01 --gpo 0.2 --gpe 0.4. The input fasta file is :

>1
GGGGGAAGGGGG
>2
GGGGAAGGGGG
>3
GGGGAAGGGGG
>4
GGGGAAGGGG
>5
GGGGGAAGGGGGG
>6
GGGGGAAGGGG
>7
GGGGGAAGGGGGA
>8
GGGGGAAGGGGGG
>9
GGGGGAAGGGGGG
>10
GGGGAAGTTG
>11
A
>12
GGGGGGAAGGGGG
>13
GGGGGAAGGGGGG
>14
GGGAGT
>15
GGGGAAGGGGG

The figure below is the alignment result of clustalw.aln. How can I adjust the parameter to get the similar alignment result as clustalw. image

Best, Neng

TimoLassmann commented 1 year ago

Hi, I added alignment types to version 3.3.4. The dna alignment type uses standard substitution, mismatch and gap penalties and will produce a sensible alignment in the case above.