How can I split the names and a. acid sequences of certain organisms into columns in the fasta file I have?
My fasta file:
>H_sapiens_ [Homo sapiens]
MANGTADVR
>P_troglodytes_[Pan troglodytes]
MANGTADV
>C_elegans_ [Caenorhabditis elegans]
MSSDSKDQ
I'm trying to get a table like this by pulling codes that are only c.elegance and human.
H_sapiens_ [Homo sapiens]. C_elegans_ [Caenorhabditis elegans]
M M
A S
N S
• 1,333 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
logbio •
•
written
by
logbio •
1 answer
Since you are not interested in main sequence, just grep line with the names:
grep "H_sapiens_" your.fasta > output1.txt
grep "C_elegans_" your.fasta > output2.txt
Then merge two files:
cat output1.txt output2.txt > output_final.txt
• 0 views
•
link
written
by
kashiff007
Log in to answer this question.
More posts like this
-
Converting .faa file into a data frame
written by Samet •
-
Extract certain fasta sequnce from multi fasta file
written by luffy •
-
Controlling for mouse and human cell-cycle aspects of data simultaneously in Seurat
written by Aaron •
-
rentrez::entrez_fetch and web API returning empty FASTA records (R)
written by alopex •
-
Resolving multiple entries to just first one in perl program.
written by kriti.awasthi23 •
-
Using BBDuk to trim degenerate primers off of amplicon reads
written by t89oconnell •
-
Get all hits from BLAST api
written by talalamin •
-
Aligning AFLP primers
written by hbw
-
Extracting promoter regions alignment
written by prosvirov.k •
-
Comparing fasta headers to set of 3 letter codes that is found in the header
written by pasanfernando87 •
cross posted: https://stackoverflow.com/questions/71118224/