Alec Ingros
February 15, 2017
Rational of today’s work: Annotating genes: 4, 7, 8, 9, 20, and 11
Tools used and/or Methods:
- HHpred
- GeneMark
- Starterator
- PhagesDB
- Phamerator
- NCBI blast
Results:
Gene 4
Start: 1620bp Stop: 2729bp FWD GAP: 20bp Gap SD Final Value: SD Score: -2.077 (Best score) Z-Value: 3.188 CP: The gene is not covered The gene can’t be extended to cover all – approximently 20 bp left uncovered. SCS: Agrees with Glimmer, Agrees with GeneMark NCBI BLAST: endolysin [Arthrobacter phage Gordon] Q:1 S:1 E-Value: 0.0 CDD: PGRP and Cpl-7 domain-containing protein E-Value: 2.2e-13 PhagesDB BLAST: Gordon_4, endolysin Q:1 S:1 E-Value: 1e-164 HHPred: Endolysin, putative lys E-Value: 2.3e-23 LO: Yes ST: Agrees with Starterator F: FS: HHpred, NCBI, PhagesDB Notes: Endolysin; PGRP and Cpl-7 domain-containing protein
Gene 4 GeneMark Results
Gene 4 PhagesDB Best Results
Gene 4 NCBI Blast Best Results
Gene 4 NCBI Blast Conserved Domain Results
Gene 4 HHpred Best Results
Gene 4 Starterator Results
Gene 7
Start: 2886bp Stop: 3269bp FWD GAP: 8bp Overlap SD Final Value: SD Score: -2.904 (Best score) Z-Value: 2.767 CP: The gene is not covered The gene can’t be extended to cover all – approximently SCS: Agrees with Glimmer, Agrees with GeneMark NCBI BLAST: hypothetical protein GORDON_5 [Arthrobacter phage Gordon] Q:6 S:39 E-Value: 2.0e-50 CDD: No good hit PhagesDB BLAST: Gordon_5, function unknown Q:6 S:39 E-Value: 6e-41 HHPred: No good hit LO: Yes ST: We cannot shorten the gene without losing the best z-score, SD value, and longest open reading frame F: NKF FS: NCBI, PhagesDB Notes:
Gene 7 GeneMark Results
Gene 7 PhagesDB Best Results
Gene 7 NCBI Blast Best Results
Gene 7 HHpred Results
Gene 8
Start: 3266bp Stop: 3793bp FWD GAP: 4bp Overlap SD Final Value: SD Score: -3.207 (Best score) Z-Value: 3.047 CP: The gene is not covered The gene can’t be extended to cover all – approximently SCS: Agrees with Glimmer, Agrees with GeneMark NCBI BLAST: hypothetical protein GORDON_6 Q:18 S:1 E-Value: 2e-74 CDD: No good hit PhagesDB BLAST: Gordon_6, function unknown Q:18 S:1 E-Value: 6e-59 HHPred: No good hit LO: Yes ST: We cannot shorten the gene without losing the best z-score, SD value, and longest open reading frame F: NKF FS: NCBI, PhagesDB Notes:
Gene 8 GeneMark Results
Gene 8 PhagesDB Best Results
Gene 8 NCBI Blast Best Results
Gene 8 Starterator Results
Gene 9
Start: 3812bp Stop: 5551bp FWD GAP: 18bp Gap SD Final Value: SD Score: -2.137 (Best score) Z-Value: 3.188 CP: The gene is covered SCS: Agrees with Glimmer, Agrees with GeneMark NCBI BLAST: terminase large subunit [Arthrobacter phage Gordon] Q:1 S:1 E-Value: CDD: Phage terminase-like protein, large subunit, contains N-terminal HTH domain [Mobilome E-Value: 3.99e-49 PhagesDB BLAST: Gordon_7, terminase, large subunit Q:1 S:1 E-Value: 0.0 HHPred: Terminase, DNA packagin E-Value: 1.1e-27 LO: Yes ST: Agrees with Starterator F: terminase, large subunit FS: NCBI, PhagesDB, HHpred Notes:
Gene 9 HHpred Best Results
Gene 9 PhagesDB Best Results
Gene 9 NCBI Blast Best Results
Gene 9 NCBI Blast Conserved Domain Results
Gene 10
Start: 5548bp Stop: 6399bp FWD GAP: 4bp Overlap SD Final Value: SD Score: -2.451 (Best score) Z-Value: 3.043 CP: The gene is covered SCS: Agrees with Glimmer, Agrees with GeneMark NCBI BLAST: hypothetical protein CAPNMURICA_8 Q:1 S:1 E-Value: 0.0 CDD: No good hit PhagesDB BLAST: CapnMurica_8, function unknown Q:1 S:1 E-Value: 1e-152 HHPred: No good hit LO: Yes ST: Agrees with Starterator F: NKF FS: NCBI, PhagesDB Notes:
Gene 10 GeneMark Results
Gene 10 PhagesDB Best Results
Gene 10 NCBI Blast Best Results
Gene 10 HHpred Best Results
Gene 10 Starterator Results
Gene 11
Start: 6396bp Stop: 6977bp FWD GAP: 4bp Overlap SD Final Value: SD Score: -4.241 (Best score) Z-Value: 2.184 CP: The gene is covered SCS: Agrees with Glimmer, Agrees with GeneMark NCBI BLAST: hypothetical protein GORDON_9 Q:5 S:4 E-Value: 1e-79 CDD: No good hit PhagesDB BLAST: Niktson_Draft_8, function unknown Q:1 S:1 E-Value: 6e-75 HHPred: No good hit LO: Yes ST: Agrees with Starterator F: NKF FS: NCBI, Phages DB Notes:
Gene 11 GeneMark Results
Gene 11 PhagesDB Best Results
Gene 11 NCBI Blast Best Results
Gene 11 HHpred Best Results
Gene 11 Starterator Results
Conclusions and Next Steps:
- We skipped over genes 5 and 6 because we thought it would be best to decide as a team what to do. The situation is that neither gene is definitively better than the other, so a decision should be made as a group. The genes we did finish were put us back on track to finish in the two week window. We will look to have another productive day like today to continue with schedule.