Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • DCseq
    Junior Member
    • Jul 2017
    • 6

    Script to make sequences equal length

    Hi!
    I would like to perform a motif analysis with my ChIP-Seq data on meme-suite.
    The program requires me to put in sequences of equal length.
    I have the nucleotide sequences from the regions where our protein of interest bound to in ChIP-Seq. These regions are of varying length.
    I would like to make them of equal length by figuring out the sequence of greatest length and adding 'N's to both ends of the remaining sequences.
    Anyone done this/has a R script?
    Many thanks for your feedback in advance!
  • Richard Finney
    Senior Member
    • Feb 2009
    • 701

    #2
    This does not "add to both ends", I think I need a little more information to do that.

    Here's a solution that appends "N"s.
    It's not R, it's an awk "one liner" (a long one liner) :

    cat yourseqfile | awk '{ a[count++] = $0; if (length($0)>l)l=length($0) } END { for (i=0; i<count; i++) { printf("%s",a[i]); for (j=length(a[i]); j<l; j++) printf "N"; printf("\n"); } }'

    example:
    bash-3.2$ cat yourseqfile
    START

    A
    ACGTACGTACGTACGTACGTACGT
    ACGTACGTACGTAC
    ACGT
    DONE

    bash-3.2$ cat yourseqfile | awk '{ a[count++] = $0; if (length($0)>l)l=length($0) } END { for (i=0; i<count; i++) { printf("%s",a[i]); for (j=length(a[i]); j<l; j++) printf "N"; printf("\n"); } }'
    STARTNNNNNNNNNNNNNNNNNNN
    NNNNNNNNNNNNNNNNNNNNNNNN
    ANNNNNNNNNNNNNNNNNNNNNNN
    ACGTACGTACGTACGTACGTACGT
    ACGTACGTACGTACNNNNNNNNNN
    ACGTNNNNNNNNNNNNNNNNNNNN
    DONENNNNNNNNNNNNNNNNNNNN
    Last edited by Richard Finney; 10-05-2018, 06:18 AM.

    Comment

    Latest Articles

    Collapse

    • mylaser
      Reply to Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
      by mylaser
      Kheloyar – Everything You Need to Know About Kheloyaar Login and Kheoyar Id
      If you are looking for an online gaming platform that offers a user-friendly experience, Kheloyar has become a name that many users search for. Whether you're interested in creating a new account, accessing your dashboard through Kheloyaar Login, or learning how to obtain a Kheoyar Id, understanding the platform's features and account process is essential.
      This guide explains everything you need to know about...
      Today, 01:13 AM
    • SEQadmin2
      Advanced Sequencing Platforms Tackle Neuroscience’s Toughest Genomics Problems
      by SEQadmin2



      Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
      ...
      07-09-2026, 11:10 AM
    • SEQadmin2
      Cancer Drug Resistance: The Lingering Barrier to Rising Survival
      by SEQadmin2



      Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

      There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
      07-08-2026, 05:17 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 07-09-2026, 10:04 AM
    0 responses
    17 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 07-08-2026, 10:08 AM
    0 responses
    10 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 07-07-2026, 11:05 AM
    0 responses
    22 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 07-02-2026, 11:08 AM
    0 responses
    31 views
    0 reactions
    Last Post SEQadmin2  
    Working...