Link to this post | posted 09 Feb, 2024 21:52

debbie

Hi Emily, Both are good questions. For the first gene, the coding potential on the GeneMark output is rather robust upstream of 37. How would you evaluate the ability of that gene to start if you are evaluating a start at base 1. What do you want to see upstream of a start to evaluate the start? How would you do that in this case? Why do you say that RBS only favors the second start? For the second gene, how is that gene going to be made? What favors the best production of that gene? Keep thinking! What did you professor have to say? best, debbie Edited 09 Feb, 2024 21:54

Link to this post | posted 09 Feb, 2024 03:37
debbie	Hi Brandon, More information is needed. What gene (please include coordinates) and genome are your inquiring about? thanks, debbie

Link to this post | posted 09 Feb, 2024 03:33
debbie	Hi ACMPhageHunters, More information is needed to help. What gene (please provide coordinates) and genome are you inquiring about? Thanks, debbie

Link to this post | posted 08 Feb, 2024 21:00
debbie	Hi Nancy, I used the posted file and was able to auto-annotate. Were other students able to annotate. sounds like something is up with the student's installation of DNA Master, or their Windows is interpreting the file as something it is not. debbie

Link to this post | posted 08 Feb, 2024 18:12
debbie	Hi Nancy, Let's see if I can help. But I need more information. Are you in DNA Master, file -> open -> FastA multiple sequence file? Is the file the .fasta file for your phagesDB phage page. Is it a text file? Can you open a different dna master file and navigate around it. Have you tried more than 1 different fasta file? debbie

Link to this post | posted 07 Feb, 2024 23:42
debbie	Adam, What do you want to call here? What is your evidence? debbie

Link to this post | posted 07 Feb, 2024 12:44
debbie	Adam, In the paper referred to the Bioinformatics Guide, the TAC 'template' is X_XX.Y_YY.Z. You have not provided enough information for me to look quickly, but check the guide to interpret the '_ 'and '.'s If your CCCTTTT lines up correctly for a -1 Frameshift, you have a canonical programmed ribosomal slippage. You also used an abbreviation 'PF' that i don't recognize. debbie

Link to this post | posted 07 Feb, 2024 03:03
debbie	Holly, Sorry for the confusion. The only program to use is DeepTMHMM. And only 1 transmembrane protein needs to be found to call a membrane protein. take care to interpret DeepTMHMM hits because some will identify signal proteins. At this time, we are not calling signal proteins. The Bioinformatics Guide is updated. debbie

Link to this post | posted 07 Feb, 2024 01:10

debbie

Hi all, If you are looking how to instal PDM utils on a M1 mac, you can follow this: On a series of MacBook Pro computers, including most recently an M1 mac, I have only ever done native installs of pdm_utils. The hardest part for me is MySQL – use an installer for 8.0.18 or later (earlier versions had funky bugs related to cascaded deletes!). Choose “Legacy Password” when prompted, and you’ll be forced to choose an 8-character or longer password (no more ‘phage’). Make sure you add mysql to the $PATH (export PATH=’/usr/local/mysql/bin:$PATH’) in your ~/.zshrc file. Other than that, if you use some flavor of conda to manage the third-party dependencies (Aragorn, trnascan-se, mmseqs2, blast) pdm_utils should be just fine. –Christian Gauthier Edited 07 Feb, 2024 01:21

Link to this post | posted 06 Feb, 2024 02:10

debbie

Hi Adam, I know very little about Streptomyces phages annotations. When I first looked at this, I most wanted to respect coding potential. Tonight, I had to ask - after your most recent questions - how many BD phages are there. Looks like there are 43 members in BD, with 11 in subcluster BD2. And additional phages in the rest of the BD subcluster. My impression is that the right arms of BD phages are filled with small genes and lots of HGT. Because of the overlap, I expected to see lots of hits to the c terminus of this overlapping gene. But if you blastn 44271- 43849, you can see what this looks like athere are 2 genes in the region. My take, with the short amount of time I have looked at this, is that there is no clear answer here. If you choose the longer overlap, please be sure to include your justification when you submit. There is definitely more digging to do. And no, this won't get answered until you look at all the data, don't leave out the functional investigation as you continue to consider this. debbie