Trim fasta sequence names

Dear all,

I have bunch of fasta files with following sequence names format:

>gi|113476179|ref|YP_722240.1| hypothetical protein Tery_2568 [Trichodesmium erythraeum IMS101]

Do you have any script, which would help me to change it to something like this:

>Trichodesmium erythraeum IMS101

Thank you in advance!

Petr

fasta

• 2,740 views • link

updated 4.0 years ago by Ram

• written 10.8 years ago by dvorikus •

1 answer

cat input.txt | sed 's/.*\[/>/' | sed 's/\]//'

• 0 views • link

updated 4.0 years ago by Ram

• written 10.8 years ago by Biomonika (Noolean)

It did exactly, what I needed. Thank you Noolean!

• 0 views • link

written 10.7 years ago by dvorikus •

Glad to help! Please also upvote my answer by clicking on the thumbs up green icon:)

• 0 views • link

written 10.7 years ago by Biomonika (Noolean)

More posts like this

How can I merge GFF files together to produce a file with gene functions from both?
written 4.5 years ago by James •

Hello all, I currently have a GFF file with functional annotations for predicted genes in a metagenome. I have used another tool to predict functions …
Fetching NCBI gene symbols from NCBI protein ids or GI identifiers
written 8.8 years ago by lakhujanivijay

I have the following information from blastx annotations of bacterial genes predicted by prodigal: Sequence name Sequence desc. Sequence length Hit desc. Hit ACC gene_1_contig_1 …
Populate data frame column X with substring from column Y using R
written 9.0 years ago by ddowlin

Hi all, I have a data frame in R with a column of gene identifiers (taken from fasta headers). These vary and include both Ensembl …
What is "^A" in NR database?
written 9.7 years ago by felix.kuo.1211 •

I have download the NR database from http://www.matrixscience.com/help/seq_db_setup_nr_gi.html And I less the file I have download. >gi|66816243|ref|XP_642131.1| hypothetical protein DDB_G0277827 [Dictyostelium discoideum AX4]^Agi|1705556|sp|P54670.1|CAF1_DICDI RecName: Full=Calfumirin-1; …
Fasta header, search and replace...?
written 10.0 years ago by Buffo

Hi everybody, I did a blastp (about 7,000 sequences), after that, I want to make a table of content, some sequences has the same hit …
Get Fasta -Protein Sequences for Given protein names(Homo sapiens)
written 11.6 years ago by genie66 •

I have a list of protein names and I wanted to get the fasta sequence of all the proteins in the list. In batch entrez …
How to safely split a fasta file with concatenated multiple fasta sequences
written 12.2 years ago by akhst7 •

Hi I have a fasta file which contains concatenated multiple ref human virus sequences from NCBI and each sequence has a usual NCBI header (starts …
How To Add An Attribute To A Sequence Header (Of Fasta) Based On A Match Off Another Fasta?
written 12.8 years ago by microbeatic

<p>I have a pair of fasta file. Both of them have common sequences and sequences that are not common. In the sequence header of sequences …
Gi Number To Taxid To Phylum Name
written 12.9 years ago by hdy

<p>I have a bunch of fasta files from NCBI, and I know to use gi_taxid_nucl.dmp to change gi into taxid numbers. But how to turn …
Multiline Fasta To Single Line Fasta
written 15.1 years ago by Palu •

I have a fasta file with following format ``` >gi|321257144|ref|XP_003193485.1| flap endonuclease [Cryptococcus gattii WM276] MGIKGLTGLLSENAPKCMKDHEMKTLFGRKVAIDASMSIYQFLIAVRQQDGQMLMNESGDVTSHLMGFFY RTIRMVDHGIKPCYIFDGKPPELKGSVLAKRFARREEAKEGEEEAKETGTAEDVDKLARRQVRVTREHNE ECKKLLSLMGIPVVTAPGEAEAQCAELARAGKVYAAGSEDMDTLTFHSPILLRHLTFSEAKKMPISEIHL DVALRDLEMSMDQFIELCILLGCDYLEPCKGIGPKTALKLMREHGTLGKVVEHIRGKMAEKAEEIKAAAD EEAEAEAEAEKYDSDPENEEGGETMINSDGEEVPAPSKPKSPKKKAPAKKKKIASSGMQIPEFWPWEEAK QLFLKPDVVNGDDLVLEWKQPDTEGLVEFLCRDKGFNEDRVRAGAAKLSKMLAAKQQGRLDGFFTVKPKE PAAKDAGKGKGKDTKGEKRKAEEKGAAKKKTKK >gi|321473340|gb|EFX84308.1| hypothetical protein …