Bowtie2 alignment using a large index (*.bt2l)
Entering edit mode
9 days ago
Abderhman • 0


I have used Bowtie2 to in the past successfully to index and align reads on R. I currently have a job that requires using large index. bowtie2_build returns successful index, but it has the unique bt2l extension (as mention in bowtie2 manual).

However, when I attempt to use this index for aligning using the bowtie2 function, the index is rejected because it's not a .bt2 extension:

-- bowtie2(bt2Index = "mm10.hg38_bowtie2", samOutput = "dspt5_1.sam", seq1 = "dspt5_1.fastq", seq2 = "dspt5_11.fastq", overwrite = FALSE)

-- Error in checkFileExist(paste0(bt2Index, ".1.bt2"), "bt2Index") : --- For argument bt2Index, file does not exist: mm10.hg38_bowtie2.1.bt2

I have attempted to find way to let bowtie2 know that I am using a large index (.bt2l), but couldn't. Anyone knows how I could get around this?

Thank you!!

Large R index Bioconductor Bowtie2 • 78 views

Login before adding your answer.

Traffic: 2973 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6