[Seqinr-forum] dist.alignment gap penalty

[PENEL SIMON]

Hi Stephanie,

thanks for spotting this
The  gap penalty option is indeed disabled for dna sequences,  I will change this and let you now as soon as this is done, 

all the best


________________________________________
De : Seqinr-forum <seqinr-forum-bounces at lists.r-forge.r-project.org> de la part de Stephanie Young <stefunnie_young at yahoo.com>
Envoyé : jeudi 9 mars 2023 05:28
À : seqinr-forum at lists.r-forge.r-project.org
Objet : [Seqinr-forum] dist.alignment gap penalty

Hi,

I noticed that the most recent version of dist.alignment() function doesn't produce any differences when gap = F vs when gap = T.

Here is a toy example:

test1 <- msa(c(test1 = "acgtg", test2 = "acgtgaa"), type = "dna", method = "ClustalOmega") |>
    msaConvert(type = "seqinr::alignment") |>
    dist.alignment(gap = F)

test2 <- msa(c(test1 = "acgtg", test2 = "acgtgaa"), type = "dna", method = "ClustalOmega") |>
    msaConvert(type = "seqinr::alignment") |>
    dist.alignment(gap = T)

> identical(test1, test2)
[1] TRUE

Is there a way to impose a gap penalty or is this function option now defunct?

Thank you,
Stephanie