部署到服务器服务器端服务器端主要是构建目录接受和存放博客内容。 eval1234567mkdir myblogcd mybloggit init --bare myblog.gitvi myblog.git/hooks/post-receive## 输入下面一行内容，目录换成自己的git --work-tree=/home/hexo/blog --git-dir=/home/hexo/hexo.g

最近在分析水稻的转录组数据，用的参考基因组是我们小组自己组装的，没有用常见的那几个参考基因组做比对，这就导致一个问题，

1conda install -c hcc aspera-cli 密钥位置： 1~/miniconda3/envs/metagenome/etc/asperaweb_id_dsa.openssh 下载NCBI nr库的原始数据： 1sudo /usr/bin/ascp -v -k 1 -T -l 400m -i ~/miniconda3/envs/metagenome/etc/asperaweb_

12345678910111213141516# 创建环境conda create -n bioperlconda activate bioperlconda install -c bioconda perl-bioperl# 查找Bioperl位置find anaconda3/* -name "Seq.pm"# 输出/home/lixiang/miniconda3/envs/bioperl/li

12345678ggdraw(p1 + theme_half_open(12)) + draw_plot(inset, .45, .45, .5, .5) + draw_plot_label( c("A", "B"), c(0, 0.45), c(1, 0.95), size = 12 ) 💌lixiang117423@foxmail.com💌lixiang11

后面要处理大量数据，明显感觉到data.frame这种格式满足不了大数据要求了，索性把data.table学一下。

1dplyr::mutate(geneid = stringr::str_split(id, "-") %>% sapply("[", 3)) 💌lixiang117423@foxmail.com💌lixiang117423@gmail.com

从GitHub克隆项目 1git clone https://github.com/eggnogdb/eggnog-mapper.git 下载数据库放在data文件夹下并解压 开始使用 需要安装diamond，版本越新越好 python：需要使用python3 i：输入文件 out：输出文件 m：使用diamond比对 1python emapper.py -i pep.fa --ou

12345678910111213141516171819202122232425import osfiles = os.listdir("G:\\33Pan_download\\下载数据\\unzip\\")all_genome = open("G:\\33Pan_download\\下载数据\\all_genome.fa","w")for i in files: print(i)

项目地址：https://github.com/fatedier/frp/ 服务端配置 123456789101112131415161718[common]# frp监听的端口，默认是7000，可以改成其他的bind_addr = 0.0.0.0bind_port = 70000# 授权码，请改成更复杂的token = pl00000000000# frp管理后台端口，请按自己需求更改das

下载 1wget https://jaist.dl.sourceforge.net/project/subread/subread-2.0.2/subread-2.0.2-Linux-x86_64.tar.gz 将bin目录加入到环境变量即可。 使用案列 1234567/srv/software/subread/bin/featureCounts -p -t exon -g transcr

1awk -F"[\t=]" 'BEGIN{OFS="\t"}{if($3=="mRNA"){print $1,$2,"gene",$4,$5,$6,$7,$8,$9"="$10 ; print $1,$2,$3,$4,$5,$6,$7,$8,$9"="$10"Parent="$10;}else{print $0}}' all.gff3 > all.gff3.new 复制粘贴使用： 1aw