1. kmer评估中杂合度为多少?高杂合是否组装为两套染色体?高杂合未区分,D高2.物种是否为多倍体?有全基因组复制事件?如果存在,D高是正常情况。3.参考近缘种,与其他人busco的结果进行比较,看差异大不大,如果近缘种都是这样的结果可能是物种本身特性。后续处理,要么用hifiasm做出来的单倍型进行后续分析,或者采用 Pu...
回答于 2025-09-17 09:39
运行代码是:Rscript $scriptdir/trait_plot.r -i $datadir/traits_gapit.txt ?格式是这个:
回答于 2025-09-15 16:50
检查输入文件./NorthAfrica-Yourawildgoat.sfs是不是有一行全是0。easySFS.py里面投影值设置根据自己样本数目进行更改
回答于 2025-09-15 16:35
问题出现在coding_change.pl脚本的,77行存在一个正则表达的匹配要求: $field[2] =~ m/^[\w\-\.\@\/]+?:([\w\.]+?):exon\d+:c.([\w\->]+)(:p.([\w\*]+))?/ or die "Error: invalid record found in exonic_variant_function fi在这一行上要求第三列的的第二部分(也就是rna名字)是([\w\.]+?)这个范围内的东西,而你的...
回答于 2025-09-11 18:01
虚拟机里面是linux系统的话使用起来没差别。压缩包上传到指定位置,然后docker load -i 加载压缩包就行了
回答于 2025-09-10 13:18