找到约 15 条结果

问题 braker.pl运行中断

...夹下的gmes.log有显示error,但是没提醒具体的问题。请问如何解决?是用的docker运行genomic-annotation这个镜像

文章 NCBI上下载参考基因组信息

...个界面,有不同组装完整的筛选,系统会自动为我们筛选比较好的基因组,一般来说我们选择level这一列中,测序水平最高的一组。就是黑色部分占据面积最多的哈~(比如下面),size最大的一个就可以,点击Assembly列链接: ...

文章 linux目录结构详细讲解

...状态,由于数据都存放于内存中,所以不占用磁盘空间,比较重要的目录有/proc/cpuinfo、/proc/interrupts、/proc/dma、/proc/ioports、/proc/net/*等。/root:系统管理员root的宿主目录,系统第一个启动的分区为/,所以最好将/root和/放置在一个...

问题 请问R/QTL这个包可以用来做果树群体的qtl定位吗?

请问R/QTL这个包可以用来做果树群体的qtl定位吗?如果可以的话,基因型该如何分型

问题 利用人类基因组进行群体遗传PSMC分析求助

...照染色体分的,一共22个常染色体22个bam文件,那么他们如何作为输入文件写进代码里面呢?我看代码输入文件是一个bam文件,我是需要把bam文件进行 samtools merge吗?

问题 进行GTDB注释时报错,报错信息为缺少参考基因组

老师您好,我在进行GTDB注释时,根据课程所给的代码下载了GTDB的数据库信息,但进行注释时,却出现了缺少参考基因组的信息,具体情况见下图。请问该如何解决?

问题 怎么将一个perl命令(或者说脚本)变成shell脚本的形式

如何将这行命令 perl bin/indel_p3in.pl QT4indels2.in genome.fa con.txt,写成一个shell脚本即*.sh的形式

文章 芯片灰度图

...低强度的信号点,刮擦的痕迹,背景太高等等,整个图像比较均一,背景均一。 基于affy包中的image()函数,对AffyBatch对象直接可以绘制出每一个样本芯片的灰度图, 以GSE66196数据为例,绘制第一个样本灰度图。 image(AffyBatch[,1])...

文章 JGI-IMG微生物功能基因组数据库介绍

...itute,JGI)于2005年创立,是综合的微生物基因组数据库及比较分析系统。IMG收录了细菌、古菌、质粒、病毒以及少量真核生物基因组数据,其数据主要来源于NCBI的RefSeq数据库,但是增添了更加详细的注释信息,例如CRISPR序列、...

问题 问一个玉米基因芯片ID转换的问题

...发现芯片文件(GPL.soft)中这些ID都很生僻(如图),不知如何转换成传统ID? 貌似那个TIGR ID已经绝版了?我在TIGR ID的对应网站无法下载到对应的转换文件(http://maize.jcvi.org/)。 请问怎么办? 谢谢! 注:这是该数据库涉及...

文章 检查kegg的map是否能和差异基因对上

...cut -d "(" -f1 > ko-gene 4. 和我们的基因-KO编号list文件进行比较 cut -d "|" -f1  KOlist | sed  "s/\t/ /g" > kegg.kocat ko |while read -r line; do  if grep -q " [[:space:]]*$line$" kegg.ko ; then echo "$line";else echo "false" ;fi ;done > result result文件里面如果...

问题 LDdecay

...0,用的是自己的数据,前一步都还正常,不知道popid.txt如何设置,这里是按照示例文件,放了所有样品的集合,但是跑出来是如下这样 the Number of subPop samples[found in VCF] is 0 sub Group Population szie is too small, SubGroup sample size: 0 ##begi...

问题 基因家族鉴定所用的基因ID与转录组基因ID不一致

老师您好,我做有参转录组的基因ID与现在做基因家族鉴定的基因ID不一致,后续做基因家族成员转录表达模式分析的时候如何将基因ID一致化?

问题 用pfam搜到的基因构建进化树报错

...建树,基因家族的鉴定是不是就不完整了?不删除的话,如何建树?

文章 perl 从数组随机抽取元素

文章简单介绍一下perl中如何随机抽取数组中的一个元素。 这里用到List::Util模块中的shuffle子程序,List::Util模块是Perl中专门用于处理列表数据的一个模块。 方法如下: use List::Util qw(shuffle);my @shuffled_nums = shuffle(1..10);   # ...