The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Did you remember to?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

Did we ask for a bug report?


Then follow instructions in Article#1894.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community
Picard 2.9.4 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

ReduceReads - No real operator (M|I|D|N) in CIGAR

Hi. When I reduce my BAM files with ReduceReads and then run HaplotypeCaller I get this warning:

Ignoring SAM validation error: ERROR: Read name 1060, No real operator (M|I|D|N) in CIGAR

Why is this? Is it OK to inore it?

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    This sounds like the aligner was sleeping on the job... If it's just a few reads, I'd ignore it, or maybe use -rf BadCigar to filter them out so we don't risk choking on them later on. If it's more than a few, by all means check the alignments throughout the file.

  • Dear Geraldine,
    I encountered the same problem when running AddOrReplaceReadGroups tools from picard.Since the software exited when it met the first bad formatted read,I had no idea about the total number of such reads.Would you please explain more about the aligner was sleeping on the job?
    Thanks and regards.

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    Hi @taojinsheng,

    "sleeping on the job" is a colloquial phrase that means "doing something wrong". In this case it means that the aligner produced invalid output.

    One easy way of knowing how many reads are affected is to run the CountReads tool with the BadCigar filter enabled, like this:

    java -jar GATK.jar -T CountReads -R ref.fasta -I reads.bam -rf BadCigar
    

    When the run finishes, there will be a line in the console output that tells you how many reads were filtered out by BadCigar.

    If the percentage of bad reads is high, then you may need to do additional QC on your data and possibly realign them from the original FastQ.

  • I got it.Thanks a bunch,Geraldine.

Sign In or Register to comment.