Biostar Beta. Not for public use.
realignement around Indels
0
Entering edit mode
23 months ago
Chadi Saad • 60
France

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 • 651 views
ADD COMMENTlink
0
Entering edit mode

Take a look at ABRA. It does realignments well.

ADD REPLYlink
0
Entering edit mode

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
1
Entering edit mode

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

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1