[BioC] ShortReadQ constructor

[Martin Morgan]

On 10/07/2013 12:28 PM, Ivan Gregoretti wrote:
> Hello,
>
> Could somebody show an example of constructing a ShortReadQ object
> instance? I can't find an example in the documentation.
>
>
> For example, lets say that I have two records with the following ids,
> sequences, and Phred qualities
>
> c("read1", "read2")
> c("ACGT", "GGGGACT")
> c("HHHH", "HHHHHHH")
>
> How would you now construct the ShortReadQ instance?

Hi Ivan

ShortReadQ(DNAStringSet(c("ACGT", "GGGGACT")),
     FastqQuality(c("HHHH", "HHHHHHH")),
     BStringSet(c("read1", "read2")))

Martin

>
> Thank you,
>
> Ivan
>
>
>
>
> Ivan Gregoretti, PhD
> Bioinformatics
>
> _______________________________________________
> 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
>


-- 
Computational Biology / Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N.
PO Box 19024 Seattle, WA 98109

Location: Arnold Building M1 B861
Phone: (206) 667-2793