I have assembled transcript files with thousands of sequences with headers as like:
>TRINITY_DN50_c0_g2_i1 len=1961 path=[0:0-1960]
>TRINITY_DN59_c0_g1_i2 len=1961 path=[0:0-1960]
But, I want to rename them into as like:
>TRINITY_1
>TRINITY_2
Just all sequences will retain with TRINITY adding chronological number. Total number sequences are 40000
4 answers
Many posters think that their problems are unique, but in most cases that's not true. Yours, in particular, is one of most frequently discussed problems. That means that searching for "rename fasta header" from the main page will give you numerous solutions.
https://www.biostars.org/post/search/?query=rename+fasta+header
Use seqkit replace, assuming your file name is trinity.fasta:
seqkit replace trinity.fasta -p "(.+)" -r "TRINITY_{nr}" > trinity.renamed.fasta
Where:
-p "(.+)"is the match pattern to match the whole header text-r "TRINITY_{nr}"is the replacement pattern, where{nr}adds the record number.
See https://bioinf.shenwei.me/seqkit/usage/#replace for more information
R version
library(Biostrings)
fa <- readDNAStringSet('your.fasta')
names(fa) <- paste0('TRINITY_',seq(fa))
writeXStringSet(fa,'your_new.fasta',format='fasta')
•
written
by
benformatics
You can also use SEDA (https://www.sing-group.org/seda). Specifically, you would use the "Rename header" operation first, to keep the "TRINITY" part of the headers using the "Multipart header" rename type.
Then, you would use the "Rename header" again, but this time with "Add prefix/suffix" rename type to add the indexes.
We will soon release a new SEDA version that comes with a CLI.
Log in to answer this question.
-
Edit a Multifasta File Through awk
written by pthom010 •
-
RNA de novo assembly - blasts - KEGG - GO
written by dimitrischat
-
Converting an output de-novo transcriptome assembled with Trinity to a .gff3 file
written by Raito92
-
parse of file multifasta in python
written by erick_rc93 •
-
cd-hit-est total seq value doesn't match number of sequences being compared?
written by c.e.chong •
-
Trinity transcriptome quality assessment
written by cwbenson1993 •
-
Please help with connecting the part of the header in a fasta file with another text file
written by seta
-
Cleaning up trinity assemblies
written by dpearton •
-
FASTA sequence extraction using list of sequence name that are not identical with the headers?
written by seta
-
Trinity assembly validation and statistics
written by giorgiocasaburi
what file format ?
what programming language?
Have you stored the information in the header in a separate location?
In addition to what Mensur said, I would also state that renaming is not recommended because the string carries meaning. You will, for example, not be able to extract the longest isoform per gene from the edited file, and it will make reproducing subsequent analysis harder. Most tools should be able to deal with the Trinity identifiers. Unless a tool definitely does not support them, I would leave them as they are.
Thanks@ Michael I appreciate this suggestions.