有时图例较多,用ggplot2画图时候会出现图例占据图的半壁江山,想让图例排成一列,解决方法如下: guides(fill = guide_legend( ncol = 1,...
有时图例较多,用ggplot2画图时候会出现图例占据图的半壁江山,想让图例排成一列,解决方法如下: guides(fill = guide_legend( ncol = 1,...
操作分类单元(Operational Taxonomic Units,OTU)是在系统发生学研究或群体遗传学研究中,为便于分析,人为给某一个分类单元(品系,种,属,分组等)设...
人类基因组中存在大量的重复区域,给定一段序列,如何知道这段序列在基因组中是否具有唯一比对,或者具有多个比对? 一种方法是通过bwa,通过AS和XS 标签判断是否有次优比对;但...
1.为什么要在上机测序之前做PCR扩增? 测到的reads数取决于flowcell上长了多少cluster(1个模板链扩增出一簇分子),有多少cluster,理论上就得到多少...
第一步:编辑bashrc 编辑 ~/.bashrc 文件,第一个变量$LOCAL_APP根据自己情况修改,后面直接复制。 #local app path LOCAL_APP=...
1. 实现目录下大批量文件的删除 当某个目录下存在大批量文件的时候,如果用rm 命令进行删除,会很慢; 解决方法: 先创建一个空目录 mkdir Blank 然后 rsync...
1. 按照时间查询文件 先touch 两个时间标记文件 touch -t 07120800 start #####07120800 代表的意思(月,日,时,分) tou...
之前下载SRA的数据,使用 prefetch -v SRA --max-size 100G 下载数据,经常会出现下到一半就断网的情况;后来尝试使用wget下载,但...
1. awk命令实现文本内容行列转化 awk -F '\t' '{for(i=1;i<=NF;i=i+1){a[NR,i]=$i}}END{for(j=1;j<=NF;j++...