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.

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • pyspark.sql模塊 模塊上下文 Spark SQL和DataFrames的重要類: pyspark.sql...
    mpro閱讀 9,892評論 0 13
  • Introduction What is Bowtie 2? Bowtie 2 is an ultrafast a...
    wzz閱讀 6,158評論 0 5
  • The Great A.I. Awakening How Google used artificial intel...
    圖羽閱讀 1,456評論 0 3
  • (一)娃他舅 建筑隊里有一個被他們喚作“娃他舅”的人,三十多歲,蓄著特朗普一樣風(fēng)吹一邊倒的長發(fā),胡子邋里邋遢的,遠...
    明望舒閱讀 475評論 0 0
  • 今天按摩的力度加強了速度也加快了。前幾天按摩時發(fā)現(xiàn)的痛點已經(jīng)記住了不用摸索,并每天都有一兩點新的痛點發(fā)現(xiàn),就像...
    請叫我筠竹閱讀 272評論 0 0

友情鏈接更多精彩內(nèi)容