micropan (version 1.0)

readBlastTable: Reading BLAST result file

Description

Reading a file produced by the BLAST+ software set up to produce tabular output.

Usage

readBlastTable(blast.file)

Arguments

blast.file

Name of file to read.

Value

The content of the file is returned as a data.frame with 12 columns and one row for each BLAST result. The columns have self-explanatory names.

Details

This function will read files produced by the BLAST+ software where the option -outfmt 6 has been invoked during its call. This option forces BLAST to produce a short tabular text output for each BLAST search. The function blastAllAll produces such files.

See Also

blastAllAll, bDist.

Examples

# Using a BLAST result file in this package
# We need to uncompress it first...
extdata.path <- file.path(path.package("micropan"),"extdata")
filenames <- "GID1_vs_GID2.txt"
pth <- lapply( file.path( extdata.path, paste( filenames, ".xz", sep="" ) ), xzuncompress )

#...then we can read it...
blast.table <- readBlastTable(file.path(extdata.path,filenames))

# ...and compressing it again...
pth <- lapply( file.path( extdata.path, filenames ), xzcompress )