Question: GATK IndelRealigner make mistakes
gravatar for whaiyu06
7 months ago by
whaiyu0620 wrote:


I used the command to realignment around indels.The command as follows:

java -jar /Software/GenomeAnalysisTK-3.6/GenomeAnalysisTK.jar \
-R $ref \
-T IndelRealigner \
-targetIntervals A01-2m-B8_merge.mark.intervals \
-I  A01-2m-B8_merge.mark.bam  \
-o  A01-2m-B8_merge.mark.realn.bam \
-known /dbsnp/Mills_and_1000G_gold_standard.indels.b37.vcfs \
-known /dbsnp/1000G_phase1.indels.b37.vcf 2>stderr/A01-2m-B8_merge.GATK_IndelRealigner.stderr

and met the error message:

SAM/BAM/CRAM file htsjdk.samtools.SamReader$PrimitiveSamReaderToSamReaderAdapter@15f11bfb is malformed. Please see more information. Error details: the BAM file has a read with no stored bases (i.e. it uses '*') which is not supported in the GATK; see the --filter_bases_not_stored argument. Offender: E00477:238:hggvjccxy:7:2217:26890:33762

What should I do ?Thanks very much

sequencing alignment • 314 views
ADD COMMENTlink modified 7 months ago by Pierre Lindenbaum112k • written 7 months ago by whaiyu0620
gravatar for Pierre Lindenbaum
7 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum112k wrote:

what is the output of

samtools view A01-2m-B8_merge.mark.bam | grep -F 'E00477:238:hggvjccxy:7:2217:26890:33762'

as stated by GATK , I suspect that it contains a 'NO-BASES/*' flag in the read sequence. GATK is not able to handle this.

ADD COMMENTlink modified 7 months ago • written 7 months ago by Pierre Lindenbaum112k

Thanks.There nothing output as your command.

ADD REPLYlink written 7 months ago by whaiyu0620

I solved this issue by adding another parameter -filterNoBases. Thanks for your advice.

ADD REPLYlink written 6 months ago by whaiyu0620
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 998 users visited in the last hour