[BioC] fastq files

Hervé Pagès hpages at fhcrc.org
Fri Mar 11 19:09:43 CET 2011


Hi David, Martin,

Actually read.DNAStringSet() has the 'nrec' and 'skip' arguments to let
you read a subset of the FASTQ records.

Cheers,
H.


On 03/11/2011 05:48 AM, Martin Morgan wrote:
> On 03/11/2011 05:42 AM, David martin wrote:
>> Hi,
>> I have the following fastq header:
>>
>> @ILLUMINA-E6F4B9_0000:5:1:1184:949#0/1
>> NGGGCTGGGTCGGTCGGGC
>> +ILLUMINA-E6F4B9_0000:5:1:1184:949#0/1
>> BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
>>
>> Does any of the R packages (e.g: shortread) can convert this fastq to
>> fasta file ?
>
> Hi David --
>
> In Biostrings
>
>    write.XStringSet(read.XStringSet(fromFile, format="fastq"), toFile)
>
> in ShortRead
>
>    writeFasta(readFastq(fromFile), toFile)
>
> both read the entire file into memory.
>
> Martin
>
>>
>> thanks,
>> david
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>


-- 
Hervé Pagès

Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M2-B876
P.O. Box 19024
Seattle, WA 98109-1024

E-mail: hpages at fhcrc.org
Phone:  (206) 667-5791
Fax:    (206) 667-1319



More information about the Bioconductor mailing list