脚本认为ID=与;之间的为gene或mRNA的ID,建议你把GFF文件里面的gene: CDS: transcript: 删除; 然后再看看fasta里面的ID与GFF里面的ID是否一致。 这样保持所有地方的ID一致,就不会因为找不到对应的ID而出错。 如果上面按照我说的都改了,脚本这里27行就可以删除gene: 获取ID信息了:
回答于 2019-11-08 10:22
这是ggplot2里面内置的色板,如果颜色不够可以使用手动设置色板scale_colour_manual():https://ggplot2.tidyverse.org/reference/scale_manual.html 颜色的获取可以参考:https://www.omicsclass.com/article/746 学习R语言基础与绘图:、R语言画图、R语言快速入门与提高
回答于 2019-11-08 09:36
检查GFF文件里面的mRNA ID是否对应一致。 参考这个:https://www.omicsclass.com/question/1941
回答于 2019-11-07 10:04
可能你的fasta文件不是标准的格式,你看看Chr2序列折叠之后长度是否一致;
回答于 2019-11-06 17:46
建立索引这一步很消耗内存,几百G以上内存的节点运行一下试试;
回答于 2019-11-06 14:46
运行时间和电脑的配置与基因的数量有关。如果电脑性能差,基因数量大,运行会慢一些,耐心等待。 2018年买的5000元的笔记本电脑,40000 vs 40000基因比对,大概运行一天时间,你可以参考一下;
回答于 2019-11-04 17:09