Modules§
Structs§
Functions§
- align_
single_ query_ to_ targets - align_
worker - build_
aligner - build_
aligner_ v2 - does not work now
- build_
bam_ record_ from_ mapping - convert_
mapping_ cigar_ to_ record_ cigar - if reverse alignment, the cigar will be reversed!
- hits2records
- mapping2str
- query_
seq_ sender - set_
primary_ alignment_ according_ 2_ align_ score - set_
primary_ supp_ alignment_ according_ 2_ align_ score - 重新设置 primary 和 supplementary。其主要在 多个参考基因序列场景使用 由于对于多个参考序列场景,每个参考序列都构建了一个aligner,所以我们会得到之多 N 个 primary。N表示参考序列个数。 primary设置逻辑:primary中,比对分数最高的作为最终的primary, 其余的primary基于 query 的 ovlp 来确定自己是 supplementary 还是 secondary supplementary设置逻辑:所有hits重新排列,按照 primary 在第一位,supplementary 按照 比对得分位列2~?。secondary排到最后 依次判断 当前 range 和 之前的 supplementary range的 ovlp。如果 ovlp>0.5 当前 supp 就降级为 secondary
- targets_
to_ targetsidx - {“target_name”: (idx, length)}