diff rmarkdown_samtools_flagstat.Rmd @ 2:77e8e7bc5b44 draft

update
author mingchen0919
date Wed, 27 Dec 2017 15:19:43 -0500
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/rmarkdown_samtools_flagstat.Rmd	Wed Dec 27 15:19:43 2017 -0500
@@ -0,0 +1,32 @@
+---
+title: 'Reads alignment stats'
+output:
+    html_document:
+      number_sections: true
+      toc: true
+      theme: cosmo
+      highlight: tango
+---
+
+```{r setup, include=FALSE, warning=FALSE, message=FALSE}
+knitr::opts_chunk$set(
+  echo = ECHO,
+  error = TRUE
+)
+```
+
+
+# Reads alignment stats on `r opt$input_bam_name`
+
+```{r 'flagstat'}
+command = paste0('samtools flagstat ', opt$input_bam, ' > ', paste0(opt$report_dir, '/report.txt'))
+system(command)
+df = readLines(paste0(opt$report_dir, '/report.txt')) %>%
+  str_replace('(\\d+\\s+\\+\\s+\\d+)\\s+', '\\1,') %>% 
+  str_split(',') %>% 
+  as.data.frame() %>% t()
+rownames(df) = NULL
+colnames(df) = c('PASS + FAIL', 'Description')
+knitr::kable(df)
+```
+