Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,708,9790TC54.8% ‑0.3 / 44.5 31E259G (GAA→GGA) bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base C (8/9);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGAAATCGACGGTTTCCGATCGCA  >  NZ_CP009273/3708845‑3709084
                                                                                                                                      |                                                                                                         
cgcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTccac                                                                                                       <  2:378944/139‑1 (MQ=255)
cgcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTccac                                                                                                       <  1:315559/139‑1 (MQ=255)
cgcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCccac                                                                                                       <  2:421857/139‑1 (MQ=255)
cgcCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCccac                                                                                                       <  1:134720/139‑1 (MQ=255)
  ccAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCccac                                                                                                       >  2:346831/1‑137 (MQ=255)
                 tctGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcc                                                                                      >  1:252663/1‑139 (MQ=255)
                   tgctgcCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATccc                                                                                     >  2:281629/1‑138 (MQ=255)
                        cccATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCAt                                                                               >  2:81703/1‑139 (MQ=255)
                          cATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCccac                                                                                                       >  1:91200/1‑113 (MQ=255)
                          cATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCccac                                                                                                       <  2:91200/113‑1 (MQ=255)
                           aTTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAg                                                                            <  1:450552/139‑1 (MQ=255)
                              cggcggTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAgcg                                                                          <  1:303841/138‑1 (MQ=255)
                                                cGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                                                       <  2:409553/139‑1 (MQ=255)
                                                cGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                                                       <  2:107117/139‑1 (MQ=255)
                                                cGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                                                       <  2:235915/139‑1 (MQ=255)
                                                cGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                                                       <  2:207656/139‑1 (MQ=255)
                                                  tCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTccac                                                                                                       <  1:278117/89‑1 (MQ=255)
                                                  tCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTccac                                                                                                       >  2:278117/1‑89 (MQ=255)
                                                  tCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCg                                                     >  1:86301/1‑139 (MQ=255)
                                                        tcGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGc                                                 <  1:346831/137‑1 (MQ=255)
                                                          gCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGccacc                                             >  2:116987/1‑139 (MQ=255)
                                                                cAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTc                                        >  1:486893/1‑138 (MQ=255)
                                                                cAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGa                                                                                         >  1:248441/1‑89 (MQ=255)
                                                                cAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGa                                                                                         <  2:248441/89‑1 (MQ=255)
                                                                          gCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGc                             >  1:449952/1‑139 (MQ=255)
                                                                                  ggCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGaaa                     >  1:419982/1‑139 (MQ=255)
                                                                                      tAACCTGTTGTTCAACGCGCCAGAAACGCTTCGTCGTGCCCACCATTTTCCCCCAGCGTTTGGGCCATCCCCAGCCTCAGCGCAACATCGGTGCCCATGTGCGGTGCCCCCCACTCCATTTTATCGCCAAAGAAATCGa                 <  1:482419/139‑1 (MQ=255)
                                                                                             tggtgCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGAAATCGACGGTTTc          >  1:249427/1‑139 (MQ=255)
                                                                                                  cAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGAAATCGACGGTTTCCGAt      <  1:281629/138‑1 (MQ=255)
                                                                                                  cAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTCCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGAAATCGACGGTTTCCGAt      <  2:86301/138‑1 (MQ=255)
                                                                                                     cgcgCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGAAATCGACGGTTTCCGATCGCa  >  2:300551/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                         
CGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTTATCGCCAAAGAAATCGACGGTTTCCGATCGCA  >  NZ_CP009273/3708845‑3709084

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑