Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273645,1890TG54.5% 26.0 / 27.6 22H253P (CAC→CCC) citEcitrate (pro‑3S)‑lyase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base G (12/0);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 1.41e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAG  >  NZ_CP009273/645100‑645277
                                                                                         |                                                                                        
ccGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGt                                                                                          <  1:557243/90‑1 (MQ=255)
ccGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGt                                                                                          <  1:697050/90‑1 (MQ=255)
  ggACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGGGa                                                                                        >  1:44086/1‑90 (MQ=255)
  ggACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGGGa                                                                                        >  1:487446/1‑90 (MQ=255)
      cGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGGGGTccc                                                                                    >  2:265328/1‑89 (MQ=255)
         ccaccaTCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCActt                                                                                 <  2:48738/90‑1 (MQ=255)
                 ttGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGGGATCCACTTCTTTCTgg                                                                         >  1:831399/1‑90 (MQ=255)
                     cGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCg                                                                     <  2:1072938/90‑1 (MQ=255)
                      gTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCgg                                                                    <  2:1230920/90‑1 (MQ=255)
                        tCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGTCCACTTCTTTTTGGGGGGGTg                                                                  >  2:210757/1‑90 (MQ=255)
                        tCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGATCCCCTTCTTTTTGGGGCgggg                                                                  >  1:720406/1‑88 (MQ=255)
                              aaaCCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTaga                                                            <  2:504998/90‑1 (MQ=255)
                                aCCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGTCCCCCTTTTTTTTGGGCGGGGCGGAGAgg                                                          >  2:990164/1‑90 (MQ=255)
                                aCCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGTCCACTTTTTTTTTGGGCGGGGGGGAGAgg                                                          >  2:119139/1‑90 (MQ=255)
                                   aCGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGTCCCCCTCTTTTTTGGGCGGGGGGGAGAGGGtg                                                       >  2:593401/1‑90 (MQ=255)
                                       cGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGGGGGGTCCACCTCTTTTTTGGGCGGGGGGGAGAGGGtgggga                                                   >  2:429476/1‑86 (MQ=255)
                                          ggCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGGGGTTCCCCTCTTTTTGGGGGGGGGCGTTGAGGGTGTggagca                                                >  2:629744/1‑90 (MQ=255)
                                                  cgAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATc                                        <  1:752496/90‑1 (MQ=255)
                                                       ggcggcTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGAcg                                   <  2:1089881/90‑1 (MQ=255)
                                                                    gCTTCTACGACGCGGCGGGCGGGATCCACTTCTTTCTGGGTCGGGGCGTAGAGGGTGTGCAGCAGATCAATCTGACGCGGGGTGATCAGc                      >  2:866322/1‑90 (MQ=255)
                                                                             acgCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCg             >  2:724698/1‑90 (MQ=255)
                                                                             acgCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCg             >  1:291712/1‑90 (MQ=255)
                                                                               gcggcgGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTc           <  1:163148/90‑1 (MQ=255)
                                                                                        gTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAg  <  1:657167/90‑1 (MQ=255)
                                                                                         |                                                                                        
CCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAG  >  NZ_CP009273/645100‑645277

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