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.

All posts created by DanRussell

| posted 03 Mar, 2016 18:47
Hey Joe,

What 1:1 means is that the first amino acid of the query sequence matches the first amino acid of the subject sequence. This implies that the start site of your query sequence matches the start site of the published subject sequence.

Or, like you said, a 1:4 alignment would mean that the first amino acid of the query matches the fourth amino acid of the subject. This means there are three more amino acids in the published subject sequence that are not in your query. This can indicate that you're choosing a different start than the published one.

I made a picture that hopefully helps show these different situations. Remember that the grayed-out bases in this picture won't appear in BLAST results, because they're not part of the alignment. But I included them in this picture to give a better sense of how the protein sequences are really lining up.



–Dan
Edited 03 Mar, 2016 18:49
Posted in: Choosing Start SitesUsing BLAST to validate start sites
| posted 03 Mar, 2016 16:15
Hi Joe,

I think the value of using BLAST when choosing start sites lies in the fact that you're comparing your choice to the final choices made by other (often experienced) people. So it's not biological evidence, but it's sort of consensus evidence. You could argue, of course, that simply because someone else has called it a certain way is no reason to call it that way too, and the skeptic in all of us resists conforming to the norm. But many of the genomes in GenBank have been very carefully QCed and pored over by those who do this work more than anyone else.

So I think of it kind of like checking your work with an expert. If you BLAST, and get 1:1 hits with many published things, it's like, "Oh, good, seems like everyone agrees with this start call. That makes me a little more confident in it." If all the BLAST hits are not 1:1, it's like "Hmmm…looks like everyone else called a different start for some possibly good reason. Maybe I should double-check and really buttress my argument on this one."

Two things to keep in mind when using BLAST results for start sites:
  1. Who submitted this genome to GenBank? Was it Debbie/Welkin/Graham? Or was it Jimmy FirstPhage? Give more weight to trusted experts.
  2. When was this genome submitted? We have more evidence now for making start site calls than we did in the past, so give more weight to newer genomes than older ones.


Hope that helps,
–Dan
Posted in: Choosing Start SitesUsing BLAST to validate start sites
| posted 03 Mar, 2016 16:05
saleadon
Hi Everyone,
In "Exploring Bacteriophage Biology", one suggested bioinformatics experiment is to design cluster-specific DNA primers. The article mentions that there is a list of existing primers on phagesdb.org under “Phorum”. However, I can't find that. Has it been moved or deleted?
Thanks,
Steve

Hi Steve,

That document needs to be updated! Yes, the PhagesDB Phorum has indeed been deleted. And I'm not sure that we'd necessarily recommend that particular experiment either, since it has proved of limited use. It might be an interesting teaching/thought experiment…like: can you locate regions that are conserved in all the members of a subcluster? A cluster? Or not? (i.e., Is it even possible to design cluster-specific primers?)

In any event, I know that Alex Peister at Morehouse College has been the de facto go-to person for cluster primers, so she would be a good resource if you're interested in the current state of those attempts.

–Dan
Posted in: Phage BiologyPCR Primers
| posted 02 Mar, 2016 21:49
Hi Denise,

Sounds like that would be enough for me to take a look! Let me know when they're posted and I'll see what I can see.

–Dan
Posted in: AnnotationLocating Terminase Gene
| posted 29 Feb, 2016 19:48
Debbie Jacobs-Sera
Hi all! I was just able to blast a gene (and repeated by Dan) in DNA Master. The Integer Overflow error is gone and my posted wait time was 31 seconds, with the blast data returned in about 1 1/2 minutes. I just started to Blast an entire genome and will keep you posted as to how that goes! We are back in business!

Hopefully this is fixed, everyone! Nice job by Debbie bugging NCBI into cooperation.

Let us know if you continue to have problems. You shouldn't need to update anything to have success, just try again.
–Dan
Posted in: DNA MasterBLAST in DNAM
| posted 25 Feb, 2016 18:34
Hi Denise,

First of all, are you sure this is a circularly permuted genome? It's always possible that the phage already has ends of its own, then you don't have to worry about choosing them yourself. I can check that if you give me the sequence file and the sequencing reads. (Posting on Dropbox or Google Drive is a good way to move those large files.)

If it truly is circularly permuted, I'm sure Debbie would be happy to take a look at potential Base 1 options. She's done this a lot!

–Dan
Posted in: AnnotationLocating Terminase Gene
| posted 25 Feb, 2016 16:03
Hi Denise,

Does this happen to be Beth or C3PO?

–Dan
Posted in: AnnotationLocating Terminase Gene
| posted 25 Feb, 2016 15:25
Sebastian Schmeier
Please the link as well Dan. Are there different links for Mac vs Win versions? I would need both.

Can you please also sent me the student and faculty user/passwd please.

Thank you very much,

Sebastian

~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~
Sebastian Schmeier
Lecturer in Bioinformatics/Genomics
Institute of Natural and Mathematical Sciences
Massey University Auckland, New Zealand
+64 9 213 6538 :: +64 9 414 0800 ext 43538
s.schmeier@massey.ac.nz :: http://sschmeier.com
~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~-~

Just sent them.

–Dan
Posted in: SEA-PHAGES Virtual MachineStudent download of 2016 VM
| posted 24 Feb, 2016 20:35
Jeffrey Lawrence is aware of the issue and is working on it now. If he comes up with a fix and/or if there's anything to report, we'll let you know as soon as we do!

–Dan
Posted in: DNA MasterBLAST in DNAM
| posted 24 Feb, 2016 14:49
Hi Keith et al,

Yes, I'm seeing some slower peformance for sure in the VM DNA Master. For some functions, it's quite quick, but for others it drags. I hit a snag when I tried to import a fasta file that was about 170 kb long. It froze, and never made it through. I tried repeating, adjusting VM memory, etc., but nothing seemed to help. That would obviously be a dealbreaker, but perhaps we can figure out some workarounds or improvements. The 2017 SEA VM will use the newer version of Ubuntu, so perhaps those upgrades will help.

I'm curious as to whether some of the problems might relate to the artificial PC file system that's set up by wine, and whether we can change some settings to improve along those lines.

–Dan
Posted in: DNA MasterBeta Test: DNA Master within the SEA VM