perl data.pl id1.txt id2.txt 现有一文件,有两列,是按照第一列内容进行展示的,如下图:A green;red;blueB yellow;greenC red;blue;grey;pink;blankD green;grey现在想将文件...
利用 NAS 对实验室数据进行存储备份
给大家介绍一个R包用于绘制Cladogram,直接附上代码: if (!requireNamespace("BiocManager", quietly = TRUE)) { install.packages("BiocManager")}BiocManager::install("microbiomeMark...
Rocky 9.1修改IP地址
通常进行SNP calling的之前我们需要对比对的bam文件进行排序、去重等工作,在这里介绍几种sort bam的方法。 1. 按照比对位置进行排序 sambamba比samtools快很多,samtools和sambamba默认按...
结构变异结果文件格式转换工具——SURVIVOR
转换成dummy dummy-variables
很多时候简单的try函数就能够增加脚本的鲁棒性,随着脚本逐渐复杂,一些较高级的错误处理函数就必不可少了,我们今天来介绍几个R中的 1. tryCatch() 函数 功能: tryCatch() 函数是一种更复杂...
叶绿体或线粒体重复序列圈图绘制
awk使用NR==FNR同时读取两个文件
ssh登陆报错
最近在使用AI软件对iTol生成的pdf文件进行修改。当AI导入文件后弹出提示窗口如下: “缺少字体。将使用替换字体显示受影响的文本。” 一般情况下可以忽略提示。 但这次的文件在忽略提示直接导...
介绍如何在大型生物数据库NCBI上下载数据
解决办法:https://www.omicsclass.com/article/2194 这里还有一个需要做的步骤,就是先将raw.vcf.tmp.idx文件找到并删除,一般情况下这个文件是空的。 然后再在命令行中添加-disable_auto_in...
细胞焦亡、铁死亡及铜死亡浅析
在这个科研内卷的年代,你没有好设备,没有好材料,没有充足经费,怎么立足?毕业、就业、发文、评职称,靠什么来实现?留给你破局的选择不多!!为此,我们组学大讲堂决定干一票大的!我们将永...
在做SV、CNV变异检测时,会使用SURVIVOR将不用样品的vcf文件合并到一起,命令如下: ls *.max.vcf >sample_files && /share/work/biosoft/SURVIVOR/SURVIVOR/Debug/SURVIVOR merge...
在生信分析时,经常会做一下蛋白多序列比对。对于比对结果中,一些不保守的列,我们可以选择剪切掉。再进行之后的分析,比如构建进化树,进化树的构建会容易很多。 剪切序列时有一个非常好用的...
InterProScan是InterPro数据库提供给用户的一个软件包,通过集成数据库鉴定蛋白质功能、所属的家族及其包含的结构域和位点