IniziaInizia gratis

Calcolare la coverage

Ora tocca a te: estendi le read fino alla lunghezza media dei frammenti per ottenere una stima della coverage più accurata. Estendere le read ti permette di combinare le read provenienti da entrambe le estremità di un frammento in un unico picco di coverage, aiutandoti a distinguere i picchi dal rumore di fondo. Il file BAM per il campione GSM1598218 che userai in questo esercizio è già stato caricato. La lunghezza media dei frammenti per questo campione è di 183 bp. Le read sono disponibili come oggetto GRanges chiamato reads_gr.

Il caricamento di tutti i dati e dei pacchetti R necessari per questo esercizio potrebbe richiedere qualche istante. Ti chiediamo un po' di pazienza.

Questo esercizio fa parte del corso

ChIP-seq con Bioconductor in R

Visualizza il corso

Istruzioni dell'esercizio

  • Estendi le read alla lunghezza media dei frammenti.
  • Calcola la coverage dei frammenti.

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

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

# Compute coverage
cover <- ___(___)
Modifica ed esegui il codice