Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273796,9950G.71.4% 73.6 / 52.8 42coding (781/1053 nt)ybhH4‑oxalomesaconate tautomerase
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base . (16/14);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGTCCCTT  >  NZ_CP009273/796857‑797118
                                                                                                                                          |                                                                                                                           
gcggAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCg                                                                                                                             <  2:295685/139‑1 (MQ=255)
         gATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTa                                                                                                                    <  1:752292/139‑1 (MQ=255)
              cGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGc                                                                                                              >  1:720396/1‑139 (MQ=255)
               gACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCg                                                                                                             <  2:771829/139‑1 (MQ=255)
                aCAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCg                                                                                                             >  1:289921/1‑139 (MQ=255)
                                   cATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAAt                                                                                                                  <  2:66455/115‑1 (MQ=255)
                                   cATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAAt                                                                                                                  >  1:66455/1‑115 (MQ=255)
                                   cATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCtt                                                                                         >  2:42154/1‑139 (MQ=255)
                                        aaTCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCa                                                                                     >  1:377868/1‑139 (MQ=255)
                                             aTCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATcgcgc                                                                               >  2:100541/1‑139 (MQ=255)
                                              tCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATcgcgcg                                                                              <  1:505841/139‑1 (MQ=255)
                                                            ggTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCgg                                                                >  1:802043/1‑139 (MQ=255)
                                                            ggTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCgg                                                                >  1:54617/1‑139 (MQ=255)
                                                              tAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtg                                                              <  1:467917/139‑1 (MQ=255)
                                                                      tGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgct                                                      <  2:300739/139‑1 (MQ=255)
                                                                        ggCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgct                                                      <  2:289921/138‑1 (MQ=255)
                                                                             aGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTc                                                 >  2:576257/1‑137 (MQ=255)
                                                                                  aTGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGCTAACCGGTGCTATTGCTATTTCCagtag                                           >  2:710412/1‑139 (MQ=255)
                                                                                  aTGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCagtagt                                          >  1:505632/1‑139 (MQ=255)
                                                                                  aTGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCagtag                                           <  1:596023/138‑1 (MQ=255)
                                                                                   tGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgcta                                                     <  1:595893/127‑1 (MQ=255)
                                                                                   tGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgcta                                                     >  2:595893/1‑127 (MQ=255)
                                                                                       aGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTtgtg                                      >  1:632152/1‑138 (MQ=255)
                                                                                                  tATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATcgcgcg                                                                              >  1:42991/1‑87 (MQ=255)
                                                                                                  tATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATcgcgcg                                                                              <  2:42991/87‑1 (MQ=255)
                                                                                                   aTCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCAc                         >  2:699018/1‑139 (MQ=255)
                                                                                                   aTCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCAc                         >  2:35/1‑139 (MQ=255)
                                                                                                           aCCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTcc                                                >  1:338535/1‑108 (MQ=255)
                                                                                                           aCCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTcc                                                <  2:338535/108‑1 (MQ=255)
                                                                                                               gtgCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgct                                                      <  1:478978/99‑1 (MQ=255)
                                                                                                               gtgCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgct                                                      >  2:478978/1‑99 (MQ=255)
                                                                                                               gtgCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGAc              >  1:344580/1‑138 (MQ=255)
                                                                                                                tgCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGAc              >  2:446960/1‑138 (MQ=255)
                                                                                                                tgCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgct                                                      >  2:295379/1‑97 (MQ=255)
                                                                                                                tgCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgct                                                      <  1:295379/97‑1 (MQ=255)
                                                                                                                tgCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACaa            <  1:699018/139‑1 (MQ=255)
                                                                                                                tgCTTATTTCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACaa            <  1:100541/139‑1 (MQ=255)
                                                                                                                  cTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACaa            <  2:377868/138‑1 (MQ=255)
                                                                                                                      tttCTCCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGt       <  2:505632/138‑1 (MQ=255)
                                                                                                                         ctcCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGTccc    <  2:614764/138‑1 (MQ=255)
                                                                                                                         ctcCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGTccc    <  2:748957/138‑1 (MQ=255)
                                                                                                                          tcCAGCGCAGAAAGGCGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGTCCCtt  >  2:177915/1‑139 (MQ=255)
                                                                                                                            cAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGTccc    <  1:446960/136‑1 (MQ=255)
                                                                                                                                          |                                                                                                                           
GCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCACCCGACAAATCGTCCCTT  >  NZ_CP009273/796857‑797118

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