Beefy Boxes and Bandwidth Generously Provided by pair Networks
Think about Loose Coupling
 
PerlMonks  

Re: Renaming headers works for some files but fails for some files

by bioinformatics (Friar)
on May 09, 2013 at 20:16 UTC ( [id://1032838]=note: print w/replies, xml ) Need Help??


in reply to Renaming headers works for some files but fails for some files

To save yourself some grief, I would simply look for > and then replace the line, regardless of what else is there (which is what you are doing anyway). You can also use  print $variable and sleep(3); in parts of the code to help see what is present in certain places, particularly if you don't get a match. We don't have the input files, so we can't tell if there are spaces or tabs in the FASTA headers.

use strict; use warnings; my @files = glob("*.fa"); for my $file ( @files ) { # open the input file open my $in, "<", "$file" or die "Cannot open "$file": $!\n"; # open the output file open my $out, ">", "$file.out" or die "Cannot open "$file.out": $!\ +n"; # reset contig number my $contig_number = 1; while ( <$in> ) { chomp; if ( $_ ~= m/^>/ ) { # it's a header print $out ">config_$contig_number\n\n"; $contig_number++; } else { # it's sequence print $out "$_\n\n"; } close $in; close $out; } } my $file_count = @files; print "Successfully processed $file_count files!\n"; exit;
Bioinformatics

Replies are listed 'Best First'.
Re^2: Renaming headers works for some files but fails for some files
by MVRS (Acolyte) on May 10, 2013 at 03:13 UTC
    Thank you very much i made changes accordingly to my script it works.......... may i know your mail id please

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://1032838]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others imbibing at the Monastery: (3)
As of 2024-04-26 04:41 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found