(一)参考基因组为minor allele
很多情况下,人类参考基因组上的碱基是major allele。但也有相当多的情况下,它们是minor allele,有时甚至是rare allele。
BRCA2基因上有一个变异位点编号为rs169547,为非同义突变。全世界98%的人都是等位基因C,只有2%的人是等位基因T。
但是参考基因组上的等位基因是T(GRCh37和38都是),这纯属侥幸。
很多情况下,人类参考基因组上的碱基是major allele。但也有相当多的情况下,它们是minor allele,有时甚至是rare allele。
BRCA2基因上有一个变异位点编号为rs169547,为非同义突变。全世界98%的人都是等位基因C,只有2%的人是等位基因T。
但是参考基因组上的等位基因是T(GRCh37和38都是),这纯属侥幸。
2023年更新:iOS 16 和 相应的 macOS 已经可以直接在日历里面添加节假日,不用人工输入URL。
iOS一般都会自带默认地区的节假日日历,例如国行版本自带中国大陆节假日日历(2022年已经有调休信息)。
下面介绍其他国家地区的日历订阅方法(如果只有iPhone:设置——邮件——账户——添加账户——其他——添加已订阅的日历。如果还有MacOS,则推荐在MacOS的日历——文件——新建日历订阅中添加,这样可以通过iCloud同步到iOS),参考来自Che’s Blog的文章,并实测有效:
在参考基因组中可能会见到诸如 chr6_apd_hap1
、chr1_gl000191_random
这样的序列,把它们去掉!否则bwa在做mapping时会认为这些区域的reads匹配不唯一,把mapping quality定为0,导致后续无法发现相应区域内的变异位点,造成假阴性!
不用什么特殊的软件或pipeline,直接使用samtools bedcov target.bed tumor.bam normal.bam
去计算每个目标区间的覆盖度,然后除一下看看比例就行(用LOG2转换一下更形象)
推荐使用COSMOS,速度比较快(5000X的大panel大约40min),无需复杂的参数,直接表格式结果,取size值高的即可(即supporting reads数目)。注意每个SV事件会列出两行。
--max-mnp-distance
参数(默认是1,可以增大比如20),但这不是万能的!拿到VCF结果之后根据坐标排序仔细核对!必要时用IGV确认一下。