# HG changeset patch
# User dktanwar
# Date 1513614006 18000
# Node ID ce3ad612a104b0cf5d666270b2b7643725cdbc71
# Parent  ed929b85b051a59e8046c9b8f48741f9efce00bd
Uploaded

diff -r ed929b85b051 -r ce3ad612a104 csaw.R
--- a/csaw.R	Mon Dec 18 10:58:39 2017 -0500
+++ b/csaw.R	Mon Dec 18 11:20:06 2017 -0500
@@ -15,6 +15,7 @@
 library("stringr")
 library("data.table")
 library("getopt")
+library("Rsamtools")
 
 
 options(stringAsfactors = FALSE, useFancyQuotes = FALSE)
@@ -29,6 +30,9 @@
 # Separate multiple input files into a list of individual files
 files <- unlist(strsplit(inputs, ','))
 
+# Index bamfiles
+indexBam(files = files)
+
 # Create windows and count reads in them ----
 Sys.time()
 windows <- windowCounts(files, spacing=150, width=200, bin=F)