若想要从sam或bam文件中提取指定区域内的reads,可以使用samtools和bedtools来实现。 首先准备一个区域信息文件。region.bed #第一列为染色体ID,第二三列分别为起始终止位置 ...
Bray-Curtis距离是以该统计指标的提出者J. Roger Bray和John T. Curtis的名字命名的主要基于OTUs的计数统计比较两个群落微生物的组成差异。与unifrac距离包含的信息完全不一样相比于jaccard距...
叶绿体是绿色植物进行光合作用的细胞器,具有合成蛋白质、淀粉、色素等功能,普遍存在于真核自养生物中,尤其是藻类和陆生植物,其基因组可自主遗传。叶绿体基因组序列对于研究植物物种起源,进化演变及不同物种间的亲缘关系等具有重要意义。
现如今,基于高通量测序得到的SNP信息无疑是海量的,然而SNP标记的验证却并不容易实现。其最好的验证方式是一代sanger测序,不过价格较高,而且通量低。如果SNP刚好位于酶切位点上,就可以将其转化为CAPS标记,利用跑胶的方式对其进行验证,既便宜又好用。
lefse分析
Plant CARE 地址:http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ PLACE 地址:https://sogo.dna.affrc.go.jp/cgi-bin/sogo.cgi?lang=en&pj=640&action=page&pag...
筛选内参基因
E value——E(Expect)值,表示随机匹配的可能性。E值越大,随机匹配的可能性也越大。E值接近零或为零时,具本上就是完全匹配了 E-value - this tells you the probability of getting this sco...
在linux上利用awk对数据进行筛选,并输出
大牛施一公:如何提高英文的科研写作能力
R 在线安装包,设置全局镜像,以免镜像不好链接中断 local({r <- getOption("repos") r["CRAN"] <- "http://mirrors.tuna.tsinghua.edu.cn/CRAN/" options(repos=r)}) 或...
Library preparation and Illumina sequencing Genomic DNA from each sample was prepared for libraries construction. The library preparations were sequenced by using an Illumina Hi...
今天有学员问 MCScanX 安装问题,我试了下,果然有错误。 make报错: g++ struct.cc mcscan.cc read_data.cc out_utils.cc dagchainer.cc msa.cc permutation.cc -o MCScanX msa.cc: In fun...
利用ggplot2绘制箱线图,并进行分面处理。
解决Bioconductor安装软件比较慢的问题
PCA pcoa分析取舍
转录组分析的3种模式及分析流程
利用EXCEl 绘制热图
查看read.table参数
MEGAN分类柱状图解释