Has anyone else ever randomly ran into the problem of misencoded quality scores with the GATK? This happened once before and now just recently again on a particular sample. The sample comes from the same batch as several other exomes. They were put through the identical pipeline. One ends up giving me an error about a misencoded quality (65) but seems to work fine if I add the fixMisencodedQuals flag. This is when using the DepthOfCoverage walker after marking Duplicates with Picard.
Pipeline is Trimmomatic -> BWA-MEM -> Picard SamFormatConverter -> Picard AddOrreplaceReadGroups -> Picard ReorderSam -> Picard MarkDuplicates
Two samples from the same batch went through the same pipeline. Anyone ever see quality scores get randomly misencoded like that before?
0 answers
No answers yet.
Log in to answer this question.
-
Does the order of SplitNCigarReads and MarkDuplicates affect RNA-seq variant calling results?
written by iamsmor •
-
No real operator (M|I|D|N) in CIGAR
written by boymin2020 •
-
Empty output from Picard's EstimateLibraryComplexity
written by rubic
-
MarkDuplicates memory issue
written by cacampbell •
-
Mark duplicates for single end reads: Why only 5'end?
written by dariober
-
Gatk Error: Samfilereader Appears To Be Using The Wrong Encoding For Quality Scores
written by Charles Warden
-
Markduplicates Creating A Loss Of Mate Pair
written by DG
-
Calculate Aligment Metrics For Ngs Data (Zero Mismatches And Uniquely Mappable Reads)
written by User6891
-
Can You Make A Bam File Picard Compatible? I Forgot -M With Bwa Mem
written by ajc8 •
-
Preprocessing The Bam For Gatk
written by Pierre Lindenbaum
