RepeatMasker输出没有.divsum文件

我运行下面的你们给的代码

RepeatMasker -e rmblast -gff -xsmall -html -norna   -source  -no_is -pa ${threads} -s -lib allRepeats.final.lib contig.fa
# -e rmblast 指定搜索引擎为rmblast,还可以选择crossmatch、abblast或者hmmer
# -species "Danio rerio" 指定物种为斑马鱼,物种名必须在NCBI物种分类数据库(Taxonomy)中,推荐使用拉丁学名
# -pa 12 使用12线程运行RepeatMasker
# -s 慢速搜索,灵敏度比高0~5%,速度慢2-3倍
# -xsmall 软屏蔽,将重复序列变成小写字母,而不是N或者X。
# -gff 输出GFF文件
# -nolow 忽略低重复的序列
# -norna 忽略Small RNA

下一步运行代码是

# 生成RepeatLandscape  ; -g 设置基因组大小计算屏蔽比例

contig_size=$(esl-seqstat --dna contig.fa| grep "Total " | tr -s [:space:] | sed 's/: /\t/g' | cut -f 2)
calcDivergenceFromAlign.pl -s contig.fa.divsum contig.fa.cat.gz
createRepeatLandscape.pl -div contig.fa.divsum -g $contig_size > contig.fa.html

我的代码全面运行都非常的正常,但是RepeatMasker输出的没有contig.fa.divsum文件!其它RepeatMasker正常输出的文件都有,请问是不是你们的过程中遗失了步骤,或者给的代码不正确?!?!

请先 登录 后评论

1 个回答

Ti Amo

文件contig.fa.divsumcalcDivergenceFromAlign.pl生成的,不是repeatmasker的结果

请先 登录 后评论
  • 2 关注
  • 0 收藏,108 浏览
  • litianhuan 提出于 2024-09-11 09:57

相似问题