天使不设防

摘要： 做数据比较的时候，由于同一个样本测序数据量不一致，需要抽取数据，控制数据量基本一致。 自己写脚本速度较慢，后面发现一个不错的工具：seqtk 原始数据抽取 如果只控制原始数据量一致，过滤低质量数据后直接使用seqtk (Version: 1.3-r106) 的子模块seq， 配合参数 -s 设定随机 阅读全文

摘要： 一、安装 我的数据库安装的是win版本，安装python后，直接命令行： 1 pip install mycli 即可。 二、使用 进入命令行后输入： 1 mycli -u root -p 888888 效果见下图 阅读全文

摘要： 今天在用R作图的时候发现配色不好看，在网上找到两个不错的取色网站，各有特色，推荐给有需要的人。 http://link.fobshanghai.com/rgbcolor.htm https://colorhunt.co/ http://www.peise.net/tools/web/ 阅读全文

摘要： 这里提供一个转置矩阵的perl脚本，R语言中的t()函数同样也能实现转置 1 use strict; 2 3 open A,"$ARGV[0]"; 4 5 my %ha; 6 my $head=<A>; 7 chomp($head); 8 my @tag=split /\s+/,$head; 9 1 阅读全文

摘要： 直接使用os模块的popen打开 import sys import os a=os.popen('/Soft/samtools-1.2/samtools flags '+sys.argv[1] ,'r') print(a.read()) 阅读全文