Ti Amo
Ti Amo

性别: 注册于 2023-07-19

事已至此,先吃饭吧

向TA求助
840金币数
6200 经验值
2个粉丝
主页被访问 7711 次

97 个回答

0 赞同

mummer比对syri进行SV检测报错

chrorder调整的时候染色体条数发生改变了。不做这一步直接根据之前的delta和coords文件做syri。但是可能存在大结构变异也会导致call变异失败

回答于 1天前

0 赞同

基因组组装BUSCO多拷贝过高

1. kmer评估中杂合度为多少?高杂合是否组装为两套染色体?高杂合未区分,D高2.物种是否为多倍体?有全基因组复制事件?如果存在,D高是正常情况。3.参考近缘种,与其他人busco的结果进行比较,看差异大不大,如果近缘种都是这样的结果可能是物种本身特性。后续处理,要么用hifiasm做出来的单倍型进行后续分析,或者采用 Pu...

回答于 1天前

0 赞同

xpclr

样本txt里面的名字和vcf的是否一致你的染色体名字之间没有空格,会认为是一个很长名字的染色体。

回答于 2天前

0 赞同

请问我在做性状绘图时出现报错,我的数据格式应如何修改

运行代码是:Rscript $scriptdir/trait_plot.r  -i $datadir/traits_gapit.txt ?格式是这个:

回答于 3天前

0 赞同

dadi

检查输入文件./NorthAfrica-Yourawildgoat.sfs是不是有一行全是0。easySFS.py里面投影值设置根据自己样本数目进行更改

回答于 3天前

0 赞同

请问我在做GWAS性状绘图时出现以下报错,我的数据格式该怎样修改...

很可能表型数据格式不符合标准。可以给出表型数据截图、写明自己操作命令以及报错。

回答于 3天前

0 赞同

老师想问一下,为什么我在myt2t文件下面创建,但是最后跑去了家...

变量workdir没赋值。export workdir=/work/myt2t

回答于 3天前

0 赞同

dadi分析

问题出现在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

0 赞同

基因组过大,能不能分割基因组做重复序列注释

可以这样加速,切割部分的注释结果可能需要重新做一下

回答于 2025-09-11 10:01

0 赞同

docker 安装

参考【docker安装】由于没有公钥,无法验证下列签名: NO_PUBKEY 7EA0A9C3F273FCD8-CSDN博客

回答于 2025-09-10 16:06