CRT home > TissueInfo: Manual: timegablast

[ synopsis ]

timegablast
--mb filename [--tiac filename] [--db database] [--evalue E-value] [--error allowed_sequencing_error] [--required-length length] [--assemble-hsps] [--best-est] [--required-length length] [--help]

[ description ]

Converts a megablast output to the tiac format. Reads a MegaBLAST output file, formatted with the -D 3 option, finds all the matches which satisfy the criteria to keep a match (i.e. --error and --required-length), and outputs these to an XML style tiac file. The --assemble-hsps option is recommended when mRNAs are used as the query sequences when running MegaBLAST.

[ options ]

--mb
MegaBLAST output filename
The input file should be a MegaBLAST output file, formatted with the -D 3 option.
Default value: none

--tiac
Output filename
This creates an output file with the filename given. If this filename does not already have a .tiac suffix, it will be automatically appended to it.
Default value: tiac

--db
Database name
The name of the database that was searched with MegaBLAST. This will be included in the description of the MegaBLAST search at the top of the tiac file, so that users can keep track of their parameters.
Default value: est

--evalue
Expectation value (double)
The expectation value used during the MegaBLAST search. As for --db, this will be included in the details at the top of the tiac file.
Default value: 1e-20

--error
Percentage sequence error (double)
Sets the maximum percentage sequence error permitted for a match to be considered.
Default value: 0.05 (i.e. 5% sequencing error tolerated)

--required-length
Minimum length (integer)
Sets the minimum HSP alignment length (between the query and the subject sequences) required for a match to be considered.
Default value: 150

--assemble-hsps

Accumulates the length of all HSPs from the same subject sequence and query sequence (which pass the sequence identity threshold).The total length of these HSPs is then used to determine if a match reaches the minimum length. This option is useful when processing output obtained from unspliced mRNA. In this case, exons in the query sequence often result in several HSPs individually shorter than the required alignment length. This option accumulates their length before applying the threshold.
Default value: false

--best-est

Instructs TImegablast to only associate an EST with a query sequence if that query sequence is the best match for that EST. It is advised that this option should be used in conjunction with --assemble-hsps.
Default value: false

--required-length
Minimum length (integer)
Minimum alignment length required to consider an hsp.
Default value: 150

--help

Print a help message.
Default value: false

News
Jul, 2009; ChIPseeqer, a comprehensive framework for analysis of ChIP-seq data developed in the Elemento lab, is now available for download. [More]
Apr, 2009; The BDVal program developed by the Campagne laboratory for MAQC-II is now available from http://bdval.org. The software supports the development and evaluation of predictive biomarker models from high-throughput data. The web site offers binary and source distributions. [More]
Jan, 2009; Twease now supports searching MEDLINE articles by Author, Journal, and Publication Year. Examples for performing these searches can be found in the updated Twease tutorial. [More]

[News Archives] [Mailing List]


Events
Dec 11th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Fabien Campagne; ICB Conference Room - Y.1301
Jan 15th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Lei Shi; ICB Conference Room - Y.1301
Feb 12th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Christopher E. Mason; ICB Conference Room - Y.1301
Mar 12th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Olivier Elemento; ICB Conference Room - Y.1301
Apr 9th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Emre Aksay; ICB Conference Room - Y.1301
May 14th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Jonathan D. Victor; ICB Conference Room - Y.1301
Jun 11th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Harel Weinstein; ICB Conference Room - Y.1301
Jul 9th; 4:00pm-5:00pm: Institute for Computational Biomedicine Research in Progress Seminar Series - Duane Hassane; ICB Conference Room - Y.1301