(ERR): "GRCh38_noalt_as" does not exist or is not a Bowtie 2 index Exiting now ...
1
0
Entering edit mode
6 weeks ago
Song • 0

Hi All,

I've been struggling with trying to align my fastq files to the Bowtie Indexed human genome. I followed this tutorial https://www.metagenomics.wiki/tools/short-read/remove-host-sequences and am now stuck where bowtie doesn't recognize the indexed database that I downloaded, which is a ready to use bowtie2 databse.

Here is my shell script:

#!/bin/bash
#SBATCH --nodes=1   # number of nodes
#SBATCH --ntasks-per-node=8   # 16 processor core(s) per node X 2 threads per core
#SBATCH --mem=500G   # maximum memory per node
#SBATCH --partition=mem   # standard node(s)
#SBATCH --job-name=%j
#SBATCH --mail-user=songkue@uga.edu   # email address
#SBATCH --mail-type=BEGIN
#SBATCH --mail-type=END
#SBATCH --mail-type=FAIL
#SBATCH --output=%j.out# job standard output file (%j replaced by job id)
#SBATCH --error=%j.err # job standard error file (%j replaced by job id) date                       

cd /project/metagenomics_of_retail_poultry/20230920/trimmed_fastq/trimmed_files

module load bowtie2/2.5.2

bowtie2 -p 8 -x GRCh38_noalt_as -1 1-446345050_R1.trim.fastq -2 1-446345050_R2.trim.fastq --very-sensitive-local --un-conc-gz
--reorder SAMPLE_host > 1-446345050_mapped_and_unmapped.sam

I keep getting this error:

(ERR): "GRCh38_noalt_as" does not exist or is not a Bowtie 2 index
Exiting now ...

I have inspected the bowtie2 index and it looked fine.

Any help will be greatly appreciated! Thank you!

alignment genomics bowtie bowtie2 • 418 views
ADD COMMENT
1
Entering edit mode
6 weeks ago
Ram 43k

-x needs the full prefix to the index files. Did you do the wget + unzip steps in that tutorial in the same folder (/project/metagenomics_of_retail_poultry/20230920/trimmed_fastq/trimmed_files) that you're in right now? If not, provide the full path to the place where you did the wget and unzip.

ADD COMMENT
0
Entering edit mode

I figured it out and found my mistake! I had to specify the folder and basename, for example GRCh38_noalt_as/GRCh38_noalt_as

ADD REPLY

Login before adding your answer.

Traffic: 2217 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6