[Bioc-sig-seq] perTile QA element in the FastqQA class

Sirisha Sunkara ssunkara at lbl.gov
Mon Apr 12 23:49:56 CEST 2010


Hi Martin,

The qa function that reads in fastq format files, doesn't seem to 
populate the perTile QA element with row information...
The row counts are zero for both the readCounts and 
medianReadQualityScore list elements of perTile.

Is this feature still work in progress..? Essentially, I am trying to 
get the TileQC plots for lanes where there was no reference genome to 
align (no export.txt files)

 >  qafq <- qa("./Contam_Screening/Run703/","s_8_1_sequence.txt", 
type="fastq")
 > qafq
class: FastqQA(9)
QA elements (access with qa[["elt"]]):
  readCounts: data.frame(1 3)
  baseCalls: data.frame(1 5)
  readQualityScore: data.frame(512 4)
  baseQuality: data.frame(94 3)
  alignQuality: data.frame(1 3)
  frequentSequences: data.frame(50 4)
  sequenceDistribution: data.frame(1663 4)
  perCycle: list(2)
    baseCall: data.frame(150 4)
    quality: data.frame(1081 5)
  perTile: list(2)
    readCounts: data.frame(0 4)
    medianReadQualityScore: data.frame(0 4)

Thank You,
Sirisha

 > sessionInfo()
R version 2.11.0 Under development (unstable) (2010-03-07 r51225)
x86_64-unknown-linux-gnu

locale:
[1] C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base    

other attached packages:
[1] ShortRead_1.5.21    lattice_0.18-3      Biostrings_2.15.22
[4] GenomicRanges_0.1.0 IRanges_1.5.74      Rmpi_0.5-8        

loaded via a namespace (and not attached):
[1] Biobase_2.7.5 grid_2.11.0   hwriter_1.2   tools_2.11.0



More information about the Bioc-sig-sequencing mailing list