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!

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!


Surround blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block.
Powered by Vanilla. Made with Bootstrap.
Picard 2.9.0 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.

index reduced bam files

blueskypyblueskypy Member Posts: 266 ✭✭

can I use samtools index to index reduced bam files? I'm using GATK version 2.7-2-g6bda569 and HC asks for the index files.

Thanks,

Tagged:

Answers

  • LaviniaLavinia Member Posts: 37

    Try it, what does it do, does HC like the resulting index?

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,388 admin

    That should be okay. Reduced bams are valid according to spec, so samtools should handle them correctly re: indexing.

    Geraldine Van der Auwera, PhD

  • blueskypyblueskypy Member Posts: 266 ✭✭

    No complaint from HC. I'm just not sure the error msgs I got from VQSR are related to this. See http://gatk.vanillaforums.com/discussion/comment/9737. btw, I don't remember HC asked for that index file before.

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,388 admin

    GATK requires that all bam files be indexed, but all GATK produce index files automatically, so usually you don't see it requested by any tools apart from the very first in your pipeline. The question is, what happened to the index file produced by ReduceReads?

    Geraldine Van der Auwera, PhD

  • blueskypyblueskypy Member Posts: 266 ✭✭
    edited October 2013

    Ah ha, what a dummy I am! I didn't realize ReduceReads also produces an index file. No wonder HC never asked for the index file before because I kept all those files at each step. But this time I thought I could delete all of them and only kept the last reduced bam file, ... :(

    Well, learn from mistakes! But I'm sorry for wasting your time, Geraldine and Lavinia!

    btw, I think it might be good to add that note to the docs of ReduceReads, so others won't make the same mistake.

Sign In or Register to comment.