\name{qa}
\alias{qa}
\alias{qa,character-method}
\alias{qa,list-method}

\title{Perform quality assessment on short reads}

\description{

  This function is a common interface to quality assessment functions
  available in \code{ShortRead}. Results from this function may be
  displayed in brief, or integrated into reports using, e.g.,
  \code{\link{report}}.

}
\usage{
qa(dirPath, ...)
\S4method{qa}{character}(dirPath, pattern=character(0), 
    type=c("SolexaExport", "SolexaRealign", "Bowtie", "MAQMap", 
           "MAQMapShort", "fastq", "BAM"),
    ...)
\S4method{qa}{list}(dirPath, ...)
}

\arguments{

  \item{dirPath}{A character vector or other object (e.g.,
    \code{\link{SolexaPath}}; see \code{showMethods}, below) locating
    the data for which quality assessment is to be performed. See help
    pages for defined methods (by evaluating the example code, below)
    for details of available methods.}

  \item{pattern}{A character vector limiting the files in \code{dirPath}
    to be processed, as with \code{\link{list.files}}. Care should be
    taken to specify pattern to avoid reading unintended files.}

  \item{type}{The type of file being parsed; must be a character vector
    of length 1, selected from one of the types enumerated in the
    parameter.}

  \item{\dots}{Additional arguments used by methods. 

    \describe{

      \item{\code{fapply}, \code{reduce}:}{Influence how evaluation
        occurs when this function is run with the \pkg{Rmpi} or
        \pkg{parallel} packages; see \code{\link{srapply}}.}

      \item{\code{Lpattern}, \code{Rpattern}:}{A character vector or
        XString object to be matched to the left end of a sequence. If
        either \code{Lpattern} or \code{Rpattern} are provided,
        \code{trimLRPatterns} is invoked to produce a measure of adapter
        contamination.  Mismatch rates are 0.1 on the left and 0.2 on
        the right, with a minimum overlap of 10 nt.}

    }
  }
}

\details{
  The following methods are defined, in addition to those on S4 formal
  classes documented elsewhere:

  \describe{

    \item{\code{qa,character-method}}{

      Quality assessment is performed on all files in directory
      \code{dirPath} whose file name matches \code{pattern}. The type of
      analysis performed is based on the \code{type} argument. Use
      \code{SolexaExport} when all files matching \code{pattern} are
      Solexa \code{_export.txt} files. Use \code{SolexaRealign} for
      Solexa \code{_realign.txt} files. Use \code{Bowtie} for Bowtie
      files. Use \code{MAQMapShort} for MAQ \code{map} files produced by
      MAQ versions below 0.70 and \code{MAQMap} for more recent output.
      Use \code{fastq} for collections of fastq-format files. Use
      \code{BAM} for \code{.bam} files; see the \code{param} argument in
      \code{\link{readAligned}} for information on how to influence
      reads included in the qa summary. Quality assessment details vary
      depending on data source.

    }
    
    \item{\code{qa,list-method}}{

      \code{dirPath} is a list of objects, all of the same class and
      typically derived from \code{ShortReadQ}, on which quality
      assessment is performed. All elements of the list must have names,
      and these should be unique.

    }
  }
}

\value{

  An object derived from class \code{\linkS4class{.QA}}. Values
  contained in this object are meant for use by \code{\link{report}}

}

\author{Martin Morgan <mtmorgan@fhcrc.org>}

\seealso{
  \code{\linkS4class{.QA}},
  \code{\linkS4class{SolexaExportQA}}
  \code{\linkS4class{MAQMapQA}}
  \code{\linkS4class{FastqQA}}
  \code{\linkS4class{BAMQA}}
}

\examples{

sp <- SolexaPath(system.file("extdata", package = "ShortRead"))
qa1 <- qa(sp)
fq <- file.path(analysisPath(sp), pattern="s_1_sequence.txt")
qa2 <- qa(readFastq(fq), basename(fq))  # name fastq files

showMethods("qa", where=getNamespace("ShortRead"))
}

\keyword{manip}