跳转至

alnkit long locate

为每个类构造查询序列(断点旁 50 bp 锚点 + 全部外源),与类组装比对后逐条 contig 出图,并把锚点/外源在 contig 上的落点还原成参考基因组坐标。

用法

alnkit long locate [全局选项] [--regions …] [--sites …] [--reference-fasta …] [--exog …] [--assembly-dir …] [--outdir …] [--threads …] [--flank-bp …]

参数

参数 必填 默认 说明
--regions 是 — breakpoint_classes.tsv
--sites 是 — breakpoint_sites.tsv
--reference-fasta 是 — 参考基因组 fasta(需 .fai)
--exog 是 — 外源 fasta,可重复
--assembly-dir 是 — 06_assembly
--outdir 是 — 输出目录
--threads 否 12 minimap2 线程数(默认 12)
--flank-bp 否 50 断点锚点长度(默认 50 bp)

全局选项(每个子命令都有)

参数 必填 默认 说明
--log-level LEVEL 否 INFO 日志级别:DEBUG / INFO / WARNING / ERROR
--quiet 否 关 只输出错误(等价 --log-level ERROR)

例子

alnkit long locate \
          --regions 05_breakpoints/breakpoint_classes.tsv \
          --sites 05_breakpoints/breakpoint_sites.tsv \
          --reference-fasta reference.fa \
          --exog exog/insertA.fa --exog exog/insertB.fa \
          --assembly-dir 06_assembly --outdir 07_ribbon_positions

退出码

0 成功(可有警告)/ 1 运行期错误 / 2 参数或前置校验失败。