Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Terminal gaps not considered as InDel #2

Open
Close-your-eyes opened this issue Dec 15, 2022 · 1 comment
Open

Terminal gaps not considered as InDel #2

Close-your-eyes opened this issue Dec 15, 2022 · 1 comment

Comments

@Close-your-eyes
Copy link

Close-your-eyes commented Dec 15, 2022

The alignment below induces terminal gaps in the subject. This is not indicated as an indel. So, from looking at the meta-data of the alignment only, one may think that it is a 'perfect' alignment.

Is that behaviour expected or desired?

aln <- pwalign::pairwiseAlignment("CCGGATCG", "ATCG", type = "global-local")
aln

Global-Local PairwiseAlignmentsSingleSubject (1 of 1)
pattern:     CCGGATCG
subject: [1] ----ATCG
score: -18.07298 

pwalign::indel(aln)

An object of class "InDel"
Slot "insertion":
IRangesList object of length 1:
[[1]]
IRanges object with 0 ranges and 0 metadata columns:
       start       end     width
   <integer> <integer> <integer>


Slot "deletion":
IRangesList object of length 1:
[[1]]
IRanges object with 0 ranges and 0 metadata columns:
       start       end     width
   <integer> <integer> <integer>
@acvill
Copy link

acvill commented Jun 9, 2023

An older, related post from the Bioconductor forum:

indels located at the edges of subject and Biostrings pairwiseAlignment function

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants