下载IGV和导入文件的方法我就不多说了，可以直接在windows平台下使用，就跟你操作QQ一样，自己摸索就好了！

著名芬兰运动员Eero Mäntyranta，他拿过七枚奥运奖牌。他的血红细胞远超正常人水平，甚至一度被奥组委误以为服用了禁药。后来经过研究发现，他的EPOR基因上的一个位点rs121918116，发生了一个G>A的变异，使得他的血氧含量达到了普通人的150%，所以他耐力惊人。

在snpPedia里面可以查看这个位点的信息：http://snpedia.com/index.php/Rs121918116 Continue reading →

前面我们在第8讲提到了公司给我的一个报告的统计表格，有人反映不会做。

本来应该只需要给我6亿条reads的（PE150测序，人30X），但是足足给了我8.9亿条！（但事实上很多paper发表的基因组高于60X的也不少）

表格里面提到了好几个概念，比如duplicate的reads，一般说的PCR造成的duplicate，在找变异的时候需要去除掉。然后是那些比对到了不同染色体的reads pair，虽然只有2.29% ，也是需要重点分析的。（前面我也讲了如何提取以及分别分析它们！） Continue reading →

看来本次直播我的基因组分析流程效果还不错，不少朋友跟着自己动手开始分析全基因组测试数据了，值得表扬。其中有好几个朋友留言向我反映公司给的统计报告里面的覆盖度的问题，如下图： Continue reading →

我博客里面有详细讲读原文查解去除PCR duplication的reads的原理和方法，还比较了samtools和picard这两个软件的区别，请点击阅看（仔细探究samtools的rmdup是如何行使去除PCR重复reads功能的），或者复制链接（http://www.bio-info-trainee.com/2003.html）到浏览器查看。 Continue reading →

分析multiple mapping的情况，还是先用公司已经提供的bam文件来吧，后面我会用自己比对得到的bam文件再重新分析一次。

公司给我的bam数据是把有效测序数据通过 BWA（Li H et al.）比对到参考基因组 (b37)，这是目前使用率很高的一个软件。在比对过程中，如果一个或一对 read(s) 在基因上可以有多个比对位置，BWA 的处理策略是从中选择一个最好的，如果有两个或以上最好的比对位置，则从中随机选择一个。这种多重比对（multiple hits）的处理对 SNP、INDEL 以及 CNV 等的检测有重要影响。 Continue reading →

这类情况仅仅针对于双端测序数据，因为根据实验原理来看，对一个DNA片段，会把它的左右两端分别测序，但是测序仪器的测序长度有限，对本次实验来说，打断的DNA片段长度在350个碱基左右(这个长度只是一个分布，并不是真实值)，理论来说测序是左右各150，加起来也就300，也就是说DNA片段中间还有50个碱基是测不到的（当然，实际上是有可能测通的）。而对这个配对的reads来说，来自于同一个DNA片段，所以理论上它们应该比对到同一条染色体的。也还是基于对sam格式的文件的理解，前面我们提到了sam文件的第3，7列指明了该reads比对到哪条染色体，以及该reads的配对reads比对到了哪条染色体(如果比对到同一条染色体，那么第7列是=符号)。所以我们只需要写脚本来提取即可！ Continue reading →

之前我们说了比对上的数据，那么会有人想到有没有没有比对上的数据呢？

既然是从我的血液里面提取到的DNA进行测序的，那么理论上测序仪出来的所有测序reads都应该是我的，也应该都可以比对到人类的参考基因组，但是实际过程中的确存在着未必对上的数据。

现有人类基因组毕竟只是个参考，也许我有某些独特的DNA序列呢？而且也不一定测序数据就都是人类的，也许我血液里面会有那么些微不纯粹呢？也有可能是某些片段变异的太多了，超过了常见的比对软件的承受能力，可以试用SHRIMP这个软件来提取一下。当然，这不是本讲的重点。 Continue reading →

昨天，我们了解了一下SAM格式的比对结果，不知道大家理解的怎么样。但是全基因组测序数据实在是太大了，即使比对后把sam文件压缩成二进制的bam文件也还有55G（如何压缩转换可查看直播十二），如果完整的导入IGV查看会略微考验计算机配置。

Continue reading →

比对工具到现在已经多如牛毛了，见列表： https://en.wikipedia.org/wiki/List_of_sequence_alignment_software 。但是能被大多数人熟知的，就是bowtie和bwa（我们在十一讲中用的才是bwa），它们把测序数据比对到参考基因组之后，都会生成一个sam格式的文件。随后的大部分分析都是基于sam格式进行的分析，虽然Jimmy多次强调这些基础知识的重要性需要大家私下自学。但是由于这个sam文件实在是太重要了！！！所以，不得不亲自抽出一讲来说说它，后面也会基于此写十多篇文章： Continue reading →

我的全基因组数据还没拿到，而且还会推迟，简单说(tu)明(cao)一下原因（还好当初为了避免广告嫌疑一直没说是哪个公司负责测序，反正用的是illumina的hiseqX10这个测序啦，所以可以尽情的吐槽）。

心烦意乱的吐槽线 Continue reading →

大家久等了，Jimmy的测序数据还没有拿到手！但是，工欲善其事必先利其器！所以jimmy在等待自己基因组的这段时间里，准备好了自己计算资源！鉴于会有不少同志们会跟着直播来自己动手分析一个公共的全基因组重测序数据(参考数据下载地址详见直播一)，我在这里和大家分享一下计算机资源及软件的基本要求。 Continue reading →