各位老师同学好,我用mitoz组装线粒体基因组,使用的默认组装器megahit,结果发现其中一条序列的蛋白质编码基因中间出现了终止密码子,查询官网后,得知原因有以下4点:(1)Ns in the PCGs;(2)assembly error;(3)wrong genetic code table for translation ;(4)nuclear mitochondrial DNA segments (NUMTs),其中密码子表的选用肯定是没问题的,但是我不知道如何去排查其他原因,我就使用mitoz另外两个组装器去组装,spades出现报错,mitoassemble组装了32个小时依然没结束。我想问问大家遇到这些情况是如何去解决的?谢谢~
以下是我使用spades和mitoassemble的代码以及图片:
spades:mitoz all --genetic_code 5 --clade Arthropoda --outprefix R1124 --thread_number 12 --fq1 R1124_L3_UDI020.R1.fastq.gz --fq2 R1124_L3_UDI020.R2.fastq.gz --fastq_read_length 150 --insert_size 300 --memory 16 --requiring_taxa 'Arthropoda' --assembler spades --skip_filter
mitoassemble:mitoz all --genetic code 5 --clade Arthropoda --outprefix R1124 --thread number 12 --fq1 R1124 L3 UDI020.R1.fastq.gz --fq2 R1124_L3_UDI020.R2.fastq.gz --fastq_read_length 150 --insert_size 300 --memory 16 --requiring_taxa 'Arthropoda' --assembler --skip_filter
以下是我使用spades和mitoassemble的代码以及图片:
spades:mitoz all --genetic_code 5 --clade Arthropoda --outprefix R1124 --thread_number 12 --fq1 R1124_L3_UDI020.R1.fastq.gz --fq2 R1124_L3_UDI020.R2.fastq.gz --fastq_read_length 150 --insert_size 300 --memory 16 --requiring_taxa 'Arthropoda' --assembler spades --skip_filter
mitoassemble:mitoz all --genetic code 5 --clade Arthropoda --outprefix R1124 --thread number 12 --fq1 R1124 L3 UDI020.R1.fastq.gz --fq2 R1124_L3_UDI020.R2.fastq.gz --fastq_read_length 150 --insert_size 300 --memory 16 --requiring_taxa 'Arthropoda' --assembler --skip_filter