[Bioc-sig-seq] Question when running ShortRead to generate a QA report
pterry at huskers.unl.edu
pterry at huskers.unl.edu
Wed Sep 9 00:09:21 CEST 2009
Dear Bioc-sig-sequencing,
My question: Referring to the ShortRead tutorial (IOQA.pdf) at Bioc2009, page 15, how to generate a 'narrative description QA file' like 'qa_080828_081110.pdf' for my data?
That is,
Using the developer version of R,
mterry at system76-pc:~/R-DEVEL/R-49517$ ./bin/R
R version 2.10.0 Under development (unstable) (2009-09-01 r49517)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
...
> library("ShortRead")
Loading required package: IRanges
...
Loading required package: Biostrings
Loading required package: BSgenome
Loading required package: lattice
> packageDescription("ShortRead")$Version
[1] "1.3.32"
>
> sp <- SolexaPath("extdata/ELAND/080828_HWI-EAS88_0003")
> rpt <- report(qa(sp), dest="reports/report_090709.pdf")
Note: p-15, of IOQA.pdf, ShortRead tutorial at BioC2009, refers to sample report 'qa_080828_081110.pdf'. Near bottom of page 15, it states: "The QA report is self-documenting, providing a narrative description of each section."
In the new "reports/report_090709.pdf/image" directory, there are 8 pairs of files plus a logo file:
-rw-r--r-- 1 mterry mterry 22143 2009-09-07 12:23 alignmentQuality.jpg
-rw-r--r-- 1 mterry mterry 15710 2009-09-07 12:23 alignmentQuality.pdf
-rw-r--r-- 1 mterry mterry 8473 2009-09-07 12:23 bioclogo-small.jpg
-rw-r--r-- 1 mterry mterry 23130 2009-09-07 12:23 perCycleBaseCall.jpg
-rw-r--r-- 1 mterry mterry 23011 2009-09-07 12:23 perCycleBaseCall.pdf
-rw-r--r-- 1 mterry mterry 20257 2009-09-07 12:23 perCycleQuality.jpg
-rw-r--r-- 1 mterry mterry 21864 2009-09-07 12:23 perCycleQuality.pdf
-rw-r--r-- 1 mterry mterry 29271 2009-09-07 12:23 perTileCount.jpg
-rw-r--r-- 1 mterry mterry 74448 2009-09-07 12:23 perTileCount.pdf
-rw-r--r-- 1 mterry mterry 17217 2009-09-07 12:23 perTileHistogram.jpg
-rw-r--r-- 1 mterry mterry 3785 2009-09-07 12:23 perTileHistogram.pdf
-rw-r--r-- 1 mterry mterry 27000 2009-09-07 12:23 perTileQuality.jpg
-rw-r--r-- 1 mterry mterry 69548 2009-09-07 12:23 perTileQuality.pdf
-rw-r--r-- 1 mterry mterry 21976 2009-09-07 12:23 readOccurences.jpg
-rw-r--r-- 1 mterry mterry 31303 2009-09-07 12:23 readOccurences.pdf
-rw-r--r-- 1 mterry mterry 19937 2009-09-07 12:23 readQuality.jpg
-rw-r--r-- 1 mterry mterry 56130 2009-09-07 12:23 readQuality.pdf
but I don't find "a self-documenting pdf file with narrative description of each section" as in the example file: qa_080828_081110.pdf? Before running this r49517 example, I had just installed 'texlive', (a TeX system) using Synaptic Package Manager onto my Ubuntu 9.04 linux laptop. Should this satisfy the requirement for 'LATEX installation' mentioned on page 15 of IOQA.PDF?
> sessionInfo()
R version 2.10.0 Under development (unstable) (2009-09-01 r49517)
x86_64-unknown-linux-gnu
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] ShortRead_1.3.32 lattice_0.17-25 BSgenome_1.13.11 Biostrings_2.13.36
[5] IRanges_1.3.66
loaded via a namespace (and not attached):
[1] Biobase_2.5.6 grid_2.10.0 hwriter_1.1 tools_2.10.0
>
P. Terry
pterry at huskers.unl.edu
More information about the Bioc-sig-sequencing
mailing list