Dna.filter_variants_consecutive

Dna.filter_variants_consecutive(proximity=[25, 50, 100, 200])

Remove nearby variants.

Remove the variants that are close to each other on the same amplicon. This is likely primer misalignment.

Parameters

proximitylist-like (int)

If i + 1 variants are within proximity[i], then the variants are removed.

Returns

keepnp.ndarray

Variants that are to be kept, i.e., variants that are close to each other are discarded.

Notes

Assumes the variants are in order.

Examples

• If >1 variants are within 25 bases of each other, remove.

• If >2 variants are within 50 bases of each other, remove.

• If >3 variants are within 100 bases of each other, remove.

• If >4 variants are within 250 bases of each other, remove.