探索 fastq 檔案
Fastq 檔案通常包含成千上萬甚至上百萬個 reads,檔案體積會非常大!在這個練習中,你會使用一個只有 500 個 reads 的小型 fastq 子樣本,能夠輕鬆載入記憶體,並可用 readFastq() 一次完整讀入。
原始的序列檔案來自阿拉伯芥(Arabidopsis thaliana),由 UC Davis Genome Center 提供。其 accession number 為 SRR1971253,下載自 Sequence Read Archive(SRA)。它包含來自葉片組織的 DNA,經過混合後使用 Illumina HiSeq 2000 定序。這些序列為單端讀取(single-read),長度為 50 個鹼基(bp)。
fqsample 是一個 ShortReadQ 物件,包含 reads、品質分數與 id 的相關資訊。輪到你來探索它了!
本練習屬於課程