超过500M的染色体,GATK 索引不支持; 1.解压g.vcf.gz 2. 解压后的vcf文件建立索引: gatk --java-options "-Xmx50g" IndexFeatureFile -I demo.g.vcf3.用解压后的g.vcf导入数据库 参考这里:https://www.omicsclass.com/question/7288
回答于 2024-05-22 10:39
按照报错信息更新 wsl,或者百度一下,或者购买我们的云服务器:https://study.omicsclass.com/index 省去安装的麻烦;
回答于 2024-05-21 17:38
没有找到成对read的比对,你这个转录组是不是单端测序,单端测序没法评估这个; 或者你hisat比对的时候,输入read1文件和read2 文件,错误的都是输入read1 文件或者都是输入的read2 文件,检查命令行是否错误;
回答于 2024-05-21 09:52
这个权重值没有统一的标准,可以根据自己的需要选择多个筛选阈值,分析hub基因,直到筛选合适的; cytoscape也有一些插件可以帮助我们筛选hub基因:MCODE
回答于 2024-05-21 09:34
DNA 那一行,是因为 属于DNA转座子但是没有子类; 其他的,比如LINE转座子,所有的都有子类;
回答于 2024-05-21 09:31
如果要去掉图例的黑框添加代码: kegg_point + theme(legend.key=element_blank())
回答于 2024-05-20 10:17