南方都市报
k-mer分类结果适合做快速筛查。Kraken类工具、Centrifuge类工具或⭐基于Mash的相似性分析可以先给出候选物种,但分类数据库的构建版本会改变结果,因此必须记录数据库来源、版本和过滤参数,并用独立比对复核。
短序列的物种判断需要考虑片段长度和保守程度。保守基因的一小段序列可能同时匹配多个哺乳动物,只有包含足够特异位点、且正反向比对结果一致时,结论才更可信。
质控工具的输出应保存为独立报告。fastp、FastQC等工具可以帮助发现质量问题,但它们只能说明序列是否适合分析,不能直接回答样本来自人、猪还是狗。
人or猪or狗的物种判定必须建立在可用读段之上,低质量碱基、接头残留和过短序列会明显增加错误匹配。短读段分析可以先检查每个样本的总读段数、平均质量、N碱基比例、接头污染和重复率⭐,再决定是否进行剪切。
人or猪or狗的正式判断应同时比较三个候选参考基因组,而不是只把数据比对到其中一个物种。把读段分别比对到人、猪、狗参考序列后,需要比较高质量唯一比对读段的比例、覆盖的染色体范围、平均覆盖深度、比对质量值和错配分布。
混合样本判断应同时检查主物种信号和次要物种信号。若人、猪、狗三个参考基因组都出现明显的高质量唯一比对读段,且次要信号分布在多个独立染色体,首先应考虑混样、交叉污✨染、样本标🎆签错误或参考数据库偏差。