根据gff信息提取基因组指定位置序列
#!/usr/bin/env python3 # -*- coding: utf-8 -*- import argparse parser = argparse.ArgumentParser(description='该脚本用于在基因组特定位置截取序列,需额外输入记录有截取序列信息的列表文件', add_help=False, usage='\np…
PCoA+boxplot
# rm(list=ls()) # # If devtools package is not installed, first install it # #install.packages("devtools") # #devtools::install_github("ChiLiubio/microeco") # pacman::p_load(m…
16s_lac_bif.sh
qiime2.slurm #!/bin/bash #SBATCH -p batch #SBATCH -N 1 #SBATCH -n 30 ./qiime2.sh #bac: V3,17, 341F, CCTACGGGNGGCWGCAG # V4, 20, 806R, GGACTACHVGGGTWTCTAAT #bif:greol:23,22: TC…
meta.sh
meta.slurm #!/bin/bash #SBATCH -p batch #SBATCH -n 30 # sh meta.sh ./meta.sh #!/bin/bash # HEADER - Do Not Modify! set -e shopt -s expand_aliases export LC_ALL=C ######## expo…
docker项目
#######npm sudo -i mkdir -p /root/data/docker_data/npm cd /root/data/docker_data/npm vim docker-compose.yml version: '3' services: app: image: 'jc21/nginx-proxy-manager:latest…
PCOA
#载入绘图所需的包 library(vegan) library(ape) library(ggplot2) library(ggrepel) #分组数较少、组内生物学重复较多 data <- read.csv("otu.txt", head=TRUE,sep="\t",row.names = 1) #载入分组文件 groups <- …
相关性热图
library(readxl) library(ggcorrplot) library(ggthemes) data <- read_xlsx("correlation2.xlsx",sheet = 2) data<-as.data.frame(data) rownames(data)<-data$id data<-data…
PCA
pacman::p_load(tidyverse,ggrepel,FactoMineR,magrittr,factoextra,RColorBrewer) df <- read_tsv("F3.xls") pca <- df %>% column_to_rownames(var="Sample_id") %>% select(-Subt…
爬取ncbi页面数据
# install.packages('RSelenium') #java -Dwebdriver.chrome.driver="chromedriver.exe" -jar selenium-server-standalone-4.0.0-alpha-2.jar -port 4444 library(readxl) library(RSeleni…