Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • amango
    Member
    • Dec 2009
    • 17

    #1

    extracting predicted gene from scaffold: end position precedes start position

    I am trying to extract sequences for a list of predicted genes from genomic scaffolds. The list of predicted genes with Scaffold IDs, start and end positions, and other info comes from published supplementary data. My script to extract the sequences doesn't work because for some genes, the start position is a larger number than the end position (fourth-to-last and third-to-last columns below). Here is an example (numbers have been changed from original):
    geneID Gene_family Class ScaffoldID start_position end_position Number_of_exons Annotation_status
    CSP1 cs Protein candidate gi|294506227|gb|GL650210.1| 61498 52100 2 intact
    CSP10 cs Protein candidate gi|294507212|gb|GL649715.1| 293074 297989 2 intact
    CSP2 cs Protein candidate gi|294507210|gb|GL650017.1| 234944 236074 2 intact
    CSP3 cs Protein candidate gi|294507295|gb|GL649612.1| 323100 323743 2 intact
    CSP4 cs Protein candidate gi|294506227|gb|GL650210.1| 41911 40888 2 intact
    CSP5 cs Protein candidate gi|294507205|gb|GL649712.1| 274408 272617 2 intact
    I am new to working with annotated genomes. Does it make sense that the some "starts" come after the "ends"? Is this because the ORF for this gene is on the opposite strand of the scaffold? If so, and if I want to obtain that sequence, what's the best way to get it--should I extract the sequence in the scaffold between the two numbers and then find the reverse complement?

    Thanks for any pointers.
  • zhidkov.ilia
    Member
    • Dec 2010
    • 25

    #2
    Some genes transcribed from opposite strand of DNA, resulting in reverse coordinates. You can add additional column (i.e. strand) adding '+' in cases when start_position < end_position and '-' start_position > end_position.

    Comment

    Latest Articles

    Collapse

    • kinpict.com
      Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by kinpict.com
      The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
      Today, 06:37 PM
    • SEQadmin2
      16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by SEQadmin2





      NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
      10-30-2026, 06:46 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 11-12-2026, 07:59 AM
    0 responses
    28 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 11-05-2026, 09:25 AM
    0 responses
    35 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 11-03-2026, 11:54 AM
    2 responses
    56 views
    0 reactions
    Last Post kinpict.com  
    Started by SEQadmin2, 10-28-2026, 09:57 AM
    0 responses
    175 views
    0 reactions
    Last Post SEQadmin2  
    Working...