This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to replace # with character in fasta file after header and space

Hi, I want to replace for multiple file # with character/name of the strain/gene name and i want to remove this line from my file -1 # ID=1_2660;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5bp;gc_cont=0.687.

>CP077971.1_2660 # 2813973 # 2814887 # -1 # ID=1_2660;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5bp;gc_cont=0.687 MexT
ATGAACCGAAACGACCTGCGCCGCGTCGATCTGAACCTGCTGATCGTGTTCGAGACCCTGATGCACGAACGCAGCGTGACCCGCGCCGCAGAGAAACTGTTCCTCGGCCAGCCGGCCATCAGCGCCGCGCTGTCGCGCCTGCGCACGCTGTTCGACGACCCGCTGTTCGTCCGTACCGGACGCAGCATGGAGCCCACCGCGCGAGCCCAGGAAATCTTCGCCCACCTGTCGCCGGCGCTGGATTCCATCTCCACCGCCATGAGTCGCGCCAGCGAGTTCGATCCGGCGACCAGCACCGCGGTGTTCCGCATCGGCCTTTCCGACGACGTCGAGTTCGGCCTGTTGCCGCCCCTGCTCCGCCGCCTGCGCGCGGAGGCGCCGGGGATCGTCCTCGTCGTGCGCCGCGCCAACTATCTATTGATGCCGAACCTGCTGGCCTCGGGGGAGATCTCGGTGGGCGTCAGCTACACCGACGAACTGCCGGCCAACGCCAAGCGCAAGACCGTGCGCCGCAGCAAGCCGAAGATCCTCCGCGCCGACTCCGCGCCCGGCCAGCTGACCCTCGACGACTATTGCGCGCGACCGCACGCGCTGGTGTCCTTCGCCGGCGACCTCAGCGGCTTCGTCGACGAGGAGCTGGAAAAATTCGGCCGCAAGCGCAAGGTGGTCCTGGCGGTGCCGCAGTTCAACGGCCTCGGCACCCTCCTGGCCGGCACCGACATCATCGCCACCGTGCCCGACTACGCCGCCCAGGCGCTGATCGCCGCCGGCGGCCTACGCGCCGAGGACCCACCGTTCGAGACCCGCGCCTTCGAACTGTCGATGGCTTGGCGCGGCGCCCAGGACAACGATCCGGCCGAACGCTGGCTGCGCTCGCGGATCAGCATGTTCATCGGCGATCCGGACAGTCTCTGA
>CP077971.1_2661 # 2815108 # 2816127 # 1 # ID=1_2661;partial=00;start_type=ATG;rbs_motif=GAGG;rbs_spacer=6bp;gc_cont=0.669  MexS
ATGTCCCGAGTGATCCGTTTTCATCAGTTTGGCCCGCCAGAGGTCCTCAAATGCGAAGAGCTGCCGACCCCGGCGCCAGCCGCAGGGGAAGTCCTGGTGCGTGTCCAGGCGATCGGCGTGAGCTGGAAGGATGTGCTCTGGCGTCAGAACCTGGCCCCGGAGCAGGCTGCGCTGCCGTCCGGTCTCGGCTTCGAACTGGCCGGCGAGGTGCTGGCGGTCGGCGCCGGCGTCGGCGACCTGCCGCTGGGTTCCCGCGTGGCCAGTTTCCCCGCCCATACCCCCGATCATTATCCGGCCTATGGCGACGTGGTGCTGATGCCGCGCGCGGCCCTGGCGGTCTACCCCGAGGTACTCACCCCGGTGGAGGCCAGCGTCTACTACACCGGCCTGCTGGTGGCCTATTTCGGCCTGGTCGACCTGGCCGGGTTGAAGGCCGGGCAGACCGTGCTGATCACCGAGGCGGCGCGCATGTACGGGCCGGTCTCGATCCAGTTGGCCAAGGCTCTCGGCGCGCGGGTGATCGCTTCCACCAAGTCCGCCGAGGAGCGCGAGTTCCTCCGCGAGCAGGGCGCCGACAAGGTGGTGGTGACCGACGAGCAGGACCTGGTCCTGGAAGTCGAGCGCTTCACCGAGGGCAAGGGCGTCAATGTCATCCTCGACGAATTGGGCGGTCCGCAGATGACCCTGCTCGGCGATGTCTCCGCCACCCGCGGCAAGCTGGTGCTGTATGGCTGCAACGGCGGCAACGAGTCGGCGTTCCCGGCCTGCGCCGCGTTCAAGAAGCACCTGCAGTTCTACCGCCACTGCCTGATGGATTTCACCGGTCATCCGGAGATGGGCCTGGAACGCAACGACGAGTCGGTGAGCAAGGCCCTCGCGCACATCGAGCAACTGACCCGCGATCGCCTGCTCAAACCGGTGGTCGACCGGGTATTCGAGTTCGACCAGGTGGTCGAGGCGCACCGCTACATGGAAACCTGTCCGAAGCGCGGCCGGGTGGTGATCCACGTCGCCGATTGA

For example-

>CP077971.1_2660 |PA_PAO1|MexT  ATGAACCGAAACGACCTGCGCCGCGTCGATCTGAACCTGCTGATCGTGTTCGAGACCCTGATGCACGAACGCAGCGTGACCCGCGCCGCAGAGAAACTGTTCCTCGGCCAGCCGGCCATCAGCGCCGCGCTGTCGCGCCTGCGCACGCTGTTCGACGACCCGCTGTTCGTCCGTACCGGACGCAGCATGGAGCCCACCGCGCGAGCCCAGGAAATCTTCGCCCACCTGTCGCCGGCGCTGGATTCCATCTCCACCGCCATGAGTCGCGCCAGCGAGTTCGATCCGGCGACCAGCACCGCGGTGTTCCGCATCGGCCTTTCCGACGACGTCGAGTTCGGCCTGTTGCCGCCCCTGCTCCGCCGCCTGCGCGCGGAGGCGCCGGGGATCGTCCTCGTCGTGCGCCGCGCCAACTATCTATTGATGCCGAACCTGCTGGCCTCGGGGGAGATCTCGGTGGGCGTCAGCTACACCGACGAACTGCCGGCCAACGCCAAGCGCAAGACCGTGCGCCGCAGCAAGCCGAAGATCCTCCGCGCCGACTCCGCGCCCGGCCAGCTGACCCTCGACGACTATTGCGCGCGACCGCACGCGCTGGTGTCCTTCGCCGGCGACCTCAGCGGCTTCGTCGACGAGGAGCTGGAAAAATTCGGCCGCAAGCGCAAGGTGGTCCTGGCGGTGCCGCAGTTCAACGGCCTCGGCACCCTCCTGGCCGGCACCGACATCATCGCCACCGTGCCCGACTACGCCGCCCAGGCGCTGATCGCCGCCGGCGGCCTACGCGCCGAGGACCCACCGTTCGAGACCCGCGCCTTCGAACTGTCGATGGCTTGGCGCGGCGCCCAGGACAACGATCCGGCCGAACGCTGGCTGCGCTCGCGGATCAGCATGTTCATCGGCGATCCGGACAGTCTCTGA

Thank you!

fasta
• 1,215 views • link
written by Neel • Neel

1 answer

The first step is pretty straightforward: sed "/^>/s/#.*$//g" will delete everything after the first # till the end of the line. This will leave you with such a file

>CP077971.1_2660
MexT
ATGAACCGAAACGACCTGCGCCGCGTCGATCTGAACCTGCTGATCGTGTTCGAGACCCTGATGCACGAACGCAGCGTGACCCGCGCCGCAGAGAAACTGTTCCTCGGCCAGCCGGCCATCAGCGCCGCGCTGTCGCGCCTGCGCACGCTGTTCGACGACCCGCTGTTCGTCCGTACCGGACGCAGCATGGAGCCCACCGCGCGAGCCCAGGAAATCTTCGCCCACCTGTCGCCGGCGCTGGATTCCATCTCCACCGCCATGAGTCGCGCCAGCGAGTTCGATCCGGCGACCAGCACCGCGGTGTTCCGCATCGGCCTTTCCGACGACGTCGAGTTCGGCCTGTTGCCGCCCCTGCTCCGCCGCCTGCGCGCGGAGGCGCCGGGGATCGTCCTCGTCGTGCGCCGCGCCAACTATCTATTGATGCCGAACCTGCTGGCCTCGGGGGAGATCTCGGTGGGCGTCAGCTACACCGACGAACTGCCGGCCAACGCCAAGCGCAAGACCGTGCGCCGCAGCAAGCCGAAGATCCTCCGCGCCGACTCCGCGCCCGGCCAGCTGACCCTCGACGACTATTGCGCGCGACCGCACGCGCTGGTGTCCTTCGCCGGCGACCTCAGCGGCTTCGTCGACGAGGAGCTGGAAAAATTCGGCCGCAAGCGCAAGGTGGTCCTGGCGGTGCCGCAGTTCAACGGCCTCGGCACCCTCCTGGCCGGCACCGACATCATCGCCACCGTGCCCGACTACGCCGCCCAGGCGCTGATCGCCGCCGGCGGCCTACGCGCCGAGGACCCACCGTTCGAGACCCGCGCCTTCGAACTGTCGATGGCTTGGCGCGGCGCCCAGGACAACGATCCGGCCGAACGCTGGCTGCGCTCGCGGATCAGCATGTTCATCGGCGATCCGGACAGTCTCTGA

This file, you could reshape accordingly with paste -d \| - - - , however I suspect that the example you provided is somehow mangled, because that MexT in a separate line is no valid Fasta. Plus, where does that PA_PAO1 in the desired output come from?

• 0 views • link
written by Matthias Zepper Matthias Zepper

Thank you so much for your time, Actually i want to add strain name also so that i can track later which gene from which strain.

• 0 views • link
written by Neel • Neel

Log in to answer this question.

More posts like this