Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,091,4350TG57.9% 14.3 / 19.7 19intergenic (+77/‑183)sodA/kdgTsuperoxide dismutase [Mn]/2‑keto‑3‑deoxygluconate transporter
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (11/0);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAT  >  NZ_CP009273/4091382‑4091521
                                                     |                                                                                      
ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAAATTTAAAAAAAAAAAAAAGaac                                                    >  2:365367/1‑88 (MQ=255)
 gCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATcc                                                   >  2:16493/1‑90 (MQ=255)
   gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGa                                                 >  2:229414/1‑90 (MQ=255)
     ataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGata                                               <  1:212882/90‑1 (MQ=255)
        ggCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAAAGATCCGataaaa                                            >  1:458224/1‑90 (MQ=255)
        ggCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGgtaaaa                                            >  2:532534/1‑90 (MQ=255)
          ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGTTTTTTTGGGATGAATTTTAAAATAAAAAAAATGATCCGataaaaaa                                          >  1:529329/1‑89 (MQ=255)
          ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGcaaaaaaa                                          >  2:170998/1‑89 (MQ=255)
           cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTGTGATCAATTTTAAAAAAAAAAAAATGATCCGataaaaaaa                                         >  1:358418/1‑88 (MQ=255)
            gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGataaaaataa                                        >  1:288270/1‑90 (MQ=255)
             cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGataaaaaaaaa                                       >  2:32046/1‑90 (MQ=255)
               tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc                                     <  1:243533/90‑1 (MQ=255)
                atCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACa                                    <  2:203281/90‑1 (MQ=255)
                     gcttaaaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGtt                               <  2:411354/86‑1 (MQ=255)
                         aaaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAATAAAAAAACGTTTCaa                           >  2:243533/1‑90 (MQ=255)
                                 aCCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtgtg                   <  2:169784/90‑1 (MQ=255)
                                     tCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTAATTGATGTGGttt               >  2:524597/1‑90 (MQ=255)
                                        ccAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGa            <  2:307660/90‑1 (MQ=255)
                                          aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATc          <  1:92895/90‑1 (MQ=255)
                                                  tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAt  <  2:413455/90‑1 (MQ=255)
                                                     |                                                                                      
GGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAT  >  NZ_CP009273/4091382‑4091521

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