Practical-3-Processing-16S-rRNA-amplicon-data.md

@@ -114,7 +114,7 @@ For now though, we're just going to use a (relatively) simple plotting function
     plotTaxon <- function(seqtab, taxtable, taxon, min_abund, taxonomicLevel) {
       seqtab <- data.frame(seqtab)
       colnames(seqtab) <- paste(taxtable[, 1], taxtable[, 2], taxtable[, 3], taxtable[, 4], 
-taxtable[, 5], taxtable[, 6], as.character(1:length(taxtable[, 1])), sep="-")
+      taxtable[, 5], taxtable[, 6], as.character(1:length(taxtable[, 1])), sep="-")
       seqtab <- seqtab/rowSums(seqtab)
       seqtab <- seqtab[, sapply(seqtab, function(x) max(x)) >= min_abund]
       seqtab$sample <- rownames(seqtab)