Hi everyone,

I'm having trouble figuring this one out. I have ~100 *.fa files with multiple lines of fasta sequences like this: file1.fa

>xyzsequence
atcatgcacac......
ataccgagagg.....
atataccagag.....
>abcsequence
atgagatatat.....
acacacggd.....
atcgaacac....
agttccagat....

The name of each sequence is delimited by a ">" and followed by a newline. I'm trying to figure out how iterate through all of my files with a ".fa" extension and create a single tab-delimited table with the name of the sequence (tab) and the name of the file it came from. Like so:
xyzsequence file1
abcsequence file1
somsequence file2
etc...

Can anyone point me in the right direction?
Many thanks,

To keep the forums high quality for all users, please take the time to format your posts correctly.

First of all, use Code Tags when you post any code or data samples so others can easily read your code. You can easily do this by highlighting your code and then clicking on the # in the editing menu. (You can also type code tags [code] and [/code] by hand.)

Second, avoid adding color or different fonts and font size to your posts. Selective use of color to highlight a single word or phrase can be useful at times, but using color, in general, makes the forums harder to read, especially bright colors like red.

Third, be careful when you cut-and-paste, edit any odd characters and make sure all links are working property.

Thank You.

The UNIX and Linux Forums

*********************************************************

A bit lengthy with a 2nd example file:

$> ll
insgesamt 20
drwxr-xr-x 2 root root  4096 2009-08-06 09:36 .
drwxr-xr-x 3 isau users 4096 2009-08-06 09:32 ..
-rw-r--r-- 1 root root   105 2009-08-06 09:20 file1
-rw-r--r-- 1 root root   103 2009-08-06 09:36 file2
$> cat file1
>xyzsequence
atcatgcacac
ataccgagagg
atataccagag
>abcsequence
atgagatatat
acacacggd
atcgaacac
agttccagat
$> cat file2
>bbbbbbbbb
atcatgcacac
ataccgagagg
atataccagag
>ccccccccccc
atgagatatat
acacacggd
atcgaacac
agttccagat
$> for FILE in file*; do tr -d "\n" < "$FILE"| awk -v file=$FILE '$0{print $0,file}' RS=">" >> outfile; done
$> cat outfile
xyzsequenceatcatgcacacataccgagaggatataccagag file1
abcsequenceatgagatatatacacacggdatcgaacacagttccagat file1
bbbbbbbbbatcatgcacacataccgagaggatataccagag file2
cccccccccccatgagatatatacacacggdatcgaacacagttccagat file2

Hi Zaxxon,

Thanks a million! I didn't want the actual sequence, just the sequence name, so I used some of your code and bits of other things that I pieced together. This is hideous and long (I know ) but it works. Next week I'll try to learn to pipe.

Code:

grep '^>' *.fa >new; sed -e 's/.fa:>/\t/g' new > new2; perl -e ' @cols=(1, 0); while(<>) { s/\r?\n//; @F=split /\t/, $_; print join("\t", @F[@cols]), "\n" } warn "\nChose columns ", join(", ", @cols), " for $. lines\n\n" ' new2 > new3; rm new; rm new2