Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122411,3390CA37.0% 56.9 / 10.1 47F24L (TTC→TTA) kdkA_13‑deoxy‑D‑manno‑octulosonic acid kinase
Reads supporting (aligned to +/- strand):  ref base C (1/28);  new base A (17/0);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.03e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGG  >  ACB122/411211‑411468
                                                                                                                                |                                                                                                                                 
ggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTAATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGcccgc                                                                                                                             >  2:126689/1‑133 (MQ=255)
   gcgcAAGCCCCAGAAACTCTACGATCGCAAACAACGCTATGGGGAGGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCGCGACGCGGTGTTCGCCCTCCAgg                                                                                                                         <  7:278119/136‑1 (MQ=255)
     gcAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCTCCAggg                                                                                                                        >  2:70077/1‑135 (MQ=255)
     gcAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCCggg                                                                                                                        >  2:296777/1‑135 (MQ=255)
                      tACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                        <  8:363810/134‑1 (MQ=255)
                       ccGAGCGCAAGCACCGCTATGGGGATGCACTCTGATGGCGAGTGGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                        <  2:82937/132‑1 (MQ=255)
                       aCGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTGCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                        <  3:51745/133‑1 (MQ=255)
                        cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCCACAGGGCGAACGGCAgacccccgcc                                                                                                     >  4:379307/1‑133 (MQ=255)
                        cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCAGGCGGCGGACTTCGGCAGCCTCGACGAGGTGTTAGCCCACCCCGGCCAACCGCTgtcccccgcc                                                                                                     >  4:192664/1‑133 (MQ=255)
                        cGATCGCAAGCAACGCTATGGGGATGCACTCTGAGGGCGAGTGGGACACTGGCGCCGGGGTACGACCCCCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgac                                                                                                     <  1:135799/135‑1 (MQ=255)
                        cGAGCGCAAGAAACGCTTTGGGGATGCACTCTGATGGCGAGTTGGACCCTGGCGCGGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgac                                                                                                     <  1:89814/135‑1 (MQ=255)
                         gATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCAGGGCCCCCGGCTgaccccccccc                                                                                                    >  4:51745/1‑129 (MQ=255)
                              caagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCt                                                                                                <  1:296777/134‑1 (MQ=255)
                               aagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                              <  3:191773/135‑1 (MQ=255)
                                 gcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCCGCGCGAACGGCGGACCCCCCCCCCCCTCAgg                                                                                            >  2:5259/1‑134 (MQ=255)
                                   aaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCAGGGCGGACGGCGGACCCCCGCCCCGCGCAggg                                                                                           >  7:205127/1‑132 (MQ=255)
                                   aaCGCTATGGGGATGCACGCGGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCGGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGa                                                                                          <  4:43989/135‑1 (MQ=255)
                                      gCTTTGGGGATGCACTCTGATGGCGAGTGGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                      <  8:279242/136‑1 (MQ=255)
                                                ggCACTCTGATGGCGAGGTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATc                                                                             <  4:295497/134‑1 (MQ=255)
                                                  cACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCTCCCGGGCGAACGGCAGACCCGCGACCCGCACAGTGAAGTGGTACGCATCGAg                                                                          >  7:234585/1‑136 (MQ=255)
                                                   aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCCGGGCGCACGGCTGACCCCCCCCCCCCGCAGTGAAGTGGTACGCCTCGAg                                                                          >  4:229049/1‑135 (MQ=255)
                                                   aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGGTTGGCCCCCCGGGGGCCCCGCGGGCCCGCGCCCCCCGCCGTGAAGGGGGTCGCaaagcg                                                                          >  5:340183/1‑130 (MQ=255)
                                                     tctGCTGGCGAGATGGACACTGGCGCCGGGCTACGAACACCTGGCGGCGGACTACGGCAGCCACGACGCGGTGTTAGCCCTCCCACGCCCTCCGCTGCCCCGCGCCCCCCTCACGGCAATCGTACGCATAGAGCgg                                                                       >  8:158447/1‑135 (MQ=255)
                                                     tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCTCCCGCCAGGGCGCACGGCTGGCCCCCGCCCCGCGGAGGGTAGTGGTTCGCGTCGAGCg                                                                        >  2:152890/1‑135 (MQ=255)
                                                     tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCCCGGGGAACGGCAGCCCCCCCGCCCCCGCCCAGAAGTGGTAGGCATcgcgcg                                                                        >  1:304700/1‑135 (MQ=255)
                                                     tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGGTAGGCCGCCCCGGCGCCCCGGGGGCCCCCGCCCCCCGCAGGGAAGTGGTTCGCCTcgcgcg                                                                        >  3:242566/1‑135 (MQ=255)
                                                     tctGAGGGCGAGTTGGACACTGGCGCCGGGGTACCAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCg                                                                        <  3:379307/135‑1 (MQ=255)
                                                          tGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                    <  6:79918/134‑1 (MQ=255)
                                                          tGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCACCAGGGCGCACGGCTGACCCGCGACCCGCGCAGTGAAGTGGTAAGCAACGAGGGTGACgg                                                                  >  7:314983/1‑136 (MQ=255)
                                                          tGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGGGTTCGCCCTCCCGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                    <  4:433422/134‑1 (MQ=255)
                                                            gCGAGGTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACggg                                                                 <  4:48512/135‑1 (MQ=255)
                                                                        cgggCGCCGGGGAACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                     <  5:157084/133‑1 (MQ=255)
                                                                        cTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                     <  4:174880/135‑1 (MQ=255)
                                                                        cTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                     <  3:74291/135‑1 (MQ=255)
                                                                        cTGGCGCCGGGGTACGAACACCGGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                     <  1:195594/135‑1 (MQ=255)
                                                                         gggCGCGGGGGTACTACCACGTGGGGGGGGACTTCGGCCGCCTCGACGGGGTGTGCGCCAGCCAGGGCTAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACg                                                    <  6:293449/134‑1 (MQ=255)
                                                                            cgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGt                                                   <  6:219478/133‑1 (MQ=255)
                                                                            cgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCTCCAGGGGCCACGGCGGACCCCCGCCCCCCGCAGTGAAGTGGTAAGCATAGAGCGTGGCGGGATCAACAACAACCTcc                                                 >  7:138252/1‑134 (MQ=255)
                                                                              ccGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCCGGGCGAACGGCTGACCCGCGCCCCGCGCAGTGAAATGGTACGCATCGAGCGGGCCGGGGTCAACTACTACGTCAAg                                               >  2:39509/1‑135 (MQ=255)
                                                                                   ttACGAACACGTGGCGGGTGAGTGCGGCAGCCTGGAGGCGGGTTTCGCCCTCCTGGGCTAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTTCGTCAAGCGCTa                                          <  1:126689/134‑1 (MQ=255)
                                                                                              tGGCGGCGTACTTCGGCAGCCTCGACGCGGTGTGCGCCCGCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCg                               <  5:262950/135‑1 (MQ=255)
                                                                                                gcggcgGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCCGGGCGCACGGCGGCCCCGCGCCCCGCGCAGTGCAGTGGTTCGCATCGCGCGTGGCGGGGTGCACAACTACGTCCAGCGCGACACCCGGGGCgg                              >  4:174997/1‑134 (MQ=255)
                                                                                                           cGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATgc                   <  7:300087/134‑1 (MQ=255)
                                                                                                            ggCAGCCTCGACGCGGTGTTCGGCCTCCTGGGCGAACGGCTGACCCGCGACCCGCGCAGTGAAGTGGTACGCATCGAGCGTGACGGGATAAACTACTACGTAAAGCGCTACACCGGGGCCGGCAAGCACATgc                   <  3:153994/133‑1 (MQ=255)
                                                                                                                           gtgtTCGCGCTCCTGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTAGGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg  <  3:192664/135‑1 (MQ=255)
                                                                                                                           gtgtTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg  <  6:438704/135‑1 (MQ=255)
                                                                                                                           gtgtGCGCCCGCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg  <  2:304700/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                 
GGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGG  >  ACB122/411211‑411468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑