This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Replace header of fasta file according to another file

Hi

I have a file T.fasta as:

>AAAA
GVNCLPQWERTYI
>BBBB
HJKSDRTYPALVNM
>CCCC
TRHJPKFLMVNCTEGH

and another file.txt with two columns: One column with header of the T.fasta file and a second column with the new name like:

CCCC     C_pr3
BBBB     B_pr2
AAAA     A_pr1

Now I want to replace the header of T.fa with the second column of the second file if first column of second file match with the header of T.fa.

The deisred output should look like this:

>A_pr1
GVNCLPQWERTYI
>B_pr2
HJKSDRTYPALVNM
>C_pr3
TRHJPKFLMVNCTEGH

**My both file are not sorted accordingly.

fasta
• 988 views • link
updated by Ram Ram • written by Nelo • Nelo

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Bash solution to replace part of FASTA headers witch contain matching values in another file
    written by clmattson • clmattson
    Hi all, this is my first time posting on BioStars! I thought I'd easily find a bash one-liner to take care of this task but …
  • Join two files by matching columns
    written by er.doug.ragnar • er.doug.ragnar
    Hi I am trying to join two files, one has 4 columns the other 3. I want to match the third column of the first …
  • Replace ID column of vcf file with another column in the text file
    written by seta seta
    Hi all, I have the standard vcf file (VCFv4.2) and would like to replace the ID column of vcf file with the column in the …
  • How to match fasta header with the first column of a text file and append the second column of the …
    written by timmers • timmers
    Aloha! I have a fasta file that looks like the following: >FFSA34B_100_M7_ID10014 ATCTAACAATGTTGCTCATGCAGGCCCTGCAGTAGATTTAACCATTCTATCCCTTCACCTAGCAGGTGTATCCTCCTTAATAGGAGCCATCAATTTTACAACTACTATTGCTAACAGACGTTTAGAAGGTATACCTACAGAAAAAATACCCTTATTTATT >ART01B_100_M7_ID10021 TGTCTTCTAATATTTCCCATGGAGGGGCCTCTGTTGATATGGCTACTTTTTCTCTCCACCTGGCTGGTGTGTCCTCTATTTTGGGTGCTATCAATTTTATAACTACTATTATTAATATGCGGAGTAATTCCTTGTCTTTTGTGCGGGTTCCTTTGTTTGT >PHR01B_100_M7_ID100236 TAAGAGGAAATGTAGCCCATAGTGGACCTGCAGTAGATTTAGGAATTTTTTCTTTACATTTAGCTGGGGTTTCTTCTATTTTAGGGGCTATTAACTTTATTACTACAATTATTAGTATACGGGAGAATATGGGGATAGAAAAAGTTACTTTATTTAGTTG I have a text file that looks like …
  • MSA ID replacement with a list ID
    written by lessismore lessismore
    Hey all, i have 2 files: 1 is an MSA >tr|B4GH81|B4GH81_MAIZE MCGILAVLGC----------SDCS--QARR-AR----ILACSRR------------------------------------LKHRGPDWSGLYQH------------------------------------------------------------------EGNFLAQQRLAIVSPLSGDQPLFNEDRTVVV-------VANGEIYNHKNVR--KQFT-GAH-- >tr|A0A1BAAANT6|A0A1BAAANT6_SORBI MCGILAVLGC----------SDWS--QARR-AR----VLACARR------------------------------------LKHRGPDWSGLYQH------------------------------------------------------------------EGNFLAQQRLAIVSPLSGDQPLFNEDRTVVV-------VANGEIYNHKNIR--KQFT-GTH--NFTTGSDCEVIIPLYEKYGENFVDMLDGVFAFVLYDTRDRT------YVAARDAIGVNPLYIGWFVVG--------------------LE-GSPDLKAAREVADYLGTIHHEFHFTV----- >tr|K3TTT2J1|K3TTT2J1_SETIT --------------------------------------------------------------------------------LRHRGPDWSGLHCH------------------------------------------------------------------QDCYLAHQRLAIVDPTSGDQLLYNEDKSVVV-------TVNGEIYNHEELK--AKL--TTH--KFQTVSDCEVIAHLYEEYGEEFVDMLDGMFAFVLLDTRDKS------FIAARDAIGICPLYMGWGLDGSVWFSSEMKALSDDCERFITFPPGHLRWYLHIKKG-SGLRRWFNLPWFL-----E--SI-PST-PYNPLLLQGMFEK---------------- It goes like that for some hundred sequences. The …
  • Separate column values for summary
    written by Jack Jack
    So I have 2 files and something peculiar occurred. My origin file is a list of **SELECT** genes like this: A1BG A2M A2MP1 And I …
  • Find and replace
    written by Nano • Nano
    I am working with a group of gene names which I want to convert to their GenBank accession numbers. I have them sorted out already: …
  • merging sequences from bedtools options -name, -tab into one file
    written by ggman ggman
    Hello all, I wrote a script below to analyze two file formats from bedtools (`-tab` option, `-name` option) so it could combine the headers if …
  • filtering reads based on read counts
    written by sumithrasank75 sumithrasank75
    I have parsed out a fastq file and have a text file like ``` read_id read-count aaaa 10 bbbb 1000 ccccc 10000 dddd 1 ``` …
  • matching and merging data frames of different lengths
    written by Diana Diana
    Hi, I have 2 files and I want to match the first column in first file with the first column in second file then if …