Team:VT-ENSIMAG/Genothreat
From 2010.igem.org
Line 13: | Line 13: | ||
Following the federal guideline, the first step is to do the six frame translation ([[Team:VT_Ensimag_2010-Biosecurity/sixframe|See More]]) in order to have the corresponding amino-acid strands. We also keep the initial nucleotid strand and its reversed DNA sequence ([[/Why?]]). So we have now 8 sequences to screen (6 amino-acid, 2 nucleotides). | Following the federal guideline, the first step is to do the six frame translation ([[Team:VT_Ensimag_2010-Biosecurity/sixframe|See More]]) in order to have the corresponding amino-acid strands. We also keep the initial nucleotid strand and its reversed DNA sequence ([[/Why?]]). So we have now 8 sequences to screen (6 amino-acid, 2 nucleotides). | ||
Then, we divide each sequence in 200bps nucleotide or 66 amino-acid subsequences. | Then, we divide each sequence in 200bps nucleotide or 66 amino-acid subsequences. | ||
- | |||
<br> | <br> | ||
Each sub-sequence is then blasted, and we extract the relevant information from the blast output. | Each sub-sequence is then blasted, and we extract the relevant information from the blast output. | ||
- | For each sub-sequence, we have to decide if it's a hit or not. | + | For each sub-sequence, we have to decide if it's a hit or not ([[Team:VT_Ensimag_2010-Biosecurity/BestMatch|See More]]). |
+ | [[Image:VTENSI_Algo.jpg|center|200px]] | ||
+ | If one sub-sequence is a hit for our algorithm, then it means that a 200bps sequence potiently dangerous has been detected. The sequence will be so flagged, and a result file containing all the select agent found and some more information will be created. | ||
=Different implementations= | =Different implementations= | ||
}} | }} |
Revision as of 19:30, 9 August 2010
GenoTHREAT
|
Global presentationThe software, GenoTHREAT takes a DNA sequence on input, and in return told the user if it's a dangerous sequence which need further investigation or if the sequence can be synthetised. It follows the algorithm given by the federal guidelines, and so uses the Blast software for sequence alignement. For managing the sequence (six-frame translation, reversed frame, extracting sub-sequence), we use Biojava, a bio-informatics toolbox ([http://www.biojava.org/wiki/Main_Page 1]). GenoTHREAT has been coded in Java, and is working on both Windows and Linux distribution. We have developped different versions in order to optimize the time execution or the cpus' utilisation according to the use (#Different implementations). AlgorithmFollowing the federal guideline, the first step is to do the six frame translation (See More) in order to have the corresponding amino-acid strands. We also keep the initial nucleotid strand and its reversed DNA sequence (/Why?). So we have now 8 sequences to screen (6 amino-acid, 2 nucleotides).
Then, we divide each sequence in 200bps nucleotide or 66 amino-acid subsequences.
If one sub-sequence is a hit for our algorithm, then it means that a 200bps sequence potiently dangerous has been detected. The sequence will be so flagged, and a result file containing all the select agent found and some more information will be created.
Different implementations |