2019-03-14

Hey, Please download raw data (https://github.com/linnabrown/error_data) to test the code.
The error of running macsypy:

Sequence number 17..Traceback (most recent call last):
  File "/home/ciillab/huangle/CRISPRCasFinder/bin/macsyfinder", line 329, in <module>
    idx.build(force = config.build_indexes)
  File "/home/ciillab/huangle/CRISPRCasFinder/macsyfinder-1.0.5/macsypy/database.py", line 118, in build
    raise RuntimeError(msg)
RuntimeError: an error occurred during databases indexation see formatdb.log

What's more, I saw the code snapets
is related to it:

#################################
# synchronization point between #
# hmmer_indexes and my_indexes  #
#################################
        if force or not hmmer_indexes:
            hmmer_indexes_proc.wait()
            if hmmer_indexes_proc.returncode == 127:
                msg = "neither makeblastdb nor formatdb can be found, check your config or install makeblastb"
                _log.critical(msg, exc_info=True)
                raise RuntimeError(msg)
            if hmmer_indexes_proc.returncode != 0:
                msg = "an error occurred during databases indexation see formatdb.log"
                _log.critical(msg, exc_info=True)
                raise RuntimeError(msg)
        self._hmmer_indexes = self.find_hmmer_indexes()
        self._my_indexes = self.find_my_indexes()
        assert self._hmmer_indexes, "failed to create hmmer indexes"
        assert self._my_indexes, "failed create macsyfinder indexes"

This is the wrong code of this sequence. The length of genome sequence may be too small compared with other genome sequences in the file of GCF_001983215.1_ASM198321v1_genomic.fna

 36723 >NZ_MOWP01000024.1 Frankia sp. CcI49 contig_23, whole genome shotgun sequence
 36724 GACCCGGTCGACCCACCCTGGGACCCGCCCTTACCGAACGACCCGGTCGACCCACCCTGAGCGCCACCCTTACCGAACGA
 36725 CCCGGTCGTCTCCTTGCCGCCAAGGTTCCGCGGTGTGGAGATATCGATGGTCACGCCTTCCTCAACCGGTCTGGTTGTTT
 36726 GGGATGACCCGCCCTGGGTTCCGTCGCCATCGAACGAGCCAGTCGACCCACCCTGAGCGCTGGACGCCGCACCGAAGGTA
 36727 CTCGTCGGGCCACCGACAATGGTCTCCGCAGTAACCGGCGCGCCAATACCGATCGTGATGCTTCCGGTCGACAGAATGGA
 36728 GAGGTTCAATCCGCTCGGAATAGCCAGGCCCTGCCCGGCATCCAAAACCCTGGGCTCCACTTCTCCGACACTCTCCTGGC
 36729 CATCTCCTACAGCGACCACGGGCGCCACCAGGTTGGCTGGAGCCGGATCAACGGCTATCGATATGTCAGCCCTGGCAGGA
 36730 GACGCCGGCAAGATAGCAGCAGCAGATATTCCACCGGCCAGAGCCAGAACGGACCAAGTCTTCTTAACAATCATGATTTC
 36731 CCCGCTCGACATCAGGATCCACCTCGTCGCAGCGGGATGGTAGCTATCAAGCCGTCGGATCAACCGCTTCCACCAGACCG
 36732 GCGTGTTGCGGAAAAATTTCGAGTGACATATTTCGCGACAGCAAAGCGTCGAACTCATCCGGACCGGATCACATGGCAAA
 36733 AGTAGATCAGATCACCAACCCGGCCCACCTGGTGTGATTGGTGCGAGTACACACCCCTATCAGGGAGATCCTCGCAACAT
 36734 CATTCCACCCAGCTCTTCATGCATGCGACAGCCCTCCGCCTCTCCCACCTGGCAGATTTCTCCAGATGGGAGAAGGGAGG
 36735 GCTGCCCGTCACACCCAAACAAAAATTAGGCGGAC

Please download all the data (faa, fna, gff) from my github to test it.

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • pyspark.sql模块 模块上下文 Spark SQL和DataFrames的重要类: pyspark.sql...
    mpro阅读 9,532评论 0 13
  • Introduction What is Bowtie 2? Bowtie 2 is an ultrafast a...
    wzz阅读 5,823评论 0 5
  • The Great A.I. Awakening How Google used artificial intel...
    图羽阅读 1,253评论 0 3
  • (一)娃他舅 建筑队里有一个被他们唤作“娃他舅”的人,三十多岁,蓄着特朗普一样风吹一边倒的长发,胡子邋里邋遢的,远...
    明望舒阅读 433评论 0 0
  • 今天按摩的力度加强了速度也加快了。前几天按摩时发现的痛点已经记住了不用摸索,并每天都有一两点新的痛点发现,就像...
    请叫我筠竹阅读 263评论 0 0