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

hello,

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 https://www.broadinstitute.org/gatk/guide/article?id=1317for 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
0
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.

Help
Access

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