SEA-PHAGES Logo

The official website of the HHMI Science Education Alliance-Phage Hunters Advancing Genomics and Evolutionary Science program.

Welcome to the forums at seaphages.org. Please feel free to ask any questions related to the SEA-PHAGES program. Any logged-in user may post new topics and reply to existing topics. If you'd like to see a new forum created, please contact us using our form or email us at info@seaphages.org.

New Features in PECAAN

| posted 16 Aug, 2019 22:54
Thanks Claire,

That is helpful and we had seen some hints in the HHpred.
The students had suggested that the only matches were all "integrases" without any further designation. I will have to take a closer look when we are back in class again. Thanks again for your time!
Best,
Heather
| posted 17 Aug, 2019 20:12
DrHHNZ
Thanks Claire,

That is helpful and we had seen some hints in the HHpred.
The students had suggested that the only matches were all "integrases" without any further designation. I will have to take a closer look when we are back in class again. Thanks again for your time!
Best,
Heather

Hi Heather - I think that is an NCBI issue. The hits used to say which they were, but I think more recently NCBI decided to use the generic "integrase" for everything so that information doesn't pop out like it used to.

Lee
| posted 18 Aug, 2019 20:03
I see you added Topcons now, too. Cool!

Steve
| posted 26 Jan, 2020 14:53
Is there a document explaining how to interpret Topcons data? I looked around the SEA-PHAGES forums, the Bioinformatics guide and found nothing. The Topcons paper is a little dense for me.

In particular I'm QC.ing AwesomeSauce - an F1 which is in Pecaan.
In that phage, gene 33 (start 28413) has that HHPred match to the Putative lactococcus lactis phage r1t holin which has no transmembrane domain when run through TmHmm. The Topcons summary shows different results depending on the program but the reliability box shows nothing (which is different from the easier to interpret membrane protein for gene 29 in this phage).

I interpreted this as no evidence for a transmembrane domain and would thus not call this the holin. Any insight on how to interpret teh Topcons summary for this gene would be appreciated!
| posted 29 Jan, 2020 03:40
TOPCONS has a very good prediction statistics compared to other prediction programs, including TmHmm. The five supporting methods, shown along with TOPCONS, aid in developing the reliability score graphic. To score high on reliability you need to have a TM domain predicted by TOPCONS and all of the other five supporting programs at the same location. The neat thing about TOPCONS is its ability to accurately predict when there is no TM domain. TOPCONS is over 95% accurate in predicting when there is no TM domain. So, if there is no TM domain found in the TOPCONS line then you can bet that it is not a membrane protein.

I looked at gene 34 (start 28413) in PECAAN and see TOPCONS evidence for two TM domains, which fits the two TM domain evidence for the phage r1t holin. At least one domain has decent support from three out of five of the other prediction programs, which is OK support. I think that there is enough evidence for the holin call.

Claire
 
Login to post a reply.