根据下面的代码,我用mummer比较两个基因组得到out.delta,用syri处理它们得到filter和coords文件,然后用syri检测sv报错:Reading coords - ERROR - aDir can only have values 1。
nucmer --mum -c 100 -b 500 -l 50 /home/zzy/Genome/Red_crucian_carp.genome.fasta /home/zzy/Genome/common_carp.genome.fna
delta-filter -1 -m -i 90 out.delta > out.filtered.delta
show-coords -THrcld out.filtered.delta > out.filtered.coords
syri -c out.filtered.coords -d out.filtered.delta -r ref.fasta -q qur.fasta
我想知道这是不是因为参考基因组是染色体水平,而另一个不是。