Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,3970AC40.5% 42.1 / 15.5 42L198F (TTA→TTCyiaN2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (17/8);  new base C (0/17);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 4.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCG  >  NC_000913/3745260‑3745530
                                                                                                                                         |                                                                                                                                     
tCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTTGCATGATGATTGGGGCAACGCTGATGCTTTCCTGGTGGGGGCGGGCCAGGCGCCTTAATTTTcc                                                                                                                                     >  2:392895/1‑140 (MQ=255)
         tttATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGGGGCAGGCCAGTCGCCTTAATTTACCCCCCCAGCa                                                                                                                            >  2:658810/1‑140 (MQ=255)
         tttATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCa                                                                                                                            <  1:979227/140‑1 (MQ=255)
                 cTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTTGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCaa                                                                                                                    <  1:648850/140‑1 (MQ=255)
                 cTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCaa                                                                                                                    <  2:268068/140‑1 (MQ=255)
                 cTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCaa                                                                                                                    <  2:678879/140‑1 (MQ=255)
                          ttCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGGCAGGCCCGTCGCCTTAATTTACCCCCCCCGCAAAAAAGAAAGAAGCAgg                                                                                                           >  2:91946/1‑140 (MQ=255)
                          ttCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGGGCAACGCTGATGCTTACCCGGGGGGGGGAGGCCAGTCGCCTTTATTTTCCCCCCCAGCAAAAAACAAAGATGCAgg                                                                                                           >  2:1000245/1‑140 (MQ=255)
                           tCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGGGCAACGCTGATGCTTACCTTGTGGTGGCAGGCCCGTCGCCTTAATTTACCCCCCCAGCAAAAAACAACGATGCAGGa                                                                                                          >  2:648850/1‑140 (MQ=255)
                                 ggATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTg                                                                                                    <  1:115191/140‑1 (MQ=255)
                                  gATTATCCAGCAGTAAGCTGTTTAGGGGGGGCATTCCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACTTGGTGGTGGCAGCCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTgg                                                                                                   <  2:757536/140‑1 (MQ=255)
                                   aTTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGc                                                                                                  <  1:392895/140‑1 (MQ=255)
                                   aTTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGc                                                                                                  <  1:952294/140‑1 (MQ=255)
                                     tATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTTCCTGGGGGTGGCAGGCCCGTCGCCTTAATTTACCCCGCCCGCCAAAAGCAACGATGCAGGAAATCTGGCAc                                                                                                >  1:644800/1‑140 (MQ=255)
                                           tCAGTAAGCTGTTTATGGCGGGCGTTGCCCCTGGCATGATGATGGGCGCAACGCTGATGTTTACCTGGTGGGGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt                                                                                          <  2:824350/140‑1 (MQ=255)
                                            cAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTg                                                                                         <  1:658810/140‑1 (MQ=255)
                                             aGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGTAGAAGGGCGGAACGCTGATGCTTTCCTGGGGGTGGCAGGGCCGGCGCCTTTATTTTCCCCCCCAGCAAAAAACAACGAAGCAGGAAAACTGGGACTCtttttt                                                                                        >  2:952294/1‑138 (MQ=255)
                                               tAAGCTGTTTATGGCGGGCATTGCCCCTGGCGTGATGATTGGCGCCACGCTGATGCTTTCCTGGGGGGGGCCGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCACGGAATCTGGGACTCTTTTTtct                                                                                      >  2:4934/1‑140 (MQ=255)
                                                gacaggTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtctc                                                                                     <  1:834696/135‑1 (MQ=255)
                                                 aGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt                                                                                          <  2:903515/134‑1 (MQ=255)
                                                 aGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGGGGCAGGCCAGGCGCCTTTATTTACCCCCCCCGCAAAAAACAAAGATGCAGGAAAACTGGCACCCTTTTTTCTcc                                                                                    >  2:749721/1‑140 (MQ=255)
                                                 aGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt                                                                                          >  1:903515/1‑134 (MQ=255)
                                                    tGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGt                                                                                        <  2:878458/133‑1 (MQ=255)
                                                    tGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGt                                                                                        >  1:878458/1‑133 (MQ=255)
                                                               ggCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCt                                                                      <  1:595469/140‑1 (MQ=255)
                                                                gCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCCGCAAAAAACAACGATGCAGGAAAACTGGCACTCTTTTTTCTCCCGTATCTGGgctgtc                                                                     >  2:834696/1‑136 (MQ=255)
                                                                         cTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCa                                                            <  2:535787/140‑1 (MQ=255)
                                                                            gCATGTTGTGGGGCGCAACGTTGGTGCTTCCGTGGTGGTGGCAGGCGAGTCGCCTTAATTTTCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGtc                                                         <  2:661580/140‑1 (MQ=255)
                                                                                 atgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTcatcat                                                    >  1:54990/1‑140 (MQ=255)
                                                                                     tgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCattatt                                                >  1:632258/1‑140 (MQ=255)
                                                                                       ggCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTgg                                              >  1:495281/1‑140 (MQ=255)
                                                                                        gcgcAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGc                                             <  1:4934/140‑1 (MQ=255)
                                                                                        gcgcAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGc                                             >  1:583364/1‑140 (MQ=255)
                                                                                              cGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGttt                                       >  2:138105/1‑140 (MQ=255)
                                                                                              cGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGttt                                       >  2:652637/1‑140 (MQ=255)
                                                                                                   aTGCTTACCTGGTGGTGGCAGGCCAGTCGCTTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCGCTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTgg                                              <  2:924636/128‑1 (MQ=255)
                                                                                                   aTGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCGCTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTgg                                              >  1:924636/1‑128 (MQ=255)
                                                                                                      cTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg                               <  2:54990/140‑1 (MQ=255)
                                                                                                       ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGc                              <  1:138105/140‑1 (MQ=255)
                                                                                                       ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGc                              <  1:652637/140‑1 (MQ=255)
                                                                                                        tACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGcc                             >  1:660241/1‑140 (MQ=255)
                                                                                                          ccTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTccctcctt                           >  2:39662/1‑134 (MQ=255)
                                                                                                          ccTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCTCCCGccac                           >  2:382321/1‑138 (MQ=255)
                                                                                                            tggtggtggCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt                         <  2:660241/140‑1 (MQ=255)
                                                                                                                   ggCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg                               <  1:805504/127‑1 (MQ=255)
                                                                                                                   ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg                               >  2:805504/1‑127 (MQ=255)
                                                                                                                   ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt                         <  1:515329/133‑1 (MQ=255)
                                                                                                                   ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt                         >  2:515329/1‑133 (MQ=255)
                                                                                                                    gCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGGCCCGTTTTTCCCtc                 >  2:295771/1‑139 (MQ=255)
                                                                                                                     cAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTTCTCCGGGCTGTTTTCCCCCAc                >  1:285781/1‑140 (MQ=255)
                                                                                                                                 cTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGc    <  1:295771/140‑1 (MQ=255)
                                                                                                                                   tAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCg  <  2:495281/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCG  >  NC_000913/3745260‑3745530

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