Consensus sequences from IRMA and read mapping/variant calling should be compared using global alignment.

[nhhaidee]

I think the consensus sequences from IRMA and read mapping/variant calling should be compared using global alignment. I think Edlib would be suitable for this task. It's fairly easy to use in a Python script (see https://github.com/CFIA-NCFAD/nf-ionampliseq/blob/76f9ea377156f508fc92b218f7663099e7ae5872/bin/edlib_align.py#L18). It'd be useful to report the differences if there are any.

We should have the number of differences between the IRMA and read mapping/variant calling consensus sequences and their closest ref seqs from BLAST and variant calling, respectively.

Originally posted by @peterk87 in #8 (comment)