文章 linux中sed在文件部分列尾添加字符

在文件中包含>字符的行尾添加-mRNA。 sed -i '/>/ s/$/-mRNA/' pep.fa 将文件中mRMA和；之间的字符删除。 sed 's/mRNA\(.*;\)/mRNA;/g' gene.gff >gene_new.gff

文章 perl 读取文件中的科学计数数字

Perl读取文件中的科学计数法数字是是以字符串读取的，需要将其转化为小数，方法如下 use Math::BigFloat;my  $i = new Math::BigFloat '1.931533e-01'; 这里用到了Math::BigFloat这个包，就可以将其转换为小数进行计算了。 更多perl语言知识...

文章 trimAL | 多序列比对中剪切掉不保守部分的参数

...信分析时，经常会做一下蛋白多序列比对。对于比对结果中，一些不保守的列，我们可以选择剪切掉。再进行之后的分析，比如构建进化树，进化树的构建会容易很多。 剪切序列时有一个非常好用的软件是trimAL: http://trimal.cgeno...

问题 老师，您好，我按照课程安装doekcr，下载了群体遗传学的镜像，在启动镜像的时候总是显示错误，请问老师该如何解决呢？

...在启动镜像的时候执行指令：docker run --rm -it -m 4G --cpus 1 -v /home/share/pop-evol:work omicsclass/pop-evol-gwas:latest 总是显示Error response from daemon：failed to create task for container: failed to create shim task: OCI runtime create failed: runc create failed: unable to ...

文章 bugbase指定表型预测

...内存，所以需要给docker更多的内存) docker run --rm -m 25 -v D:\project\bugbase:/work -it omicsclass/ampliseq-q1:v1.2 2.生成指定表型的预测文件（运行默认为16S数据） make.user.table.r -i costom_pathway/ 参数说明 -i    准备好指定表型的文件...

问题 泛基因家族分析PAV_gene.list中基因数量，少于gene.pep.fa中序列的数目

老师好！PAV结果中，为什么PAV_gene.list中基因数量（204个）少于02.gene.family/check下gene.pep.fa中序列的数目（218个)？ 如何判断少的这14个属于PAV_gene.list中Fam_N列下的哪个基因？ 或者要重新分析？

问题 老师您好！

老师您好！我用你们提供代码在GDC里面下载数据时候出现错误，错误代码如下： GDCdownload(query = query, +             directory = DataDirectory,files.per.chunk=6, method='client') Downloading data for project TCGA-CHOL 试开URL’https://gdc.cancer.gov/files/...

文章 Linux中tree命令显示目录名字遇到中文乱码问题

Linux中tree命令显示中文问题 有中文名字的文件，使用tree显示会有\xxx的内容，如： 网上有说安装新的tree可以解决的这个问题。后来发现使用tree的-N选项就可以正常显示中文文件：

问题 生存分析中的生存曲线如何绘制？

问题 circos中热图部分画不出来

问题 挂载过程中遇到以下问题需要怎么处理

文章 vcftools去除或保留vcf文件中的样品

vcftools可以去除或保留vcf文件中的样品，用到的选项为 --indv 和 --remove-indv  ，指定要从vcf文件中保留或删除的样品。 可以多次使用此选项来指定多个样品。 具体用法如下： 下图为原始vcf文件。 只保留1和10号两个样品，执...

问题 拟南芥UVR8有7个RCC1蛋白保守结构域，在之前一些文献的图片中可以直观的看到，为什么现在输入序列，结果却是长长的一条，怎么才能做出文献中那种图？求解，谢谢

问题 新手小白，刚开始学习docker，想问一下，做泛基因家族分析最开始，怎么将老师服务器中使用的加载环境变量用到的env.sh（还有老师课程中一些其他必须要用到脚本）粘贴到docker中呢

问题 qRT-PCR的结果中Ct值在25-32之间

如题，qRT-PCR的结果中Ct值在25-32之间，请问这样的数据能用吗？