local realignment without known indels

gilgigilgi Posts: 19Member
edited July 2012 in Ask the GATK team

Dear GATK team,

I understood from the best practice documentation that if I have no known indels list for the genome I can skip the realignment step. Is this correct, or shall I call Indels and then do realignment and then re-call SNPs and Indels with the UnifiedGenotyper?

Thanks for the help!

Best Answer

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 8,173Administrator, GATK Dev admin
    Answer ✓

    Hi there,

    You can run those tools on a merged bam file, or by passing several separate bam files in your command line (see the documentation if you're not sure how to do this).

    There are no hard and fast numbers for file sizes -- it depends enormously on your dataset and on your computing hardware. I can only recommend you test a few different configurations and see what works for you.

    Geraldine Van der Auwera, PhD

Answers

  • gilgigilgi Posts: 19Member

    I understand now that I need to run the RealignerTargetCreator and then the IndelRealigner tool.
    I am working on yeast barcoded targeted genome.
    The RealignerTargetCreator and the IndelRealigner are not read group aware - right?
    Would you suggest to run them on each sample separately - or on a merged file?

    One more question regarding the work on merged bam files - what is the maximal size of merged bam you recommend to work with?

    Many thanks for all the help!

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 8,173Administrator, GATK Dev admin
    Answer ✓

    Hi there,

    You can run those tools on a merged bam file, or by passing several separate bam files in your command line (see the documentation if you're not sure how to do this).

    There are no hard and fast numbers for file sizes -- it depends enormously on your dataset and on your computing hardware. I can only recommend you test a few different configurations and see what works for you.

    Geraldine Van der Auwera, PhD

  • gilgigilgi Posts: 19Member

    Thanks!

Sign In or Register to comment.