找到约 15 条结果

问题 PASA预测基因数目过少

老师好,我在运行PASA_assembly时,因是真菌,best_candidate.gff3有6514个基因,有编码潜力只有588个,最后筛选high_quality_gene.gff3只有1538,过滤掉距离太短只剩下994。您说2500个基因左右即可,但是我预测出来基因远低于...

问题 hmmer搜索结构域

想修改WRKY_hmmerOut_final.txt中结构域18、19列也就是结构域起始、终止位置,不知道具体在哪一列,还有上一次默认做出来图中发现后面普遍锌指结构不全,应该是要整体后移以下,如图:,结构域缺失较多基因,CDD数据库中...

文章 火山云 | tos:// 数据下载 | windows/linux下载方式

...网址:https://www.volcengine.com/docs/6349/148777 01-2:双击下载安装包.exe,即可安装打开。 01-3:在登录页面选择“共享连接登录”,输入 共享连接 ,即可登录下载。 01-4:使用共享链接登录 TOS Browser 后,只拥有对应文件或文件...

文章 绘制venn图报错Unexpected parameter length for "category"

在进行韦恩图绘制时候出现报错: 出现这个报错原因是在输入数据应当是一个有名字list,要求是这个list中每个元素都有名字。 venn <- list(A,B,C)names(venn) <- c("A","B","C") 这样之后再进行绘制就没有问题了 附:...

问题 dadi

这是我用命令 cd $dir/sfs awk '{ print $2}' $data/africa.group.txt | sort | uniq -c  # 无外群,生成foldfs easySFS.py -i$data/Africadadi.vcf.gz -p $data/africa.group.txt -a --proj=74,56,60 -o./ ## 生成了3个群体、以及三个群体两两之间sfs文件 cd dadi/ mkdi...

问题 老师我在构建泛基因列表,合成last.txt文件,结果只有两列

老师,按照您给指令(for i in `find ../01.jcvi -name "pangene_filter.*.last"`;do     name=`basename ${i} | cut -d "." -f2`;     grep -v "^#" ${i} > ${name}.last;     python3 $script/find_1vs1_v2.py ${name}.last ${name}.txt && /bin/rm ${name}.last; Done cat *....

问题 根据课程基因组重测序脚本,进行片段inner size,片段选择是否异常和测序深度与覆盖度统计均未能生成pdf文件,直接弹出如下对话框和如下命令

 Exception in thread "main" java.awt.AWTError: Can't connect to X11 window server using 'localhost:12.0' as the value of the DISPLAY variable.

问题 麻烦求助,在使用ncbi_gff_2_Ensembl_gff.pl这个脚本时出现报错显示:cant find the string terminator "Usage end"...in line 10

问题 老师从gff文件中用mRNAID_to_geneid.pl提取gene和mRNA之间对应关系时报错Protocol error是什么原因,脚本没有自己更改,用别人电脑可以,在自己电脑操作就不行。

问题 请问,首次使用MCscan (Python version)可视化作图时,没有mcscan_seqids和mcscan_layout文件,该怎么创建这两个配置文件呢?自己创建配置文件运行后,结果文件有问题

问题 请问老师有没有get_gene_position.pl,get_gene_bed.pl,get_fa_by_id_from_bed.pl这个三个脚本下载地址?

老师您好! 我是课程购买者,但是这三个脚本实在没找到下载地址,请问哪有?谢谢!

问题 eggNOG 保留蛋白序列时提示同一 ID 被多个一级特征使用,而子特征 Parent 属性没有得到相应更新,导致注释关系不明确问题,该怎么解决

...arning提示,gff和fasta都是非模式物种ncbi数据库染色体水平组装,ga_protein_coding.gff3文件正常生成了,请问warning对结果有影响吗,应该如何解决?谢谢! agat_sp_filter_feature_by_attribute_value.pl --gff $gff --attribute gene_biotype --value protein_...

问题 老师,我在利用perl ../script/get _ene _ exon_ from _ gff.pl in1 时,老是提不出文件,他说十几行不合适但是这个和视频里面一样,

试了一天了,就是不合适

问题 老师,转录组课程中一个基因提取一个蛋白序列作基因代表序列用于同源注释,我152个基因,但是只提取出138个蛋白序列,没提取出来在总蛋白库里都有,想问下啥提不出全部

问题 您好!我是R语言新手,在用hchinamap程序绘制地图时,出现cannot coerce type 'closure' to vector of type 'character',出现这类错误怎样解决呢?请指教,谢谢!

步骤如下: install.packages('hchinamap') download.file('https://mdniceczx.oss-cn-beijing.aliyuncs.com/chinadf.rda', file.path(dir, 'chinadf.rda')) load(file.path(dir, 'chinadf.rda'), verbose = TRUE)chinadf