ExpansionHunter位点定义与重复估计 - SpiralKernel/bioinfo-wiki GitHub Wiki
ExpansionHunter位点定义与重复估计
更新:2026-10-15;状态:可用;验证范围:官方文档与方法核对;专业命令为模板,本轮未执行该流程。
ExpansionHunter按目标位点的序列结构估计重复,位点定义与证据范围决定可解释性。
最小catalog结构
官方格式使用JSON数组。以下只演示一个人工定义,没有配套BAM,不是可直接运行的数据集:
[
{
"LocusId": "toyRepeat",
"LocusStructure": "(CAG)*",
"ReferenceRegion": "chrToy:100-130",
"VariantType": "Repeat"
}
]
ReferenceRegion使用0起始半开坐标;LocusStructure表达位点结构。它与TRGT BED不是同一格式,不能仅改扩展名。
输入与模板
准备匹配的短读长比对、参考、catalog和适用样本信息;本轮未安装或执行:
ExpansionHunter --reads sample.bam --reference reference.fa --variant-catalog repeats.json --output-prefix sample.repeat
跨越、侧翼和重复内部reads能提供不同信息,短reads/片段长度可能限制估计;特殊RareRepeat设置需要相应依据,不能把所有目标都改成该类型来追求更长数值。输出重复次数、区间和失败状态按头部与版本核对。
阅读串联重复STRVNTR与中断、STR与重复扩增检测及提取;工具与许可以官方仓库为准。