Thank you!!!!!. It works.
• 1 views
•
link
written
by
kari_vo3 •
How To Rename FASTA Headers
Hello, i want to change the fasta header of this input file :
>M04631:312:000000000-C6V6K:1:2107:11495:1734 1:N:0:ACTGAGCG+TTATGCGA
>M04631:312:000000000-C6V6K:1:2107:13059:1785 1:N:0:ACTGAGCG+TTATGCGA
In this fasta header:
>adjH001
>adjH002
>adjH....
>adjH099
>adjH100
what script do I have to use??
• 3,505 views
•
link
updated
by
5heikki
•
written
by
kari_vo3 •
•
written
by
kari_vo3 •
1 answer
If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
• 1 views
•
link
written
by
WouterDeCoster
Log in to answer this question.
More posts like this
-
Using shell command in Nextflow DSL=2 for BWA mem aligner
written by rohitsatyam102
-
Word size CD-HIT-EST
written by Nathan •
-
How to get script for inserting information in a table to work?
written by Mimmi Ahlmén •
-
Extract Reverse Reads Based on Header Information of Forward Reads
written by dennis.guettler •
-
Labelling sequences within fasta files according to sample name.
written by Mitra •
-
PROVEAN input script
written by s_herrera •
-
Editing header of a fasta file
written by DVR •
-
EDIT:SOLVED:changing the fasta header in a multi-fasta file with biopython.
written by skbrimer
-
Exclude organisms in remoteblast with bioperl
written by lizabe •
-
Help Needed For Formatting Fasta Headers
written by redspider19800915 •
Hello kari_vo3!
Questions similar to yours can already be found at:
We have closed your question to allow us to keep similar content in the same thread.
If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.
Cheers!
Looks like you have MiSeq data, please provide the actual context about what you're doing. It should include the words "MiSeq" and "fastq". It's highly unlikely that your actual question will end up actually renaming fasta entries.
You need to explain your naming convention further. There’s no information in those headers to generate what you’re asking for.
If you just want all the sequences named
adjHfollowed by some number, how high should these numbers go? How many sequences do you have? Since you want them 0-padded this matters.Are the names arbitrary in order?
i have 300000 reads and then i'll renamed them to a metagenomic analysis :
That's not clear. We still don't know where the
adjH001identifiers are based on.You need to tell us how:
>adjH001pairs with>M04631:312:000000000-C6V6K:1:2107:11495:1734 1:N:0:ACTGAGCG+TTATGCGAor any other sequence. Is it purely based on the order in which they occur?
You may be interested in SEDA (https://www.sing-group.org/seda/), a user-friendly GUI-based tool with a powerful option to rename FASTA headers (https://www.sing-group.org/seda/manual/operations.html#rename-header).