Give summary stats without merging

Author: bensutherland

01_scripts/account_reads.R

@@ -0,0 +1,18 @@
+# Clear space
+# rm(list=ls())
+
+## Install and load packages
+# NA
+
+# Set working directory
+setwd("~/Documents/04_HAB/eDNA_metabarcoding_HAB_18S_ngsfilter")
+
+# Read in data
+my.data <- read.csv2(file = "07_results/input_fastq_read_count.txt", header = F, sep = ",")
+
+head(my.data)
+
+# Summarize
+sum(my.data$V1)
+summary(my.data$V1)
+sd(my.data$V1)

01_scripts/account_reads.sh

@@ -0,0 +1,4 @@
+#!/bin/bash
+# Count number of reads per sample in 02_raw_data 
+for i in 02_raw_data/*R1_001.fastq ; do grep -cE '^\+$' $i ; done > 07_results/input_fastq_read_count.txt
+