Future plans
Faster pairwise alignment function (for d = 1 only)
Removing the Needleman-Wunsch function

Using a fix d = 1 value allows to use aggressive sequence comparison
strategies. The new pairwise comparison function introduced here is
more than 40 times faster than the Needleman-Wunsch function. The
timing is now 0.75 s for the test dataset.
Functional version
Cleaning
Faster microvariant function
Slightly faster microvariant function (now returns a set)
Simplification of the microvariants function
Simplify the nucleotides tuple
Cleaning
Using tuples instead of strings does not speed up
Add a -b option for post-processing

Add a -b option to produce the input data for the swarm breaker python
scripts. Add a step that sorts the amplicons by decreasing
abundance. Additional code cleaning and commenting.
Cleaned version