To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

Indexing not recognized by alignment scripts

Hello,

I've followed the exact commandline instructions for two separate reference files (hg.fa and Homo_sapiens_assembly19.fasta). I am definitely generating the files, but other tools (bwa mem) find no index when I use the reference as an argument. Full example below:

$ bwa index -a bwtsw -p Homo_sapiens_assembly19 Homo_sapiens_assembly19.fasta $ samtools faidx Homo_sapiens_assembly19.fasta $ picard CreateSequenceDictionary \ REFERENCE=Homo_sapiens_assembly19.fasta \ OUTPUT=Homo_sapiens_assembly19.dict
(all of these finish without an error)

$ ls Homo_sapiens_assembly19* Homo_sapiens_assembly19.amb Homo_sapiens_assembly19.bwt Homo_sapiens_assembly19.fasta Homo_sapiens_assembly19.pac Homo_sapiens_assembly19.ann Homo_sapiens_assembly19.dict Homo_sapiens_assembly19.fasta.fai Homo_sapiens_assembly19.sa

However

$ bwa mem Homo_sapiens_assembly19.fasta filtered.fastq [E::bwa_idx_load_from_disk] fail to locate the index files

Any ideas? I could swear I've run this exact same pipeline before with the same files without any issues...

I'm running bwa-0.7.12 with picard-1.140 and samtools-1.2 on a Macbook Pro OSX El Capitan. Every command executed within the same directory.

Thanks for any help!

Tagged:

Answers

  • jfb325jfb325 IthacaMember

    Resolved. The -p flag in the bwa index line should have the same text as the first argument of bwa mem. Working after using:

    bwa mem Homo_sapiens_assembly19 filtered.fastq
    instead of
    bwa mem Homo_sapiens_assembly19.fasta filtered.fastq

    or using
    bwa index -a bwtsw Homo_sapiens_assembly19.fasta
    instead of
    bwa index -a bwtsw -p Homo_sapiens_assembly19 Homo_sapiens_assembly19.fasta

Sign In or Register to comment.