SUBREAD Aligning reads with reference genome - Error in align The number of input file names is different from the number of output file names.
0
0
Entering edit mode
saira • 0
@021c4eb8
Last seen 2.6 years ago
United Kingdom

Enter the body of text here

Code should be placed in three backticks as shown below


reads1 <- list.files( path = "path", pattern = "*_1.fastq$" )
reads2 <- list.files( path = "path", pattern = "*_2.fastq$" )

align(index = "ht", readfile1 = reads1, readfile2 = reads2, input_format = "FASTQ", output_format = "BAM", nthreads = 16)   **#this line  brings up error****

Error in align(index = "ht", readfile1 = reads1, readfile2 = reads2, input_format = "FASTQ", : The number of input file names is different from the number of output file names.

align Alignment subread • 879 views
ADD COMMENT

Login before adding your answer.

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