Question: GATK IndelRealigner make mistakes
gravatar for whaiyu06
21 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 • 781 views
ADD COMMENTlink modified 21 months ago by Pierre Lindenbaum124k • written 21 months ago by whaiyu0620
gravatar for Pierre Lindenbaum
21 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum124k 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 21 months ago • written 21 months ago by Pierre Lindenbaum124k

Thanks.There nothing output as your command.

ADD REPLYlink written 21 months ago by whaiyu0620

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

ADD REPLYlink written 20 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: 1934 users visited in the last hour