inst/testScripts/bwa/11.bwa,index,level3.R

In HenrikBengtsson/aroma.seq: High-Throughput Sequence Analysis using the Aroma Framework

############################################################################
#
############################################################################
library("aroma.seq");

# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# Indexing a reference genome
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
path <- "annotationData/organisms/Lambda_phage";
fa <- FastaReferenceFile("lambda_virus.fa", path=path);
print(fa);

is <- buildBwaIndexSet(fa, verbose=TRUE);
print(is);


############################################################################
# HISTORY:
# 2012-09-24
# o Created.
############################################################################