Python Script to Calculate Total Number of genes
This is the HCV Genome Sequence:
>GU294484.1 Hepatitis C virus isolate PK-1, complete genome
ACCTGCCTCTTACGAGGCGACACTCCACCATGGATCACTCCCCTGTGAGGAACTACTGTCTTCACGCAGAAAGCGTCTAGCCATGGCGTTAGTATGAGTGTCGTGCAGCCTCCAGGACCCCCCCTCCCGGGAGAGCCATAGTGGTCTGCGGAACCGGTGAGTACACCGGAATTGCCAGGACGACCGGGTCCTTTCTTGGATAAACCCGCTCAATGCCTGGAGATTTGGGCGTGCCCCCGCAAGACTGCTAGCCGAGTAGTGTTGGGTCGCGAAAGGCCTTGTGGTACTGCCTGATAGGGTGCTTGCGAGTGCCCCGGGAGGTCTCGTAGACCGTGCAACATGAGCACACTTCCTAAACCTCAAAGAAAAACCAAAAGAAAACCCATCCGTCGCCCACAGGACGTCAAGTTCCCGGGTGGCGGACAGATCGTTGGTGGAGTATACGTGTTGCCGCGCAGGGGCCCACGATTGGGTGTGCGCGCGACGCGTAAGGCTTCTGAACGGTCACAGCCTCGCGAACGACGACAGCCTATCCCCAAGGCGCGTCGGAGCGAAGGCCGGTCCTGGGCTCAGCCTGGGTACCCTTGGCCCCTCTATGGTAATGAGGGCTGCGGGTGGGCAGGGTGGCTCCTGTCCCCCCGCGGCTCCCGTCCATCTTGGGGCCCAAACGACCCCCGGCGAAGATCCCGCAACTTGGGTAAAGTCATCGATACCCTTACGTGCGGATTCGCCGACCTCATGGGGTACATCCCGCTCGTCGGCGCTCCCGTAGGAGGCGTCGCAAGAGCCCTCGCGCATGGCGTGAGGGCCCTTGAAGACGGGATAAATTTTGCGACAGGGAACTTGCCCGGTTGCTCCTTTTCTATCTTCCTTCTTGCTCTACTCTCTTGCTTAATTCATCCAGCAGCCAGTCTAGAGTGGCGGAATACGTCTGGTCTCTATGTCCTTACCAACGCCCGTTCCAACAGCAGTATAGTGTACGAGGCCGACGACGTTATCCTGCACACACCCGGCTGTATACCTTGTGTTCAGACCGGCAACACATCCAAGTGCTGGACCCCAATGACACCCACGGTGGCAGTTAAGTATGTCGGAGCAACCACCGCTTCGATACGCGGTCATGTGGACCTGTTAGTGGGCGCAGCCACGATGTGTTCTGCGCTCTACGTGGGTGATGTGTGCGGAGCCGTCTTCCTCGTGGGGCAAGCCTTCACGTTCAGGCCGCGACGCCATCAAACGGTCCAGACCTGCAACTGCTCGCTGTACCCAGGCCATCTCACAGGACATCGAATGGCTTGGGATATGATGATGAACTGGTCCCCTGCTGTTGGCATGGTGGTGGCGCACATCTTACGCCTACCCCAGACCCTGTTTGATATAATAGCCGGGGCCCATTGGGGCGTCTTGGCGGGTCTAGCCTACTATACCATGCAGGGCAACTGGGCCAAGGTCGCAATCATCATGGTTATGTTCTCAGGGGTCGATGCCGTTACGTACATCACTGGGGGCACTGCAGCTCGTGGGGGCCAAGGGCTGGCTAGCCTAATCGTCCGGGGGCCTGAGCAGCGCCTGGAGCTGATCAACACCCATGGCTCGTGGCACATCAACAGTACTGTCCTCCACTGCAATGAGTCCATAAACACAGGGTTTATAGCTGGGTTGTTTTATTATCATAAGTTCAACTTACTGGATGTCCCGAAGGCTCAGCAGCTGCAAGCCCATCACTTTCTTCAGGCAGGGGTGGGGCCCCTTGACAGATGCCAACATCCACCGGCCCTTCTGATGACAACCGTACTGCTGGCATACGCACCTAGACCTTGTGACAGCGTAAAGCAGCACGTGTCTCCGGTCCTGTGTATGCTTCCACACCATCGCCCAGTGGTGGTAGGCACTACTGATCCTAAGGGCGCTCCCACCTATAACTGGGGCGAGAATGAGACAGACGTGTTCCTGCTGAATCCCTGCGGCCTCCTAGTGGTCGGTGGTTTGGGTGGCACGTGGGAGGAACTCCACCGGGGTTTGTCAAGACGTGCGGAGGTTCCCCCTTGTGACATCTATGGGGGTGGGGGGGAGATCCACCAATGGTTCAGACCTCTTCTGCCCCACCGACTGCTTCAGGAAACATCCCGAGGCCACATACAGCCGGTGCGGCTCGGGGCCCTGGTTGACACCTCGATGCATGGTCGACTATCCATACCGGCTTTGGCATTACCCATGTACAGTCAATTTTACACTGTTCAAGGTGAGGATGTTTGTGGGTGGGTTTGGCATCGGTTTACCGCCGCTTGCAACTGGACTAGGGGGGAGCGCTGCGATATCGAGGATCGTGACCGCAGCGAGCAACATCCCCTGCTGCATTCAACAACTGAGCTTGCCATACTGCCTTGCTCTTTCACGCCCATGCCCGCATTGTCAACAGGGTTAATACACCTCCACCAAAACATCGTGGATGTCCAATACCTTTATGGCGTTGGATCTGGCATGGTGGGATGGGCGTTGAAATGGGAATTTGTCATCCTCGTTTTCCTCCTCCTAGCAGATGCACGCGTGTGCGTTGCCCTTTGGCTGATGCTGATGATATCACAAGCAGAAGCAGCCTGGAGAACTTGTCACGCTGAACGCCGTCTCTGCTGCCGGGACACATGGTATCGGCTGGTACCTGGTAGCATTTTGCGCGGCGTGGTACGTGCGGGGAAACTCGTCCCGCTGGTGACCTACAGCCTGACGGGTCTTTGGTCCCTAGCATTGCTCGTCCTTCTACTCCCCCAGCGGGCGTATGCTTGGTCGGGTGAAGACAGTGCCACCCTCGGCGCTGGGATCTTGGTCCTCTTCGGCTTCTTTACCCTGTCACCTTGGTATAAGCACTGGATCAGCCGCCTCATGTGGTGGAACCAGTACGCCATATGTAGGTGTGAGTCTGCTCTCCAAGTATGGGTCCCCCCCCTACTTGCCCGCGGGAGTAGGGACGGTGTTATCCTGCTAACAAGCCTGCTTTATCCATCATTAGTTTTTGACATCGCTAAGCTGCTGATAGCCGTAATAGGCCCATTATATCTAATACAGGCCGCCATCACTACTACCCCCTACTTTGTGCGTGCGCATGTTCTGGTCCGCCTTTGCATGTTCGTGCGCTCCGTGACGGGGGGAAAGTACTTCCAGATGGCCATACTGAGCGTCGGCAGATGGTTTAACACCTACCTATATGACCACCTTGCACCGATGCAACACTGGGCCGCAGCAGGCCTCAAAGACCTGGCAGTAGCCACTGAACCTGTAATATTCAGTCCCATGGAAATCAAGGTCATCACTTGGGGCGCGGACACGGCAGCTTGCGGAGATATCCTATGCGGGCTGCCCGTCTCTGCACGATTAGGCCGTGAGGTGTTGTTGGGACCTGCTGATGACTATCGGGAGATGGGCTGGCGTCTGTTGGCCCCGATTACAGCATACGCCCAGCAAACTAGGCGTCTTTTTGGGACTATTGTGACCAGCTTGACTGGCAGGGACAAGAACGTGGTGGCCGGCGAAGTGCAGGTGCTTTCTACGGCTACCCAGACCTTCCTAGGTACAACATTGGGAGGGGTTATGTGGACTGTTTACCATGGAGCAGGTTCGAGAACACTTGCGGGCGTCAAACATCCTGCGCTCCAAATGTACACAAATGTAGATCAGGACCTCGTTGGATGGCCAGCTCCTCCGGGGGCTAAGTCTCTTGAACCGTGCACCTGCGGGTCTGCGGACTTGTACTTGGTTACCCGCGAAGCTGATGTCATCCCTGCTAGACGCAGGGGGGACTCCACAGCGAGCTTGCTCAGTCCTAGGCCTCTCGCCTGTCTCAAAGGTTCCTCTGGAGGTCCTGTTATGTGCCCTTCGGGCCACGTAGCGGGGATCTTTAGGGCTGCTGTGTGCACCAGAGGTGTAGCAAAAGCCCTACAGTTCATACCAGTGGAAACCCTTAGCACACAGGCTAGGTCTCCATCCTTTTCTGACAATTCAACTCCTCCTGCTGTTCCACAGAGCTATCAAGTAGGGTACCTTCATGCCCCGACCGGCAGCGGTAAGAGCACAAAGGTCCCGGCCGCTTATGTAGCACAAGGATATAATGTTCTCGTGTTGAATCCATCAGTGGCGGCCACACTAGGCTTCGGCTCTTTCATGTCGCGAGCTTATGGGATCGACCCCAACATCCGCACCGGGAACGGCACGGTTACAACTGGTGCTAATCTGACCTATTCCACCTATGGTAAGTTTCTCGCGGACGGGGGTTGCTCGGGGGGAGCATATGATGTGATTATCTGTGATGAGTGTCATGCCCAAGACTCTACTAGCATACTGGGTATAGGCACGGTCCTAGATCAGGCTGAAACGGCTGGGGTGAGGCTGACGGTTTTAGCAACAGCAACTCCCCCAGGCAGCATCATTGGGCCCCATTCTAACCTCAAAGAAGTGGCCCTTGGTTCTGAGGGGGAGATCCCTTTCTTCGGCAAGGCCATACCGCTAGCCCTGCTAAGGGGGAAAGGCACCTTATTTTTTTCCATTCCAAGAAAAAATGTGATGAGATGGCATCCAAACTCAGAGGCATGGGGCTCAACGCTGAAGGAGTACTACAGGGGTCTTGATGTGTCCGTCATACCAACATCAGGAGACGTTGTAGTTTGCGCTACTGACGCCCTCATGACTGGATTCACCGGAGACTTCGACTCTGTCATAGATTGCAACGTGGCTGTTGAACAGTACGTTGATTTCAGCTTGGACCCCACCTTTTCCATTGAGACTCGCACTGCTCCCCAAGACGCGGTTTCCCGCAGTCAACGTCGTGGCCGTACGGGCCGAGGTAGACTCGGCACGTACCGATATGTCACCCCCGGTGAAAGACCGTCTGGGATGTTTGACTCGGCTGTTCTCTGTGAGTGCTATGACGCGGGCTGCTCGTGGTACGACTTGCAGCCCGCCGAGACCACAGTCAGACTAAGAGCTTACTTGTCCACGCCGGGGTTACCTGTCTGCCAAGACCACTTGGAATTTTGGGAGAGCGTCTTCACTGGACTAACTCACATAGATGCCCACTTTCTATCACAGACCAAGCAGCAGGGACTCAACTTCCCATACCTAGCTGCCTACCAAGCCACTGTGTGCGCTCGCGCGCAAGCTCCTCCCCCAAGTTGGGACGAGACATGGAAGTGTCTCGTGCGGCTTAAGCCAACACTACATGGACCTACACCCCTTCGATATCGGCCGGGGCCTGTCCAAAATGAAACCTGCTTGACACACCCCATCACAAAATACCTCATGGCATGCATGTCAGCCGATCTGGAAGTAACCACCAGCACCTGGAGCACCTGGGTGTTGCTCGGAGGGGTCCTCGCGGCCCTGGCAGCCTACTGCTTGTCGGTCGGCTGCGTAGTCATTGTGGGCCACATTGAGCTGGGGGGCAAGCCGGCGCTCGTTCCTGACAAAGAAGTGTTGTATCAACAATACGATGAGATGGAGGAGTGCTCACAAGCTGCCCCATATATCGAACAAGCTCAAGTAATAGCCCACCAGTTCAAGGAAAAAGTCCTTGGATTGCTACAGCGAGCTACCCAACAACAAGCTGTCATTGAGCCCATAGTAGTTACCAACTGGCAAAAGCTTGAGGCCTTCTGGCACAAGCACATGTGGAACTTTGTGAGTGGGATTCAGTACCTAGCAGGTCTCTCCACTTTGCCCGGCAACCCCGCTGTGGCGTCTCTTATGGCGTTCGCTGCTTCAGTCACCAGTCCCCTGACGACCAATCAAACTATGTTTTTTAACATACTCGGGGGATGGGTTGCTACTCATTTGGCAGGGCCCCAGAGCTCTTCCGCATTCGTGGTAAGCGGCTTGGCCGGCGCTGCCATAGGGGGCATAGGCCTGGGCAGGGTCTTACTTGACATCCTGGCAGGATACGGAGCTGGTGTCTCAGGCGCCTTGGTGGCTTTCAAAATCATGGGGGGGGAACTCCCCAATGCCGAGGACGTGGTCAATCTGTTGCCCGCCATACTATCTCCGGGTGCTCTCGTCGTCGGGGTGATATGCGCTGCCCTACTACGTCGGCACGTGGGACCTGGGGAGGGAGCGGTACAGTGGATGAACAGGCTCATCGCGTTCGCATCCCGGGGCAACCACGTCTCACCGACGCACTATGTTCCCGAGAGCGATGCTGCGGCAAGGGTCACCGCATTGCTGAGTTCTCTAACTGTCACAAGTCTGCTCCGGCGGTTACACCAGTGGATCAATGAAGACTACCCAAGCCCTTGTAGCGACGATTGGCTACGTACCATCTGGGACTGGGTCTGCATGGTGTTGCTCGACTTCAAGACATGGCTGTCTGCTAAGATCATGCCATTGCTCCCTGGGTTGCCCTTCATTTCCTGTCAAAAGGGATATAAGGGCGTTTGGCAGGGGGACGGCGTGGTGTCCACTCGCTGTCCTTGCGGAGCAGTGATAACCGGTCATGTGAAGAACGGGTCCATGCGGCTTGCAGGACCACGTACATGTGCTAACATGTGGCACGGCACCTTCCCCATCAACGAGTACACCACCGGACCCAGCACACCTTGCCCATCACCCAACTACACTCGTGCACTGTGGCGCGTGGCTGCCAACAGCTACGTCGAAGTGCGACGGGTGGGAGACTTCCACTACATCACGGGGGCCACAGAAGATGAGCTCAAGTGTCCGTGCCAAGTGCCGGCTGCTGAGTTCTTTACTGAAGTGGATGGGGTGAGACTTCACCGTTACGCCCCTCCATGCAGGCCCCTGTTGAGGGATGAGATCACTTTCGTAGTAGGGCTGAATTCTTACGCGATAGGATCCCAACTCCCTTGTGAGCCCGAACCGGACGTCTCTGTGCTGACCTCGATGTTGAGAGACCCTTCCCATATCACCGCCGAGACGGCAGCGCGCCGCCTTGCACGCGGGTCCCCTCCATCAGAGGCAAGCTCATCCGCCAGTCAACTATCGGCTCCATCGTTGAAGGCCACTTGCCAAACGCATAGGCCTCATCCCGACGCGAGCTGGTGGACGCCAACTTGTGTTTGGCGACAAGAGATGGGCAGCAACATCACACGGGTAGAGTCCGAAACAAAGGTTGTGATTCTTGACTCATTCGAACCTCTGAGGGCCGAGACTGATGACACCGAGCTCTCGGTAGCAGCAGAGTGTTTCAAGAAACCTCCCAAGTATCCTCCAGCCCTCCCTATCTGGGCTAGGCCAGACTACAACCCTCCACTGTTGGATCGTTGGAAATCACCGGATTATGAACCACCAATTGTTCATGGGTGCGCCTTACCACCACAGGGTACTCCACCGGTGCCTCCCCCTCGGAGGAAAAGAACAATCCAGCTGGACGGCTCCAATGTGTCCGCGGCGCTAGCTGCGCTAGCGGAAAAATCATTCCCGGCCTCAAAACCGTTGGAAGCGGGTAGCTCATCCTCAGGGGTCGATACACAGTCCAGCACTACTTCCAAGGTGCCTCCCTCTTCGGAGAGAGAGTCCGACACAGAATCGTGCTCGTCCATGCCTCCTCTCGAGGGGGAGCCGGGCGATCCAGACTTGAGTTGCGACTCTTGGTCCACTGTTAGTGACAGCGAGGAGCAGAGCGTGGTCTGCTGCTCTATGTCGTATTCTTGGACCGACGCCCTGATAACACCATGTAGTGCTGAGGGAGAGAACTGCCCATCAGCCCACTCAGCAATCTTGGAGAGACATCACAACCTAATCTATTCAACGTCGTCTAGAATCGCTTCTCAACGTCAGAAGAAGGTCACCTTCGACAGGCTGCAGGTGCTCGACGACCATTACAAAACTGCATTAAAGGAGATAAAGGAGCGAGCGTCAAGGGTAAAGGCTCGCATGCTCACCATCGAGGAAGCGTGCGCGCTCGTCCCTCCTCACTCTGCTCGGTCAAAGTTCGGGTATAGTGCGAAGGACGCTCGCTCCCTGTCCAGCAAGGCCATTAACCAGATCCGCTCCGTCTGGGAGGACTTGCTGGAAGACACCACAACTCCAATTCCAACCACCATCATGGCGAAGAGCGAGGTTTTTTGTGTGGATCCTACTAAAGGAGGCCGTTTTTTGCTCGTCTCATTGTCTACCCTGACCTGGGGGTGCGCATCTGTGAGAACGTGCCCTATATGGCGTGATACAGAAGTGGGGAGTGGGGACGATGGGTCCTGCCTATGGATTCCAATACTCGCCTCAACAGCGGGTCGAACGTCTGCTGAAGATGTGGACCTCAAAGAAAGCCCAGTTGGGGTTCTCGTATGGTACCCGCTGCTTTGGCTCGACTGCCACTGGACAGGACATCAGGGTGGAAGAGGAGATATACCAATGCTGGAGCCTTGGACCGGAGGCCAGGAAAGTGATCTCCTCCCTCACGGAGCGGCTTTACTGCGGAGGCCCTATGTTCAACAGCAAGGGGGCCCAGTGTGGTTATCGCCGTTGCCGTGCCAGTGGAGTTCTGCCTACCAGCTTCGGCAACACGATCACTTGTTACATCAAGGCCACAGCGGCTGCAAAGGCCGCAAACCTCCGGAAGCCTGGCTTCTTGTTTGGGGAGGATGGATCCTGGTCGTATTACCTGAGGACCGAATGGGGTTCGAATGGAGATACGGGCAGTCCTGGGAGAGCCTTCACCGGAGGCTATGGACCAGGTATTTCTTGCTTCCACCCGGAGATGGCCCCCACAGGCCAACCCTACGACCTTTGGGCTCATTACATCTGGCTCCTCCAACGTCTCCGTGGCGCGGGACGATACGGGGAAGAGGTATTATTACCTCACTCGTGGTGCCACCACCCCCCTGGCCCGTGCTGCTTGGGAGACAGCTCGTCACACTCCAGTTACCTCCTGGCTGGGGAACATCATCATGTACGCGCCTACTATTTGGGTGCGCATGGTGGTGGTGGCACACTTTTTCTCCATACTCCAATCCCAGGAGATACTTGGTCGCCCCCTTGGCTTTGGAATGTACGGGGCCACTTACTCTGTCACTCCGCTGGATTTACCAGCAATCATTGGAAGACTCCATGGTCTACGCGCGTTTACGCTCCATATTTACTCTCCAGCAGAGCTCAATACGGTCGCGGGGACACTCAGGAAGCTTGTGCCCCCCCCTACGAGCTTGGAGACATCGGGCACGAGCAGTGCGCGCTATGCTTATCGCCCAGGGAGGGAAGGCCAGGATTTGTGGGCTTTATCACTTCAATTGGGCGGTACGCACCAAGACCACCCTCACTCCACTGCCAGCCGCTGGCCAGTTGGATTTATCCATCTGGTTTACGGTTGGTGTCGGCGGGAACGACATTCTCGCAGCGTGTCACGCGCCCGAACCCGCCATTTGCTGCTTTGCCTACTCCCTACTAACAGTAGGGGTAGGCATCTTTCTCTTGCCAGCTCGATGAGCTGGTAAGATAACACTCCATTTCTTTTTTGTTTTTTTTTTTTTTTTTTTT
I want to write Python script to Calculate Total Genes present in this 9474 bp sequence. The Start Codon is (ATG) and I want to calculate the genes from all 3 reading frames
• 3,739 views
•
link
1 answer
So we can say that finding genes through python script is not possible
• 0 views
•
link
Log in to answer this question.
Sounds like a school project, could you show what you have tried to get a starting point to work on ?
That is not how you utilise BioPython. Take a look at the link I gave you.
Some points :
101010to format your text as code, it is much more readable.txtentension if it is a fasta file ?printa file handle like you did withprint(genome)(you can... but it is not readable)Try something like this :
Then to go forward and recover informations about ORF, follow the link that jrj.healey gave you
i already used that i want to calculate total number of genes not ORF
Then use a proper gene finding tools such as GeneMark or Glimmer.
Writing a python script to identify which ORFs are legitimate genes is not going to be simple.
Alternatively, use an annotation pipeline like PGAP/PROKKA, and then use biopython to count the number of
geneorCDStags that are detected in one of the annotated files.How can i use gene annotation pipeline
https://github.com/ncbi/pgap/wiki/Quick-Start
or
https://github.com/tseemann/prokka#installation
Is this just a 'toy' exercise? If so, this might get you started:
https://rody.blog/2015/11/14/orf-finder-for-fasta-files-with-dna-sequences-using-python/
If not, you should know that gene finding is more subtle than simply finding all
ATG...STOPsequences.i want to calculate total number of genes present in my sequence not ORF
I want to calculate genes from all 3 reading frames
Genes cannot be "calculated". They can be at best "found". The algorithm involved is complex, and cannot be explained over a forum discussion in my opinion.
We know that by now. But what you are asking for is a complex task. Not just a simple script. Please read some literature on gene prediction.
What Algorithm can i used?
Pick a suitable one from this list, then read their paper and replicate their algorithm.
EDIT:
Please note that your past posts on this topic have been closed because this is not a subject matter any forum can help you on - it's a project that you will need to work on (if you wish to not use the tools available already). If you continue to open questions for this topic, I'm afraid you risk a stronger response from the forum.
Sir it is my right to ask any question on this forum until i can get the suitable answer
FYI, you're going to be ignored now by all of the regulars on the site, which means you will receive no further help here.
sorry sir but i can find any suitable answer of my question so thats why i am asking here again and again
We have already told you “suitable answers”.
If you did not get an answer you liked the first time around, why would you assume the 2nd or 3rd time would be any different?
Sir my question is can i calculate total no of genes present in HCV genome through simple python script and in the reply one person gives me a code from which orfs can be found easily but i want genes to be calculated and then some one said that use GTF file , some said use linux to calculate genes , some said use annotation pipelines and someone said its required complex algorithm But my question is that it is possible or not without algorithm?
It is not possible without a more complicated algorithm. This has been explained thoroughly already. Further comments from you about this will be deleted as they are not productive.
Which Algorithm can i used for this project?
Correct, but it is not your right to ask the same question over and over again, all whilst ignoring the correct advice you have been given, because it didn’t answer your exact question (because what you want to do is not sensible).
Hello anasjamshed1994!
We believe that this post does not fit the main topic of this site.
The OPs question is beyond what is reasonably able to be assisted on the forum, and this is the second or third iteration of this same thread.
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!