Lib.rs

count_constant_sites

给定一个包含核苷酸多重序列比对的FASTA文件，此工具计算比对中恒定位点的数量。输出是适合用于IQTREE的 -fconst 的行，即用逗号分隔的4个数字，表示A、C、G和T的数量。

恒定位点是指整个比对列都是单个核苷酸的位置。此工具不区分大小写。仅考虑A、C、T和G（即不考虑间隙和模糊核苷酸）。

待办事项

• 扩展以支持蛋白质字母表