This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Can we use JELLYFISH for amino acid seqeunces?

I have a very large file of 1000's of amino acid sequences. I would like to identify ALL overlapping k-mers.

I was hoping to use JELLYFISH, but it only seems to work with nucleic acid sequences. Is there any way to use the tool to read and parse the fasta file of amino acids?

jellyfish genome sequencing
• 2,577 views • link
updated by Biostar Biostar • written by marksingh1982 • marksingh1982

This is a question; you're not posting about a tool. I've made the required change now, please be more careful in the future.

• 0 views • link
written by Ram Ram

1 answer

Jellyfish cannot deal with the protein sequences. You can use skbio iter_kmers() for this. BioPython solution:

from Bio import SeqIO
from skbio import Sequence

myfile=SeqIO.parse('test.fa','fasta')
for record in myfile:
    sequence=Sequence(str(record.seq))
    for kmer in sequence.iter_kmers(4, overlap=True):
        print(str(kmer))
• 0 views • link
written by Sej Modha Sej Modha

Log in to answer this question.

More posts like this