Question: realignement around Indels
0
gravatar for Chadi Saad
22 months ago by
Chadi Saad60
France
Chadi Saad60 wrote:

Hello,

Is it possible to realign reads and force the aligner to put Indels at a certain position, instead of aligning reads to the left ? I found GATK realigner, but it demand an vcf file, while I have just the positions (without the REF / Alt fields in VCF).

alignment • 626 views
ADD COMMENTlink written 22 months ago by Chadi Saad60

Take a look at ABRA. It does realignments well.

ADD REPLYlink written 22 months ago by genomax67k

thanks, but how can I specify my prefered Indel Positions ? (Imagine I have an homopolymer of 5 A , with insertion of 1 A , and I need the aligner to put it in the middle of the homopolyper...)

ADD REPLYlink written 22 months ago by Chadi Saad60
1

If the realigner aligns all indels to the same position, why do you need it to be on a particular position? Any position will be arbitrary anyway.

ADD REPLYlink written 22 months ago by h.mon25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1013 users visited in the last hour