Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • mcrawford
    Member
    • Jun 2015
    • 10

    #1

    Need my reference sequence to be in a specific format

    I need my reference genome (S. cerevisiae strain S288c) to be in a specific format for some analysis but can't find anything matching what I need online.

    i.e. the data I have is in this format:
    >I [org=Saccharomyces cerevisiae] [strain=S288C] [chromosome=I]
    CCACACCACACCCACACACCCAC...etc

    But I need it to be in a table like:

    chr pos seq
    1 1 C
    1 2 C
    1 3 A
    1 4 C
    etc.

    I'm hoping that one of the following will be true, please let me know if you can help... Thanks!
    1. The data already exists somewhere online
    2. There's a program that can make this file for me out of my reference sequence (possibly PySAMstats?)
    3. Some hints to write my own R script to make the file?
    Last edited by mcrawford; 09-22-2015, 07:53 AM.
  • GenoMax
    Senior Member
    • Feb 2008
    • 7142

    #2
    There will always be cooler ways of doing this. What follows happens to be one way. I am assuming chromosome 1 is followed by 2 etc in the genome file.
    Use at your own risk

    Save the code below in a file (e.g. table.pl).

    Code:
    use strict;
    use warnings;
    
    my $datafile = $ARGV[0];
    my $outfile = $ARGV[1];
    
    open (IN, $datafile) or die "can't open the datafile: $datafile\n";
    open (OUT, ">$outfile") or die "can't open the outputfile: $outfile\n";
    
    my $count = 1;
    my $chrom = 0;
    
    while(my $line = <IN>){
            if ($line =~ /^>/){
                    $chrom++;
                    $count = 1;
                    }
            else {
                    chomp $line;
                    my @nuc = split(//, $line);
                    while (@nuc){
                            my $a = pop(@nuc);
                            print OUT $chrom." ".$count." ".$a."\n";
                            $count++;
                            }
                    }
    
    }
    close IN;
    close OUT;
    Run like this

    Code:
    $ perl table.pl genome.fa table.txt

    Comment

    • mcrawford
      Member
      • Jun 2015
      • 10

      #3
      Originally posted by GenoMax View Post
      There will always be cooler ways of doing this. What follows happens to be one way. I am assuming chromosome 1 is followed by 2 etc in the genome file.
      Use at your own risk
      Thanks so much GenoMax, this worked perfectly and produced exactly what I wanted! You're the best

      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
      27 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...