problem with adding numerical sequence at the end of line

[garethboy]

Hi,

Anyone has any idea how to get this:

>no_name
TATGCATCGATGCACATATGCTAGTGCGCTAGTGTCGAGGCTAGCTACG
>no_name
GACGTACGTAGCATGCATGCATGCGTAGCTGTAGCTAGC
>no_name
GCTAGCTAGGTAGGTCATGTAGTAGGTGCACTGAGCTAGCTAGCTAGCTAGCAGC
>no_name
GCTAGCATGCTAGCTAGCTAGCACTAGCTAGCTAGCTAGCTAATGCATCATC
>no_name
GCTACGTAGCATGCTAGCGGATCATGCATGCATGCTAGCATCGATGCTAGCATGCAT

become this:

>no_name_1
TATGCATCGATGCACATATGCTAGTGCGCTAGTGTCGAGGCTAGCTACG
>no_name_2
GACGTACGTAGCATGCATGCATGCGTAGCTGTAGCTAGC
>no_name_3
GCTAGCTAGGTAGGTCATGTAGTAGGTGCACTGAGCTAGCTAGCTAGCTAGCAGC
>no_name_4
GCTAGCATGCTAGCTAGCTAGCACTAGCTAGCTAGCTAGCTAATGCATCATC
>no_name_5
GCTACGTAGCATGCTAGCGGATCATGCATGCATGCTAGCATCGATGCTAGCATGCAT

[atcghelix]

Here's one way using Perl. Save the text in a file named numbers.pl (or whatever). Usage would be:

perl numbers.pl --in file_to_change.fasta --out revised_file.fasta

#!/usr/bin/perl

use strict;
use warnings;
use Getopt::Long;

my $inFile;
my $outFile;

GetOptions  ("in=s"      => \$inFile,
             "out=s"      => \$outFile);

if (!$inFile or !$outFile) {
    die "Must supply both infile and outfile as command line arguments.\n";
}

open(my $inFH, "<", $inFile) or die "couldn't open infile for reading.\n";
if (-e $outFile) {
    die "Output file $outFile already exists--aborting so you don't overwrite.\n";
}
open(my $outFH, ">", $outFile) or die "couldn't open outfile for writing.\n";
    
my $counter = 1;
while (my $line = <$inFH>) {
    chomp $line;
    if ($line =~ /^(>.*)/) {
        print $outFH $1 . "_$counter\n";
        $counter++;
    } else {
        print $outFH "$line\n";
    }
}

[Jeremy]

Heres another way: R

library(seqinr)
read.fasta("fastafile.fa")->fa
write.fasta(fa,names=paste(getName(fa),1:5,sep="_"),file.out="fa_new_name.fa")

where you swap '1:5' with '1:n', n being the number of sequences you have.

[garethboy]

Anyone know how to use AWK to do this task?

[garethboy]

Thanks. I am pretty weak in Perl. Do you have any idea using AWK to do this?

[atcghelix]

What version of Awk are you running/what operating system?

[garethboy]

Running is UNIX

[atcghelix]

This work? (It assumes all sequence strings are on a single line)

awk '{if($0 ~ /^>/){print $0"_"(NR+1)/2}else{print $0}}' input.fasta > changed.fasta

[Kennels]

try this

paste - - < input.fa | awk ' { print $1"_"NR"\n"$2 } ' > output.fa

make sure to have spaces between the hyphens for 'paste'

[garethboy]

Thank you everybody. I have done my task. =)