3 个回答
0 赞同
MCSCANX共线性图代码报错
老师，我是用您们的biolinux, 也cd到MCScanX的安装目录执行，检查了GFF文件，现在运行GFF文件没有报错，但family_circle_plotter报了其他行的错误,253行， loc2=gene_feat.get(gene2.get(i)).toString(); 另外MADS.circle.PNG是自动生成，还是要先手动创建这个文件？是不是文件夹或文件读取权限有需要设置，我是win10 64...
回答于 2018-10-15 15:15
1 赞同
提取蛋白编号信息
不好意思，第一次提问，下次一定写的详细来问您 我是这样运行的： 先创建了pep.fa文件 然后执行 [share] perl get_fa_by_id.pl AT.gff TAIR10.pep.all.fa pep.fa AT.gff是之前提取的基因位置信息，没有问题 运行完没有报错，但 pep.fa文件是0字节
回答于 2018-10-15 13:35
1 赞同
替换蛋白编号信息的问题
@omicsgene 老师，按照视频的脚本仔细核对没有错误，编译执行后也有没有报错，但pep.fa文件一直是空，文件名字没有空格， 如果脚本有问题，编译起来就会报错，还得请教老师来解决，谢谢 use Bio::SeqIO; use Bio::Seq; my$in = Bio::SeqIO->new(-file => "$ARGV[1]" , -format => 'Fasta'); my$out = B...
回答于 2018-10-15 00:08