Search
Question: QuasR error usign preprocessReads
0
gravatar for Paul Ríos
2.3 years ago by
Paul Ríos0
Perú
Paul Ríos0 wrote:

I have fastq files from a paired end sequencing, so I would like to prepare the input sequence(remove adapters and N bases) before to do the alignment. So I ran the next commands but it shows a error:

infiles <- system.file(package="QuasR", "extdata",
                       c("nsc0001_S1_L001_R1_001.fastq.gz", "nsc0001_S1_L001_R2_001.fastq.gz"))

outfiles <- paste(tempfile(pattern=c("output_1.fastq","output_2.fastq")),sep="")

preprocessReads(filename=infiles[1],
                outputFilename="output_1.fastq",
                filenameMate=infiles[2],
                outputFilenameMate=outfiles[2],
                nBases=0, complexity=0.6)

Error in preprocessReads(filename = infiles[1], outputFilename = "output_1.fastq",  : 
  unsupported file format (must be one of 'fasta' or 'fastq')

Could someone help me with this?

Thank you

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by Paul Ríos0
1
gravatar for Mike Smith
2.3 years ago by
Mike Smith2.7k
EMBL Heidelberg / de.NBI
Mike Smith2.7k wrote:

Admittedly the error message is misleading, but the files you're trying to read don't exist.

> library(QuasR)
> infiles <- system.file(package="QuasR", "extdata", c("nsc0001_S1_L001_R1_001.fastq.gz", "nsc0001_S1_L001_R2_001.fastq.gz"))
> infiles
[1] ""

You more likely want something like:

infiles <- file.path("/here/is/my/data", c("nsc0001_S1_L001_R1_001.fastq.gz", "nsc0001_S1_L001_R2_001.fastq.gz"))

where "/here/is/my/data" is the path to the folder containing the fastq files you're trying to analyse.

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by Mike Smith2.7k
0
gravatar for Paul Ríos
2.3 years ago by
Paul Ríos0
Perú
Paul Ríos0 wrote:

Thank you for the clarification! That works.

ADD COMMENTlink written 2.3 years ago by Paul Ríos0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 136 users visited in the last hour