文件contig.fa.divsum是calcDivergenceFromAlign.pl生成的,不是repeatmasker的结果
我运行下面的你们给的代码
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正常输出的文件都有,请问是不是你们的过程中遗失了步骤,或者给的代码不正确?!?!