找到约 15 条结果

问题 进化树构建

老师您好,利snp构建进化树时,在vcf文件格式转换为supergene.phy这一步,生成的supergene.phy文件很多位点是0(左图),导致后续利iqtree2构建进化树报错(右图),想请教老师问题出在哪?

问题 braker注释结果有问题

braker.gff3文件第三列未出现gene、mRNA等层级这个问题该怎么解决

问题 docker GetOrganelle的批量组装

如何修改课程使的docker代码,以实现GetOrganelle的批量组装

问题 微生物16s分析

OTU分析按照视频方法一(pick_de_novo_otus.py)做的,显示usearch版本不对

文章 pfam 隐马尔可夫模型序列最新下载方法

HMMER 隐马尔可夫模型序列最新下载方法: 方法一: 利Pfam网站下载,此网站可以到2023年1月,之后将停,现在如果使PFam网站下载,方法如下:(此方法仅适合在2023年1月之前使): 打开 ② 以上下载的是结果...

文章 如何预测同一基因不同转录本对应蛋白的功能域差异

在真核生物,由于转录后可变剪接的存在,赋予了基因更多的功能,一个基因会产生多个mRNA,不同mRNA所翻译出来的蛋白功能是有很大差异的,那么了解不同转录本翻译的蛋白各是什么功能是一个非常必要的工作,今天我就给...

文章 Linux命令处理文件合并

...编程语言处理会非常慢,而且容易出错。 可以使Linux 的paste和awk处理: 命令如下: paste file1.txt file2.txt | awk 'BEGIN {OFS="\t"} {print $1, $2, $3, $4 + $11, $5 + $12, $6, $7}'  > result.txt 使paste命令将两个文件合并,再使awk命令计...

问题 基因串联重复

老师,您好! 请问老师,基因组普遍存在的基因串联重复导致的基因多拷贝,形成的原因是什么呢?它们的形成机制是什么呢? 比如有些基因重复,是由于转座子的介导,全基因组复制导致的,那基因串联重复又是由于什...

文章 eclipse下搭建shell脚本编辑器--安装开发shell的eclipse插件shelled

... 通过添加的方式设置shell解释器,直接选定解压后目录的sh.exe即可。 然后写一个shell测试一下是否能运行: shell文件最终是要在linux运行的,在windows的eclipse运行shell最多只能来调试,解除一些基本的BUG,大家在这里开发...

文章 2024最新版《基因家族分析实操》直播培训开始报名了!

...加组学大讲堂直播课程,两天时间实现SCI自由! 我们使关键词“genome wide gene family”粗略统计了近3年基因家族分析类SCI数量:2021年1912篇,2022年2111篇,2023年1890篇。 考虑到从2023年开始不少基因家族文章的题目缺少了传统...

问题 TCGA生存分析

cutoff为什么是2500?是从这个图里得出来的吗?视频里提到2500是位生存时间,但是生存曲线是在ROC之后做的; 怎样处理表达矩阵样品重复的问题?代码只对clinic进行了unique

文章 天宫一号曾为农业做出的贡献!

...衰减期。北京时间2018年4月2日早 8 时 15 分左右,在太空飞行了六年半的天宫一号目标飞行器已再入大气层,再入落区位于南太平洋部区域,绝大部分器件在再入大气层过程烧蚀销毁,结束它的历史使命。 作为一名学...

问题 WGCNA相关问题

老师您好,我按照课程给的代码读取WGCNA基因表达数据,数据格式与示例数据相同,为什么每次dim之后就只有0行和1列呢,后面的筛选也筛选不到基因,聚类树也画不出来,求解答

问题 基因家族分析-基因家族成员鉴定

老师, 在NCBI-CDD, 如果显示结构域不完整, 是否将这些基因去除?

问题 GATK 分染色体call 变异时出现错误

如下命令:  for i in $(cat $workdir/data/data.txt) ; do  for Chr in Chr1 Chr2 Chr3 Chr4 Chr5 Chr6 Chr7 Chr8 Chr9 Chr10 Chr11 Chr12; do   echo "gatk --java-options "-Xmx100g" HaplotypeCaller \     -R $REF \     -I $workdir/3.map/result/${i}.sorted.dedup.bam \   ...