sanger测序如果遇到indel杂合突变,两条allele的测序峰交织在一起会导致峰图几乎不可读。
Indigo是一个用于解决该问题的在线工具,可以直接上传abi文件,拆分两条allele。
但是,Indigo运行速度经常很慢,解决方法是使用它的命令行版本:tracy.
在参考基因组中可能会见到诸如 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确认一下。症状:用优盘安装系统,安装成功了,但是重启之后就进不去啦!不管你是否把优盘拔掉。
能用于定位oligo(或primer)的工具很多,比如BioEdit、BLASTN都可以,但是它们不能处理简并碱基的情况,只会把简并碱基当成错配。
现在列一个能够用于定位含简并碱基的oligo的工具:usearch中的search_oligodb命令。
命令格式为:
usearch -search_oligodb [your/target/seq.fa] -db [your/probes.fa] -alnout or -blast6out [your/result]
可用参数有:
-strand: plus or both
-maxdiffs: [int]