Local realignment around indels
1
0
Entering edit mode
19 months ago
Elisa • 0

Hi everyone,

I'm new in bioinformatics field. I'd like to know how to and which are the best tools for local realignment around indels. I have fastq files coming from MiSeq (Illumina) as pipeline input and I'm performing the bioinformatic pipeline using python (and running tools from comand line). I should perform local realignment only after PCR duplicates removal? And which tool should I use to first identify indels around which local realignment will be done?

Thanks in advance, Elisa

local realignment • 476 views
ADD COMMENT
0
Entering edit mode
19 months ago

I should perform local realignment only after PCR duplicates removal?

if you're calling the variant with GATK4, you don't need to realign the indel because haplotypeCaller will do it on the fly. Realignment disappeared in gatk4

ADD COMMENT
0
Entering edit mode

I was calling the variants with bcftools. Do you suggest me to use GATK4 instead? (I have a targeted sequencing problem with tumor samples)

ADD REPLY
0
Entering edit mode

when publishing it's (sadly?) better to say you used the "GATK best practices" to call your variants...

ADD REPLY

Login before adding your answer.

Traffic: 1212 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6