pssm.search from Bio.motifs not working with amino acid sequences
I am trying to find a motif out of a list of amino acid sequences.I am using the package Bio.motifs. I followed the documentation as in theBiopython Tutorial and Cookbook chapter 14.
Everything go smoothly and I generate the PWM and the PSSM out of my sequences. Then I wanted to use this PSSM to look for that motif in a new sequence. But when I try to use the method search (as in 14.6.2) with my PSSM I get this error:
ValueError: PSSM has wrong alphabet: IUPACProtein() - Use only with DNA motifs
So I have two question:
- Is the package
Bio.motifssuitable for proteins (from documentation is not clear anymore to me)? And can I trust the PSSM I generated? - How can I now use this PSSM to look for the motif? Is there another program or python package that allow me to do so? I searched online and I could not find anything so far.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Mapping PWM of transcription factor motifs to the human genome
written by John •
-
Bio.motifs.threshold module
written by nyanovsky •
-
Search for specific motif in MEME analysis
written by Laura •
-
Error when searching for matches using the PSSM score with proteins (.motif BioPython)
written by Agenor Neto •
-
Biopython Bio.motifs: How to create a motif object with aligned sequences
written by kinetic •
-
Convert binding site formats from different databases and compare the motifs - any tools available?
written by JJ
-
Predicted target genes - which tools are the best to use?
written by JJ
-
How To Use The Biopython Motifs Module To Check Whether A Motif Is Present In My List Of Sequences?
written by ericmajinglong
-
Wrong Alphabet In Biopython.Motifs When Trying To Pssm.Calculate
written by Jelena_bioinf •
-
Import Fasta Sequences To A Motif
written by RossCampbell
