Entering edit mode
4.3 years ago
devarts
▴
40
Hi All. I'm trying to an alignment for multiple files, paired end, using Rsubread. I set up the readfile vector with file names as follows.
read1 <- list.files(path= "path to/Yarrowia_lipolytica/Genewiz-us-ngs-00_fastq", pattern = "*_R1_001.fastq.gz$")
read2 <- list.files(path= "path to/Genewiz-us-ngs-00_fastq", pattern = "*_R2_001.fastq.gz$")
That seems to work, in example reads1 appears to be chr [1:24] "F353-B-gku_R1_oo1.fasta.gz"....
But, when I try to run align, as follows, I get an error that is a mystery to me.
align(index = "index", readfile1 = reads1, readfile2 = reads2, output_file = "./Rsubread_alignment.BAM", input_format = "FASTQ", output_format = "BAM")
Error in normalizePath(files, mustWork=mustWork)
: path[1] = "F353-B-gku_R1_oo1.fasta.gz".: No such file or directory
I wonder if the file names need the whole path attached somehow?
Any help would be much appreciated.
Thanks, David
Please use the formatting bar (especially the
code
option) to present your post better. You can use backticks for inline code (`text` becomestext
), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.Thank You. I'm new.
Thank You. That worked.
Please only add answers when you're answering the question on the top-level post. This psot should have been a comment on my answer. I've moved it to a comment on the top-level post for now.
Please accept my answer (the green check mark on the left) so the post is marked as solved.
Whoops. My bad. I am learning.
That's OK. Welcome to biostars!