How to deal with no calls from GATK Unifiedgenotyper for indels

[audqf]

I am working with 454 single-end long-read data for a region on chr6, using hg19 and dbSNP build 135 as references. GATK Unifiedgenotyper with the option -glm INDEL returns no calls on realigned+recalibrated bam file and on recalibrated bam file. Using additional arguments, such as --output_mode EMIT_ALL_SITES, or setting call/emit thresholds, does not make a difference, either.

Earlier posts on this forum (http://seqanswers.com/forums/showthread.php?t=15634) indicate that a high sequencing error rate may be the cause.

So I wonder:
1. If resetting the sequencing error rate is the solution, how do I do it? Is this error rate reported in the bam file?

2. Are there other ways of tweaking GATK to deal with this problem?

Thanks for your help!

[ulz_peter]

I once tried to twaek GATK to output Indels in 454 data, but I think they turned that utility off as soon as you specify the PL:454 (or maybe it was PL:ROCHE, I don't remember) tag. In case you find a solution, I'd be happy to it

[audqf]

In case ulz_peter and others are interested, it turns out that GATK UG does not make indel calls for 454 data.

I also posted the same question here:

http://getsatisfaction.com/gsa/topics/how_to_deal_with_no_calls_from_gatk_unifiedgenotyper_for_indels?utm_content=reply_link&utm_medium=email&utm_source=reply_notification#reply_7912928

(btw, this is a pretty good place to have GATK-related questions answered)

Now the wiki page for GATK UG clearly states this limitation:

http://www.broadinstitute.org/gsa/wiki/index.php/Unified_genotyper

Hope this helps. I'm trying samtools for indel calling now.