Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,991,5290AT52.6% 12.7 / 25.5 19T63S (ACC→TCC) uvrDDNA helicase II
Reads supporting (aligned to +/- strand):  ref base A (0/9);  new base T (10/0);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGC  >  NZ_CP009273/3991440‑3991606
                                                                                         |                                                                             
gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa                                                                               <  2:612377/90‑1 (MQ=255)
  ggTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGACGTTTTcc                                                                             >  1:117233/1‑90 (MQ=255)
          gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGCGGTTTTCAACaaaac                                                                     >  1:951154/1‑88 (MQ=255)
          gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGAGGGTTTTCAACaaaaa                                                                     >  1:461070/1‑88 (MQ=255)
          gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGACGGTTTTCAACaaaaa                                                                     >  1:96261/1‑88 (MQ=255)
            gcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGCGGTTTCCaaaaaaaagg                                                                   >  2:230886/1‑86 (MQ=255)
                     tGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAAcagcgggggggg                                                          >  1:192690/1‑87 (MQ=255)
                        aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc                                                       <  2:699530/90‑1 (MQ=255)
                        aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc                                                       <  2:297253/90‑1 (MQ=255)
                             aTCGCCTGGTTGCTGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGGGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCAt                                                  <  2:679042/90‑1 (MQ=255)
                                gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGTTTTCCAACAAAACGGCGGCGGGGGGGCGTCATCGt                                               >  1:299799/1‑90 (MQ=255)
                                    ggTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACGGGGGGGGGGGGGGGGCCTCGTATCg                                           >  1:423589/1‑90 (MQ=255)
                                        gatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCa                                       <  2:210374/90‑1 (MQ=255)
                                          tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAAc                                     <  2:252445/90‑1 (MQ=255)
                                                    aaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc                           <  1:175020/90‑1 (MQ=255)
                                                    aaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc                           <  2:459040/90‑1 (MQ=255)
                                                           tCGCCATACTCGATTATGGCGGTGACGTTTTCCAACAAAACGGCGGCGGGGGTGGGGCATCGTATCGGGGGACAGATGGGGGCGAGcccg                    >  1:943300/1‑88 (MQ=255)
                                                              ccATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGc                 <  2:555894/90‑1 (MQ=255)
                                                                        ttATGGCGGTGACGTTTTCCAACAAAACGGGGGGGGGGGTGGGGCCTCGTATCGGGGGACCGATGGGGACGAGCCCGGGCGGGGTGTggg       >  1:888293/1‑90 (MQ=255)
                                                                             gCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGc  <  2:173101/90‑1 (MQ=255)
                                                                                         |                                                                             
GTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGC  >  NZ_CP009273/3991440‑3991606

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