qcShortRead <- function(fastqs,reportOutDir=getwd()){ qaSamples <- lapply(fastqs,qa) names(qaSamples) <- gsub("\\.fastq.*|\\.fq.*|\\.fq|\\.fastq","",basename(fastqs)) reportOutFiles <- lapply(names(qaSamples),function(x){ report(qaSamples[[x]],dest = file.path(reportOutDir,x)) }) links <- file.path(reportOutDir,names(qaSamples)) links <- paste0("<a href='", (file.path(links,"index.html")),"'>", names(qaSamples),"</a>") qaSampleFrame <- do.call(ShortRead::rbind,qaSamples) fqQCTable <- data.frame(SampleNames=links, as.data.frame([email protected]$readCounts )) return(list(FQQC_Table = fqQCTable,ShortReadQC=qaSampleFrame)) }