Re^2: renaming 1000's of FASTA files


Do you know where your variables are?
	PerlMonks

Re^2: renaming 1000's of FASTA files

by garyboyd (Acolyte)

on Jul 12, 2011 at 08:50 UTC ( #913862=note: print w/replies, xml )

Need Help??

in reply to Re: renaming 1000's of FASTA files
in thread renaming 1000's of FASTA files

Yes that's correct, in the .txt file there are id's with 4 different endings, so the script would scan the multiple fasta file and generate (in this example) 4 files with headers:

F0Z7V0F01A03EB_210
F0Z7V0F01A03EB_180
F0Z7V0F01A03EB_136
F0Z7V0F01A03EB_362
[download]

and the sequence would be the same for all files, ie the sequence from the fasta file F0Z7V0F01A03EB

Comment on Re^2: renaming 1000's of FASTA files Download Code

Replies are listed 'Best First'.

Re^3: renaming 1000's of FASTA files
by CountZero (Bishop) on Jul 12, 2011 at 13:18 UTC

You keep a hash with key the ID without the ending of "_xxx" and value the ID with the ending. But a hash-element can have only ONE value, so you overwrite the older elements with the newer, unless you take care and store the IDs in an array and then store the arrayref as a value in the hash-element.

CountZero

A program should be light and agile, its subroutines connected like a string of pearls. The spirit and intent of the program should be retained throughout. There should be neither too little or too much, neither needless loops nor useless variables, neither lack of structure nor overwhelming rigidity." - The Tao of Programming, 4.1 - Geoffrey James

[reply]

In Section Seekers of Perl Wisdom

Domain Nodelet^?

www.com | www.net | www.org

Node Status^?

node history
Node Type: note [id://913862]
help

Chatterbox^?

How do I use this? | Other CB clients

Other Users^?

Others rifling through the Monastery: (4)

As of 2022-08-12 05:11 GMT

Sections^?

Information^?

Find Nodes^?

Leftovers^?

Voting Booth^?

No recent polls found

Notices^?