MulaiMulai sekarang secara gratis

Menghitung coverage

Sekarang saatnya Anda memperluas read hingga panjang fragmen rata-rata untuk memperoleh estimasi coverage yang lebih baik. Memperluas read memungkinkan Anda menggabungkan read dari kedua ujung suatu fragmen menjadi satu puncak coverage, sehingga membantu membedakan puncak dari derau latar. Berkas BAM untuk sampel GSM1598218 yang akan Anda gunakan dalam latihan ini sudah dimuat untuk Anda. Panjang fragmen rata-rata untuk sampel ini adalah 183 bp. Read tersedia sebagai objek GRanges bernama reads_gr.

Memuat semua data dan paket R yang diperlukan untuk latihan ini mungkin memerlukan waktu. Mohon bersabar.

Latihan ini adalah bagian dari kursus

ChIP-seq dengan Bioconductor di R

Lihat Kursus

Petunjuk latihan

  • Perluas read hingga panjang fragmen rata-rata.
  • Hitung coverage fragmen.

Latihan interaktif praktis

Cobalah latihan ini dengan menyelesaikan kode contoh berikut.

# Extend reads to the average fragment length of 183 bp
reads_ext <- ___(reads_gr, width=___)

# Compute coverage
cover <- ___(___)
Edit dan Jalankan Kode