天使不设防

摘要： 基本思路是，将key对应的value设置为list，将对应的值append进去。 示例： f=open("a1.txt") ha={} for i in f: i=i.strip().split() print(i[0],i[1]) for k in i[1:]: ha.setdefault(i[0 阅读全文

摘要： vim编辑shell脚本自动插入文件头部信息，将下面的代码写入home目录xia .vimrc 文件即可。 shell 文件头： 1 autocmd BufNewFile *.sh exec ":call Setcomment()" 2 func Setcomment() 3 call append 阅读全文

摘要： Jvarkit : Java utilities for Bioinformatics ：一个java写的生物信息工具包：http://lindenb.github.io/jvarkit/ 阅读全文

摘要： github地址：https://github.com/agordon/datamash 阅读全文

摘要： linux 下fasta序列按碱基显色 1 use strict; 2 use warnings; 3 4 my $black = "\033[0;30m"; 5 my $dary_gray = "\033[1;30m"; 6 my $light_gray = "\033[0;37m"; 7 my 阅读全文

摘要： 地址： https://github.com/jespermaag/gganatogram 阅读全文

摘要： 切片ik通过索引访问，然后为每个map分配内存； 切片jk通过获得切片内每个元素的拷贝来分配内存，并未成功为切片内每个map分配内存，使用时赋值也就失败了 1 package main 2 3 import fmt 4 5 func main(){ 6 ik := make([]map[int]in 阅读全文

摘要： GATK4 检测的SNP标记，有些位点会在检测过程中完成 phasing，在后续做基因型填充的时候有坑。 GATK4 phasing 结果的缺失位点不是 ./. 也不是 .|. 而是直接变成一个单独的点；下图黄线标记出来的部分，上面是原始结果，下面是修改后结果。 基于此后续做 phasing 才能顺 阅读全文

摘要： 打开多个文件时必备啊 1 use FileHandle; 2 3 my (%index,%fh); 4 # 创建句柄,这里利用gzip 压缩下 5 foreach my $k(keys %index){ 6 $fh{$k}=FileHandle->new("| gzip -c > ${k}_1.tx 阅读全文

摘要： substr EXPR,OFFSET,LENGTH,REPLACEMENT substr EXPR,OFFSET,LENGTH substr EXPR,OFFSET Extracts a substring out of EXPR and returns it. First character is 阅读全文