I have 29 fasta files (.fa as extension) named and stored sequences according to their genes.
(Example: ribosomal protein L1, ribosomal protein L6P/L9E,...)
There were a total of 722 species existed among all those 29 fasta files. Each sequence has their genes and species name labelled in the first line and second line was filled with its sequence.
There will be more than 1 gene sequence for 1 species.
I want to transfer the 722 species from 29 fasta files sorted according to genes into separate 722 files (sort them under species instead of genes).
The name of the species in the parent file were enclosed by square brakets [ ]
.
How can I use for loops to extract the 722 files and name the files according to its sequence name?
Example from Ribosomal Protein L1.fa
:
>gi|103486926|ref|YP_616487.1| 50S ribosomal protein L1 [Sphingopyxis alaskensis RB2256]
MAKLTKKQKALEGKVDAQKLHGVDEAIKLVRELATAKFDETLEIAMNLGVDPRHADQMVRGVVTLPAGTGKDVKVAVFAR
Example from Ribosomal Protein L6PL9E.fa
:
>gi|410479108|ref|YP_006766745.1| ribosomal protein L6P/L9E [Leptospirillum ferriphilum ML-04]
MGFTHTVEFTLPSLIKASIEKQTIITLSSPDKELLGQFAADVRSIRPPEPYKGKGIKYSGEKILRKEGKTGKK
For the first example,
Species name: Sphingopyxis alaskensis RB2256
Gene sequence: MAKLTKKQKALEGKVDAQKLHGVDEAIKLVRELATAKFDETLEIAMNLGVDPRHADQMVRGVVTLPAGTGKDVKVAVFA
I wanted to name the file as Sphingopyxis alaskensis RB2256.fa
and insert all sequences with this species name into this file.
I am using bash shell to do this. I can use grep
to have things done:
grep -A+1 "Sphingopyxis alaskensis RB2256" *.fa >> Sphingopyxis alaskensis RB2256.fa
But I will need to do it 722 times to get my sequences sorted according to species.
Is grep in for loops can be used to simplify the work? Or there are alternative methods to do so?
>
is part of the file. It's a requirement for the fasta format.>
and so your question will be asked of people who are familiar with the file formats you deal with.