Kevin Meza Landeros

Script; obtencion de anotacion de zonas de interes

......@@ -31,7 +31,7 @@ mcols(homoS) = mcols(homoS)[,c("source","type","ID","Name")]
Despues de haber trabajado nuestros datos con bowtie, utilizamos el archivo resultante .bam
#cargamos nuestros datos de las enfermedades monogenicas.
bamFile <- "/home/aschafer/Documentos/Genomicas/Semestre_4/Genomica_humana/Enfermedades_monogenicas/sequences_aligned_A_sort.bam"
bamFile <- "/home/aschafer/Documentos/Genomicas/Semestre_4/Genomica_humana/Enfermedades_monogenicas/sequences_aligned_sort.bam"
informacion <- c("rname", "strand", "pos", "qwidth")
informacion <- ScanBamParam(what = informacion)
bam <- scanBam(bamFile, param = informacion)
......