Extract first and last column of fasta-header
Hi everyone,
I have a multi-fasta file name multi.fasta with the following structure:
>A 124 B
ATCGTA...
>C 567 D
GTCAG...
My goal is to create a new file, with the new fasta-headers containing only the first and last column. If I use
awk -F" " '/>/ {print $1,$(NF)}' multi-fasta > modified_multi-fasta
This will print the fasta headers with the first and last columns, but won't print the nucleotide sequences. Can you guys please help me out?
• 2,073 views
•
link
updated
by
Hugo
•
written
by
genomes_and_MGEs •
•
written
by
genomes_and_MGEs •
3 answers
awk '/^>/ {print $1,$(NF);next;} {print;}'
• 0 views
•
link
written
by
Pierre Lindenbaum
You can use the "Rename header / Multipart header" operation of SEDA (https://www.sing-group.org/seda/manual/operations.html#multipart-header), it is very useful for this kind of FASTA headers.
Log in to answer this question.
More posts like this
-
Append assembly accession to nucleotide accession number in RefSeq Genbank file
written by genomes_and_MGEs •
-
Sequence length of fasta sequences with special characters removed
written by genomes_and_MGEs •
-
Split multi-fasta file and keep structure
written by genomes_and_MGEs •
-
Problem renaming headers in multi fasta file using certain scripts
written by lukhanyomakhabane •
-
Append fasta header to corresponding fasta filename
written by genomes_and_MGEs •
-
Batch rename protein fasta headers
written by genomes_and_MGEs •
-
Add strain name directly on fasta header
written by genomes_and_MGEs •
-
Extract fasta-headers that have the same word
written by genomes_and_MGEs •
-
Removing everything after the last underscore in the header of a fastafile
written by timmers •
-
Tool: Retrieve a subset of FASTA from large Illumina multi-FASTA file
written by hcwang •
