Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • coding regions

    Hello,

    how can I identify correct coding regions in a sequence--with the help of blast

  • #2
    With one of the translated blast searches.

    From http://www.ncbi.nlm.nih.gov/books/NBK1734/

    There are three varieties of translated BLAST search; “tblastn,” “blastx,” and “tblastx.” In the first variant, “tblastn,” a protein sequence query is compared to the six-frame translations of the sequences in a nucleotide database. In the second variant, “blastx,” a nucleotide sequence query is translated in six reading frames, and the resulting six-protein sequences are compared, in turn, to those in a protein sequence database. In the third variant, “tblastx,” both the “query” and database “subject” nucleotide sequences are translated in six reading frames, after which 36 (6 × 6) protein “blastp” comparisons are made. Protein sequences are better conserved than their corresponding nucleotide sequences. Because the translated searches make their comparisons at the level of protein sequences, they are more sensitive than direct nucleotide sequence searches. A common use of the “tblastn” and “blastx” programs is to help annotate coding regions on a nucleotide sequence; they are also useful in detecting frame-shifts in these coding regions. The “tblastx” program provides a sensitive way to compare transcripts to genomic sequences without the knowledge of any protein translation, however, it is very computationally intensive.

    Comment


    • #3
      Ok. so in annotation, how do I know whether my region is at 5' or 3' end? partial or complete?

      Comment


      • #4
        Co-ordinates of the blast "hit" should provide that information. Unless I am not understanding your question.

        Comment


        • #5
          In case like the one below, is the CDS both 5' and 3' partial or its only 5' partial? and what is the correct start position?

          Range 1: 673 to 1681GenBankGraphics
          Next Match
          Previous Match
          Alignment statistics for match #1 Score Expect Identities Gaps Strand
          1857 bits(1005) 0.0 1008/1009(99%) 1/1009(0%) Plus/Minus

          Query 12 AGTTCCCCCTTATCATTTTTCACATAAAATACTCGTTCCATCAGTCCACGTTCCAGGTTC 71
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1681 AGTTCCCCCTTATCATTTTTCACATAAAATACTCGTTCCATCAGTCCACGTTCCAGGTTC 1622

          Query 72 ACCGCAGAGTTGTTGTGGACCCCAAATTCATACAGGTTGCCCAAACCTGCAATCCGGAAC 131
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1621 ACCGCAGAGTTGTTGTGGACCCCAAATTCATACAGGTTGCCCAAACCTGCAATCCGGAAC 1562

          Query 132 ACTCTGCGGGTTTTGTGTCTCGTGTCTTTAAAGTGGCTACGGAGTTTGGGGTGCCGCGCA 191
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1561 ACTCTGCGGGTTTTGTGTCTCGTGTCTTTAAAGTGGCTACGGAGTTTGGGGTGCCGCGCA 1502

          Query 192 CGTGTGCGCCGAGACTGCACCCCAAGCCACTCTTCAAGACACCCCTATTTCAGCTCCTGG 251
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1501 CGTGTGCGCCGAGACTGCACCCCAAGCCACTCTTCAAGACACCCCTATTTCAGCTCCTGG 1442

          Query 252 AATAGGGCTGGGTCATCCAGCCCACCCAGGTGCCTAACCCAATTACGCCATGCTTTTCCG 311
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1441 AATAGGGCTGGGTCATCCAGCCCACCCAGGTGCCTAACCCAATTACGCCATGCTTTTCCG 1382

          Query 312 CTCAGTGGGTTGTTGAATACTCGCCGTAGGGGTGTCTCCTCCACAAATTGTGCGCGATAT 371
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1381 CTCAGTGGGTTGTTGAATACTCGCCGTAGGGGTGTCTCCTCCACAAATTGTGCGCGATAT 1322

          Query 372 TGGCGACAGATCTCAATGGATTCATCGCGCAGTGGAATTGAAGCTCTTCCAACATCAGTA 431
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1321 TGGCGACAGATCTCAATGGATTCATCGCGCAGTGGAATTGAAGCTCTTCCAACATCAGTA 1262

          Query 432 CGGCACGGCATAAACACATATGGCAGGGCCATCGAAATGGCAGTACGGGAGTCAGTTGGA 491
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1261 CGGCACGGCATAAACACATATGGCAGGGCCATCGAAATGGCAGTACGGGAGTCAGTTGGA 1202

          Query 492 TTCATATTTGCTTGCTCACACTTCCGTACCAGGAACCGCCAAACCGTGAGTTCGTTGGCA 551
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1201 TTCATATTTGCTTGCTCACACTTCCGTACCAGGAACCGCCAAACCGTGAGTTCGTTGGCA 1142

          Query 552 GTGGTTTTCTTTGGAATCCCGCCGAATTCTAGTTTTGCGAGAGTCACTAGGTTTCTGACG 611
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1141 GTGGTTTTCTTTGGAATCCCGCCGAATTCTAGTTTTGCGAGAGTCACTAGGTTTCTGACG 1082

          Query 612 GCAGGGCCAAACTCACGGATGAGAAATGTGGACTGCTCTTTTTGTGTCACTGTGCCTTCT 671
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1081 GCAGGGCCAAACTCACGGATGAGAAATGTGGACTGCTCTTTTTGTGTCACTGTGCCTTCT 1022

          Query 672 TCGTCCACAATCTCTGGAACCAGAACGTGTTTAGTCAGACAGTCCTGAAGGGATAGAATA 731
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 1021 TCGTCCACAATCTCTGGAACCAGAACGTGTTTAGTCAGACAGTCCTGAAGGGATAGAATA 962

          Query 732 GGGGCCGTCGGACCATTGAGTTTCATAACC-TCATCGCTCCATTCTGTTCCATAGCCTGT 790
          |||||||||||||||||||||||||||||| |||||||||||||||||||||||||||||
          Sbjct 961 GGGGCCGTCGGACCATTGAGTTTCATAACCATCATCGCTCCATTCTGTTCCATAGCCTGT 902

          Query 791 TCCACCTGGGAGATTTGTGACTGTTCCCGAGTTGGGAGTTCCTGATCGGATTTTAGCTCC 850
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 901 TCCACCTGGGAGATTTGTGACTGTTCCCGAGTTGGGAGTTCCTGATCGGATTTTAGCTCC 842

          Query 851 CACCTGTAGACATACTCTTTCAATGCCTTCCGGGCGAGGACACATGAAAAATATGCCCCC 910
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 841 CACCTGTAGACATACTCTTTCAATGCCTTCCGGGCGAGGACACATGAAAAATATGCCCCC 782

          Query 911 GCGGCCAGCAGGGAGAAGCGAACTGCTAGCCGAACGTACCGTGACGTTTGTTGAATCACT 970
          ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 781 GCGGCCAGCAGGGAGAAGCGAACTGCTAGCCGAACGTACCGTGACGTTTGTTGAATCACT 722

          Query 971 GAGCCCCCAGTCTGAGTCAACTTGTTTCGAGTAGACATCCACAGACCAT 1019
          |||||||||||||||||||||||||||||||||||||||||||||||||
          Sbjct 721 GAGCCCCCAGTCTGAGTCAACTTGTTTCGAGTAGACATCCACAGACCAT 673

          Comment


          • #6
            There is not enough context here as to what this alignment is about/from. If you are interested in "coding" regions you should be using a translated search for confirmation.

            Comment


            • #7
              in the above case, I have done a blastn so as to detect any hits (several are obtained) but showed only one. my aim in this is first to identify the orientation of the alignment, whether its partial or complete at 5' or 3' or both. Then after that identify reading frames using ORFfinder.

              Comment


              • #8
                Make a note of the co-ordinates in the blast result above. Your query is going from 12-1059 but the hit is actually showing co-ordinates that decrease from 1681-673. That is an indication that the "hit" in this case is on the opposite strand (also indicated by the Q/S = "Plus/Minus" in the header). Depending how how long your query was you can make a judgement of how complete the coverage on the hit is.

                Comment


                • #9
                  Hi,

                  I believe you would get faster and easier-to-read results if you directly go with blastx...

                  s.

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Techniques and Challenges in Conservation Genomics
                    by seqadmin



                    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                    Avian Conservation
                    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                    03-08-2024, 10:41 AM
                  • seqadmin
                    The Impact of AI in Genomic Medicine
                    by seqadmin



                    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                    02-26-2024, 02:07 PM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, 03-14-2024, 06:13 AM
                  0 responses
                  33 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-08-2024, 08:03 AM
                  0 responses
                  72 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-07-2024, 08:13 AM
                  0 responses
                  81 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-06-2024, 09:51 AM
                  0 responses
                  68 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X