I am using the prody package to do some co-evolution analysis on my multiple sequence alignment (MSA).
When I parse the MSA file using prody I get a Bus error: 10
Code:
from prody import *
msa = parseMSA('test_aln.fasta')
This outputs the following in my terminal:
>>> from prody import *
>>> msa = parseMSA('cbpa_aln.fasta')
Bus error: 10
I can't figure out why i am getting this error. I've tried restarting my machine and re-installing prody but to no avail.
Test Data:
>CbpA
--------------------------------------------------
------------ASALKDYYAIMGVKPTDDLKTIKTAYRRLARKYHPDVS
KEP---DAEARFKEVAEAWEVLSDEQ--RRAEYDQMWQH---RN-DP--Q
FNRQFHHGDG----------------------------------------
---------QSFNAEDFDDIFSSIFGQHA----RQ---------------
---------------SRQRPATRGHDIEIEVAVFLEETLTEHKRTISYNL
PVYNA-FG---------------MIEQEIPKTLNVKIPAGVGNGQRIRLK
GQGTPGENG-----GPNGDLWLVIHIAPHPL-FDIVGQDLEIVVPVSPWE
AALGAKVTVPTLK-ESILLTIPPGSQAGQRLRVKGKGLVSK---KQT--G
DLYAVLKIVMPPKP-DENTAALWQQLADAQSSFDPRKDWGKA--------
---
>gi|422828735|ref|ZP_16876905.1|
--------------------------------------------------
-------------MELKDYYAIMGVKPTDDLKTIKTAYRRLARKYHPDVS
KEP---DAEARFKEVAEAWEVLSDEQ--RRAEYDQMWQH---RN-DP--Q
FNRQFHHGDG----------------------------------------
---------QSFNAEDFDDIFSSIFGQHA----RQ---------------
---------------SRQRPATRGHDIEIEVAVFLEETLTEHKRTISYNL
PVYNA-FG---------------MIEQEIPKTLNVKIPAGVGNGQRIRLK
GQGTPGENG-----GPNGDLWLVIHIAPHPL-FDIVGQDLEIVVPVSPWE
AALGAKVTVPTLK-ESILLTIPPGSQAGQRLRVKGKGLVSK---KQT--G
DLYAVLKIVMPPKP-DENTAALWQQLADAQSSFEPRKDWGKA--------
---