Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • bioman1
    Member
    • May 2012
    • 80

    GFF editing-Perl, awk

    I would like to get the help of computer experts, in editing my GFF file

    Code:
    ##fsdfsd
    ##sdd-ver gen 5.5.7
    Xm Gen CDS     1       148     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     149     193     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     194     279     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     280     412     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     413     499     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     500     702     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen extracted region        1       148     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="437225 <- 437372";ID=Bm
    Xm Gen extracted region        149     193     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="436969 <- 437013";ID=Bm
    Xm Gen extracted region        194     279     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="435418 <- 435503";ID=Bm
    Xm Gen extracted region        280     412     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="435209 <- 435341";ID=Bm
    Xm Gen extracted region        413     499     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="434376 <- 434462";ID=Bm
    Xm Gen extracted region        500     702     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="434084 <- 434286";ID=Bm
    Replace (Xm Gen CDS) rows with values present in (Xm Gen extracted region) rows. i.e first row ($4 column:1 is replaced with 437225 value and $5 column:148 is replaced with 437372, in row 2 ($4 column:149 replaced with 436969, $5 column:193 replaced with 437013 and so on) and print output like below

    Code:
    ##gff-version 2
    ##source-version geneious 5.5.7
    Xm Gen CDS     437225       437372     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     436969     437013     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     435418     435503     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     435209     435341     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     434376     434462     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen CDS     434084     434286     .       +       .       Name=;created by=User;modified by=User;ID=Bm
    Xm Gen extracted region        1       148     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="437225 <- 437372";ID=Bm
    Xm Gen extracted region        149     193     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="436969 <- 437013";ID=Bm
    Xm Gen extracted region        194     279     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="435418 <- 435503";ID=Bm
    Xm Gen extracted region        280     412     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="435209 <- 435341";ID=Bm
    Xm Gen extracted region        413     499     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="434376 <- 434462";ID=Bm
    Xm Gen extracted region        500     702     .       +       .       Name=Extracted region from gi|371442828|gb|JH557032.1|;Extracted interval="434084 <- 434286";ID=Bm

    I have used grep & sed commands as below.

    head -2 file
    join <(grep "Xm Gen CDS" file | cat -n) \
    <(grep "Xm Gen extracted region" file | cat -n) | \
    sed 's/^[0-9]* //;s/CDS [0-9]*\s[0-9]*\(\s.*interval="\([0-9]*\)\s<-\s\([0-9]*\)\)/CDS\t\2\t\3\t\1/;s/ Xm Gen extracted.*//'
    grep "Xm Gen extracted region" file


    But I would like convert this to shell script like get input from the user and give as desired by the user.
    Any perl script is appreciated.

Latest Articles

Collapse

  • SEQadmin2
    From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
    by SEQadmin2


    Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


    The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
    ...
    06-02-2026, 10:05 AM
  • SEQadmin2
    Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends
    by SEQadmin2


    With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.


    Introduction

    Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
    05-22-2026, 06:42 AM
  • SEQadmin2
    Environmental Genomics in the Age of NGS: From Microbes to Conservation Strategies
    by SEQadmin2

    Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.


    Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...
    05-06-2026, 09:04 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 06-02-2026, 12:03 PM
0 responses
20 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 11:40 AM
0 responses
14 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 05-28-2026, 11:40 AM
0 responses
29 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 05-26-2026, 10:12 AM
0 responses
31 views
0 reactions
Last Post SEQadmin2  
Working...