...想请教一下,在鉴定基因家族成员时,如果此类基因家族的典型结构域有2个比如说A,B,我看文献上已鉴定出来的此类家族一般成员同时含A,B;或者含A;或者B就说明是属于这个基因家族的,但是鉴定的时候有时候会出现A,B,C这...
请问一下,多序列比对之后的结果可以直接用于构树,不进行gap删除吗?因为我的序列比较多(几百条),用gblocks提取的保守序列很短,不到10bp,太短了构不了树。
...享文件夹,但是我在输入ll查询文件后,发现红框里显示的东西在我电脑的文件夹里并没有。而且我在电脑文件创建一个测试文件后,在终端输入ll后也不会出现我所新创建的文件。 我是用的mac,按照组学课堂docker教学课程学习...
clinvar是一个开放的数据库,每个研究机构都可以向其提交数据,对于提交的信息,会有专家团队进行审核评级。对于数据库中的位点,根据注释信息的可靠性,分成了1到4个不同的星级,星级越高,可信度越高。 数据库地址...
在使用GATK CALL 变异的时候非常的慢,一天一个样本还没有结束,能不能分染色体CALL 变异?原先的命令如下: gatk --java-options "-Xmx50g" HaplotypeCaller -R $REF \ -I $workdir/3.map/result/${i}.sorted.dedup.bam \ -O ${i}.g.vcf.gz --max-alterna...
补充:使用的是是组学大讲堂的镜像,用singularity调用的,前一步顺利运行,这一步报错的 麻烦老师尽快再回答一下吧!!卡在这里好几天了!!! 或者能不运行mcmctree,用别的什么结果运行cafe5吗 —————————— ...
上次的问题不太明确,这次重新提问一下,谢谢各位老师! 读入cucumberprotein.fasta文件,创建一个id和序列对应的hash, 遍历该hash,并输出到一个新创建的文件中
用自己做的物种蛋白序列单独构建进化树和加入其他参考物种(包括自己的物种)构建进化树的分组不一样,希望老师看到能解答下这种情况该如何处理?
我想做圈图,由于在一个位置上有12个基因,直接用标签没办法显示,想问问老师,这个图中的红色标注的折线是怎么实现的呢?
请问画圈图时,用热图表示,用颜色的深浅来表示数值的大小,这个要怎么设置,图例要怎么设置。还有怎么改变其中一个圈的宽度