天使不设防

摘要： C语言小练习：计算非压缩fastq格式的GC含量 1 #include <stdio.h> 2 #include <stdlib.h> 3 #include <string.h> 4 #define buff 1024 5 6 typedef unsigned long long int u_llo 阅读全文

摘要： http://methodspopgen.com/methods-to-infer-populations-history/ 阅读全文

摘要： 使用copy，直接改变原片的值，而不是先创建一个副本。 阅读全文

摘要： getline() 函数无论一行多长，动态分配内存读入行 1 #include <stdio.h> 2 #include <stdlib.h> 3 #include <string.h> 4 5 int main(int argc,const char *argv[]) 6 { 7 FILE *fp 阅读全文

摘要： 需要安装包：ggforce，下面以R自带数据做局部放大演示。 require(ggplot2) require(ggforce) require(reshape2) data(CO2) co2<-melt(CO2,id.vars=c("Plant","Type","Treatment"),varia 阅读全文

摘要： set.seed(1) h1<-hist(rnorm(1000,100,5)) h2<-hist(rnorm(1000,99,5)) plot(h2,col=rgb(255,0,0,50,maxColorValue =255),border = NA) plot(h1,col=rgb(225,225 阅读全文

摘要： 一、芯片数据 此次拿到的illumina芯片数据并不是原始的数据，已经经过GenomeStudio软件处理成了finalreport文件，格式如下： 之前没处理过芯片数据，对于这种编码模式（Forward，top AB）的基因型数据很疑惑，查了很多资料，收效甚微。看过建明大神对芯片这块儿的介绍，发现 阅读全文

摘要： 实现：eval 1 a="indv1" 2 indv1="Sus1" 3 4 eval tmp='$'$a 5 echo $tmp //这里 echo 返回值为Sus1 阅读全文

摘要： 去除基因组序列中的未定位的scaffold、Contig序列和线粒体序，主要针对NCBI提供refseq基因组序列，组装到染色体级别的物种基本都通用。将所有碱基统一成大写字母，并计算每条染色体长度，每80个字符换行。处理脚本如下： 1 use strict; 2 open A,"$ARGV[0]"; 3 open B,">$ARGV[1]"; 4 open C,">$ARGV[2]... 阅读全文

摘要： 目标如题，有多个fasta文件和一个文件名列表，将文件名列表中包含的文件匹配出来并提取第一条序列合并成一个fa文件。这个采用perl实现，用法和代码如下： 1 #!/usr/bin/perl -w 2 use strict; 3 4 sub usage{ 5 die "usage:perl $0 \n",unless(@ARGV==3); 6 } 7 usage(); 8 ... 阅读全文