This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BWA: Importing a third-party indexed genome

I guess my question is going to sound pretty basic, but nevertheless I need to understand clearly how to use bwa .

I am not able to run a bwa indexing of the GRCh38 reference genome on my laptop, due to memory usage requirements (~5GB needed, according to the bwa man page).

The 1000 Genomes Project offers what seems to be an indexed version of it on its FTP servers.

Can I download the entire directory and use bwa locally for my alignments, without running a brand new indexing?

bwa alignment grch38 1000genomes
• 2,425 views • link
updated by Biostar Biostar • written by GD • GD

1 answer

If you are not able to create an index then you would not be able run the alignments either since that also typically requires 5G of memory. See this note from Heng.

But the answer for your original question is yes. It is possible to use pre-formatted indexes that you get from someplace else.

• 1 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this
  • Which 1000 genomes 30x files should I use for imputation?
    written by Apprentice Apprentice
    I would like to do an imputation using data of 1000 genomes phase 3 30x (GRCh38) as a reference panel. 1000 genomes phase 3 30x …
  • BCFtools consensus
    written by Peerzada • Peerzada
    I am trying to incorporate variants from multisample vcf file using bcftools consensus as :- First download variant and reference files as ``` bcftools view …
  • Extract a gene region from a whole genome cram file from 1000 genomes
    written by michael.flower.14 • michael.flower.14
    I'm trying to download whole genome sequencing data from the 1000 genomes project to pad out some GVCF files I have for use in variant …
  • 1000 genomes missing SNPs on chr X in GRCh38
    written by b_ha • b_ha
    I am trying to get a reference together for imputation using the 1000 genomes data set. The GRCh38 version of the data set was made …
  • differences between Illumina Platinum Variant Calls and NIST variant calls
    written by kspata kspata
    Hi All, I am validating a bioinformatics pipeline for SNP and INDEL calling. For this purpose I mapped the reads from Illumina Platinum Genome ([https://www.ebi.ac.uk/ena/data/view/ERR194147][1]) …
  • 1000 genomes project mitochondrial DNA sequences
    written by bbb • bbb
    Help! How (where) to get mitochondrial dna sequences from "1000 genomes" project? This is what I found: http://www.1000genomes.org/announcements/phase3-mitochondrial-chromosome-variants-now-available-2015-10-06/ http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/README_chrMT_phase3_callmom.md I need fasta files of individual …
  • Human genome full chromosomes build
    written by priess1991 • priess1991
    Hi! I need to use some data of the 1000 human genome project. On 1000genomes.org's ftp-server I found all data in form of 14gb big …
  • News: 1000 Genomes data on GRCh38
    written by fairley • fairley
    The low coverage 1000 Genomes sequence data has been realigned to GRCh38. Reads were aligned to the full assembly, including the GRC maintained alternate loci …
  • 1000 Genomes - Sample Subset Of A Vcf File
    written by pablo.riesgo • pablo.riesgo
    <p>Hi there,</p> <p>I'm trying to retrieve genotype data for a given sample from the 1000 genomes FTP repository, following their guidelines: <a href="http://www.1000genomes.org/faq/how-do-i-get-sub-section-vcf-file">http://www.1000genomes.org/faq/how-do-i-get-sub-section-vcf-file</a></p> <p>So I …
  • Problems Accessing 1000G Data
    written by Ahdf-Lell-Kocks • Ahdf-Lell-Kocks
    <p>I am just copy+pasting from the 1000g instructions here: <br/> <a href="http://www.1000genomes.org/faq/how-do-i-get-sub-section-bam-file">http://www.1000genomes.org/faq/how-do-i-get-sub-section-bam-file</a></p> <pre><code>samtools view -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data/HG00154/alignment/HG00154.mapped.ILLUMINA.bwa.GBR.low_coverage.20101123.bam 17:7512445-7513455 </code></pre> <p>But I can't get it to work:</p> …