tommyau / bamclipper

Remove primer sequence from BAM alignments by soft-clipping
MIT License
31 stars 10 forks source link

clipping is not working if deletions follows the primer site #15

Open MarieLataretu opened 1 year ago

MarieLataretu commented 1 year ago

Hi there,

it seems like bamclipper misses bases if a deletion follows the primer site:

igv_snapshot_22290_covpipe2-clip-viralrecon_annotated

The corresponding bedpe line is this:

MN908947.3  22262   22290   MN908947.3  22626   22650

From my understanding, it should definitively be caught with -u 8 -d 8, so I guess the deletion is the issue.

Best Marie