I hope someone could help me with a simple line command.
I have a protein FASTA file:
>Protein1 sequence >protein2 sequence
I want to extract the sequences that has specifics peptides (ERLE). I could do that using grep -A, however it prints the sequence and the header of the below protein. For example, let say the sequence is on protein 1, when I use grep it prints the sequence 1 whit the header of protein 2:
Is there a way to solve that?
Thank you in advance