Entering edit mode
7.2 years ago
shashank.gupta
▴
10
Hi, I am trying to change the fasta header with a given list. I have a fasta file which looks like-
>P1212_1 M04I6TDD
GACGAACGCTGGCGGCGTGCTTAACACATGCAAGTCGAACGAGAGGAGAAGGAAGCTTGCCTCTTTCAAATCT
>P1213_2 M04JTWVX
GATAAACGCTGGCGGCGCACATAAGACATGCAAGTCGAACGAACTTAACCATTAGTTTACTATTGGAGCGGTTAGTG
And I have another file in text format (tab delimited), which looks like -
P1212 SHA1111
P1213 SHA2222
P1214 SHA3333
I want to rename the header of above fasta file according to corresponding name, the output file looks like-
>SHA1111_1 M04I6TDD
GACGAACGCTGGCGGCGTGCTTAACACATGCAAGTCGAACGAGAGGAGAAGGAAGCTTGCCTCTTTCAAATCT
>SHA2222_2 M04JTWVX
GATAAACGCTGGCGGCGCACATAAGACATGCAAGTCGAACGAACTTAACCATTAGTTTACTATTGGAGCGGTTAGTG
I am trying to use scripts mentioned here Renaming fasta headers according to a matching name list , but the output file is not correct. I'm so grateful for helping. Regards.