Extract specific peptides from a FASTA file
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:
sequence
>Protein2
Is there a way to solve that?
Thank you in advance
• 1,840 views
•
link
written
by
ijvechetti •
1 answer
Log in to answer this question.
More posts like this
-
Making pairwise matrix using ChIP-Seq peak binding matrix in R
written by Ankit
-
How to get Protein names and their ID for a given list of peptide sequence (using Python)
written by Mathew •
-
Replace names in FASTA headers from the first white space
written by bionix •
-
How to make blast go on to the next protein
written by marcelofolhas •
-
Extracting the coordinates for a specific base from multifasta
written by Louis Kok •
-
Calculate the coverage of a protein having a list of its peptides
written by arronar
-
how to write a script to grep information from another file without losing location information(in …
written by Anny •
-
calculate the length of a sequence after adding the length of previous sequences
written by User 6777 •
-
Any tool to generate protein similarity matrix
written by Zealseeker •
-
How To Find Paralogs Using Blast
written by ngsgene
Thank you so much. Worked really well