Question: RealignerTargetCreator inputs multiple Bams??
0
gravatar for micro32uvas
3 months ago by
Pakistan
micro32uvas0 wrote:

Hello Everyone,

I have been working with RealignerTargetCreator and in its pages, RealignerTargetCreator. It states that it takes multiple inputs as Bam. But its not taking it. when I checked Galaxy.it stated that GenomeAnalysisTK: RealignerTargetCreator accepts an aligned BAM input file.

I want to give all my Bam files for indel realigner intervals.

Here's the command that I am using; java -jar GenomeAnalysisTK.jar -T RealignerTargetCreator -R reference.fasta -I input.bam -known indels.vcf -o forIndelRealigner.intervals

Your help is appreciated. Qurat

ADD COMMENTlink modified 3 months ago by Devon Ryan70k • written 3 months ago by micro32uvas0
1
gravatar for Pearl
3 months ago by
Pearl30
UK
Pearl30 wrote:

For using multiple bam files as inputs to RrealignerTargetCreator, you can use "-I input1.bam -I input2.bam -I input3.bam" in the command.

ADD COMMENTlink written 3 months ago by Pearl30

This is in Galaxy, micro32uvas can't issue commands.

ADD REPLYlink written 3 months ago by Devon Ryan70k

You just solved my problem Pearl :) I'll get back here, if there's still some errors

ADD REPLYlink modified 3 months ago • written 3 months ago by micro32uvas0
0
gravatar for Devon Ryan
3 months ago by
Devon Ryan70k
Freiburg, Germany
Devon Ryan70k wrote:

There's a separate biostar site for Galaxy. Anyway, I would strongly encourage you NOT to use GATK in Galaxy, it doesn't follow GATK's best practices. Use freeBayes or samtools instead.

BTW, what you're trying to do isn't possible in the Galaxy wrapper.

ADD COMMENTlink modified 3 months ago • written 3 months ago by Devon Ryan70k

I never prefer Galaxy anyways, and yes its very limited as well. Its just I was missing the -I switch everytime, which was creating the problem. I just checked in galaxy for the given options there, if i was missing any.

Thanks Devon

ADD REPLYlink written 3 months ago by micro32uvas0

Ah, if you're already using the command line then use the HaplotypeCaller and then skip this step entirely.

ADD REPLYlink written 3 months ago by Devon Ryan70k

Can you kindly elaborate this Devon? Actually i am working on Haplotypecaller but along with this preprocessing phase. i.e. after marking duplicates, now i am up with realign target creator, indel realigner and then base calibration. Then would be the haplotype calling using haplotype caller.

Your help at this stage would be epic for my work

P.S. I am already on command line.

ADD REPLYlink modified 3 months ago • written 3 months ago by micro32uvas0

If you're using the HaplotypeCaller then you don't need to align around indels. There's also no need to perform base quality recalibration any more. This will save you from wasting a LOT of time for no benefit.

ADD REPLYlink written 3 months ago by Devon Ryan70k

i cannot understand how would this omit realignment. wouldn't realigner give further improvement in terms of mapping and enhanced calibration. If you could precisely share the command that would cover all of it, it would be highly appreciated

ADD REPLYlink written 3 months ago by micro32uvas0

Why realignment isn't needed requires explaining how the HaplotypeCaller and indel realignment work. I'll just link to this article on the GATK site, which goes over the process in more detail than would fit in a comment. Suffice it to say, indel realignment is already part of the HaplotypeCaller.

Regarding base recalibration, the head to head comparisons that Brad Chapman has shown (have a look at his blog, he documents everything there) that this isn't really beneficial any more. I know it's still in the GATK "best practices", but those seem to only be slowly updated.

ADD REPLYlink written 3 months ago by Devon Ryan70k
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: 1278 users visited in the last hour