compareToBlast: Compare results of SSAHA database search to results of the...

In janmatula/GenDP: Tools for genotyping data processing

Description

Compares top results of SSAHA database search to top results of BLAST database search. Outputs the DSC parameter which can be used to test the validity of SSAHA search results as compared to BLAST.

Usage

compareToBlast(SSAHAresults, blastResultsFile, results = 10)

Arguments

SSAHAresults	list of top 10 results of obtained by using the searchHashTableFunction.
blastResultsFile	.txt file containing results of BLAST database search obtained using the -outfmt 0 command.
results	number of top results of the SSAHA metod for each query sequence that will be compared to BLAST search results

Value

Vectorcontaining the DSC percentage parameter for each query sequence.

Examples

DSC <- compareToBlast(s11_rev_SSAHAresults, "Blast_s11_rev.txt", results =10)

janmatula/GenDP documentation built on Jan. 29, 2021, 7:29 a.m.