run bwa on rstudio
0
0
Entering edit mode
6 days ago
Nemo • 0

Hello,

I am looking to find a way to run bwa (Burrows-Wheeler Aligner) on my fastq files on rstudio. I have uploaded the bwa zip file (downloaded from Fast and accurate long-read alignment with Burrows-Wheeler transform) in rstudio and tried to run the following command:

bwa.path <- "/folder1/user/bwa/"

system(paste0(bwa.path," bwa index -a bwtsw /folder1/user/ref1.fasta"))

However, it keeps getting me the error of No such file or directory. Is this the standard way to call bwa command on rstudio? (

fastq Aligner bwa rstudio index • 235 views
ADD COMMENT
0
Entering edit mode

Remove the whitespace in front of bwa in the system command but seriously, run this in terminal or via a workflow manager. This is nothing to be done from inside R. Non-standard and nothing to get used to.

ADD REPLY
0
Entering edit mode

Thanks @ATpoint for your reply. Removing white space does not solve the problem. (I am asked to do it via R.)

ADD REPLY
1
Entering edit mode

/folder1 does not look like any standard path in the root dir in either Linux or macOS. That is probably wrong. cd into the folder with the bwa binary and type pwd, use that path. Same for the fasta.

ADD REPLY

Login before adding your answer.

Traffic: 1735 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