# copy example data to current working directory
file.copy(system.file(package="QuasR", "extdata"), ".", recursive=TRUE)
# create alignments
sampleFile <- "extdata/samples_chip_single.txt"
genomeFile <- "extdata/hg19sub.fa"
auxFile <- "extdata/auxiliaries.txt"
proj <- qAlign(sampleFile, genomeFile, auxiliaryFile=auxFile)
proj
# alignment statistics using a qProject
alignmentStats(proj)
# alignment statistics using bam files
alignmentStats(alignments(proj)$genome$FileName)
alignmentStats(unlist(alignments(proj)$aux))
Run the code above in your browser using DataLab