Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,527,9260AC53.3% 5.7 / 14.4 15G219G (GGT→GGGnanMN‑acetylneuraminate epimerase
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base C (0/8);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGA  >  NZ_CP009273/4527841‑4528012
                                                                                     |                                                                                         
aaatattaCCGGTGAAAGCAAGTGCAAATACGGCATCAGCCCGCAATCCTGGTTTCGCTGCGCCATAAATAATCCAGGTTTTACAACCtt                                                                                      <  1:484648/87‑1 (MQ=255)
aaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATGAGCCAGGTTTTATCACCtt                                                                                      <  1:60221/90‑1 (MQ=255)
    attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                  >  2:364237/1‑90 (MQ=255)
                            tACGGCACCCGTTCGCAATCCTGTTTTGGCTTCGCCATTAATAAGCCAGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTcc                                                        <  2:160641/90‑1 (MQ=255)
                                              cccTGTTTTGGTTTCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                      <  2:163353/89‑1 (MQ=255)
                                               ccTTGTTTGGCCTCCCCCTTAAAAAACCCGGTTTTTACCCCTTTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                      <  1:12309/90‑1 (MQ=255)
                                                  ggTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                  >  2:175419/1‑90 (MQ=255)
                                                    tttGCCTTCCCCATAAATAAGCCAGGTTTTATCCCTTTTTTTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                <  2:79985/90‑1 (MQ=255)
                                                    tttGCCTCCCCCATAAAAAACCCGGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                <  2:151334/90‑1 (MQ=255)
                                                     ttGGCTTCCCCATTAAAAAGCCCGGTTTTTTCCCCTTTTACCCACCCCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                <  1:271031/90‑1 (MQ=255)
                                                       ggCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                             >  2:596647/1‑90 (MQ=255)
                                                            cGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                        >  2:148515/1‑90 (MQ=255)
                                                                aTTAAAAAGCCGGTTTTTTCCCCCTTTATTCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                    <  2:437568/90‑1 (MQ=255)
                                                                        gCCAGGTTTTTTCCCCTTTTTTCCCACCCCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg             <  1:7344/90‑1 (MQ=255)
                                                                                  aTCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa  >  2:154552/1‑90 (MQ=255)
                                                                                     |                                                                                         
AATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGA  >  NZ_CP009273/4527841‑4528012

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: