这篇文章出现在我上高一时的语文书上（人教版2003版）。原作者是托马斯·刘易斯（Lewis Thomas），来自他的著作《细胞生命的礼赞——一个生物学观察者的手记》（1974年）（The Lives of a Cell: Notes of a Biology Watcher）。翻译者是湖南科学技术出版社《第一推动丛书》。该文章的英文原名叫做：The Music of This Sphere，中文译名：《这个世界的音乐》。难点在文章的最后一段，我当时的语文老师直言：怎么看都看不懂。
英文原文发表在著名的《新英格兰杂志》上（New England Journal of Medicine）。中文原文由于官方不公开，我只好从第三方随便搞一个。
下面我就把课本上的中文译文和英文原文对照起来一起展示。
GTF虽然比较古老，但确实有些软件（即使是最近的，例如SQANTI）会用到。转换的方法可能有各种各样，下面给出我实践过可行的方法：
gffread old.gff3 -E -F -T -o new.gtf
# -F preserve all GFF attributes (for non-exon features)
# -E show all warningsgffread old.gtf -E -F --keep-genes -o new.gff3可能有人会问GFF3转成GTF为什么不用gt gff3_to_gtf？主要问题在于这样子会把ID编号全部丢掉（我手头上的v1.5.8），所以不用了。
PharmGKB有一个隐藏得很深的基因下载页面，在这里可以下载获得各种excel格式的Haploid定义表、功能表、频率表等。
URL是：https://www.pharmgkb.org/page/cyp2a6RefMaterials，把其中的cyp2a6替换成自己要的基因即可。
全部的基因列表是：https://www.pharmgkb.org/page/pgxGeneRef
最近（2020年以来），Aspera（ascp）无法从NCBI下载SRA了，但是可以用ENA下载。现在以（SRR10609482）为例：
访问https://www.ebi.ac.uk/ena/browser/view/SRR10609482，获取准确链接如下：
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR106/082/SRR10609482/SRR10609482_1.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR106/082/SRR10609482/SRR10609482_2.fastq.gz可以看出中间出现了一个奇怪的三位数082。
很多情况下，人类参考基因组上的碱基是major allele。但也有相当多的情况下，它们是minor allele，有时甚至是rare allele。
BRCA2基因上有一个变异位点编号为rs169547，为非同义突变。全世界98%的人都是等位基因C，只有2%的人是等位基因T。
但是参考基因组上的等位基因是T（GRCh37和38都是），这纯属侥幸。
在参考基因组中可能会见到诸如 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确认一下。