Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,497,0240TG20.6% 60.8 / 16.4 34intergenic (+32/+33)ypdK/alaCinner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base G (3/4);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCA  >  NC_000913/2496889‑2497169
                                                                                                                                       |                                                                                                                                                 
ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGttt                                                                                                                                      <  2:119690/149‑1 (MQ=255)
    tAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttttca                                                                                                                                  <  3:119573/149‑5 (MQ=255)
          tAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttttc                                                                                                                                   >  8:82219/1‑139 (MQ=255)
          tAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttttc                                                                                                                                   <  7:82219/142‑4 (MQ=255)
          tAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTga                                                                                                                            >  3:78256/1‑149 (MQ=255)
                   ggAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtgt                                                                                                                   <  2:159694/149‑1 (MQ=255)
                         aTGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCa                                                                                                                                                 >  5:95675/1‑113 (MQ=255)
                         aTGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCa                                                                                                                                                 <  6:95675/113‑1 (MQ=255)
                                   aaaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCg                                                                                                   >  1:142200/1‑149 (MQ=255)
                                    aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGttt                                                                                                                                      >  4:55628/1‑113 (MQ=255)
                                    aaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGttt                                                                                                                                      <  3:55628/113‑1 (MQ=255)
                                         ttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTc                                                                                                             >  6:53193/1‑133 (MQ=255)
                                         ttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTc                                                                                                             <  5:53193/133‑1 (MQ=255)
                                                   gCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgctgc                                                                                   >  7:52000/1‑149 (MQ=255)
                                                      tttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCAt                                                                                                                                                   >  3:121014/1‑82 (MQ=255)
                                                      tttCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCAt                                                                                                                                                   <  4:121014/82‑1 (MQ=255)
                                                               tGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgtagt                                                                        >  3:49212/1‑148 (MQ=255)
                                                                       gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTc                                                                <  4:49212/148‑1 (MQ=255)
                                                                       gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTc                                                                <  8:88226/148‑1 (MQ=255)
                                                                       gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCg                                                               <  2:142200/149‑1 (MQ=255)
                                                                                             tgaggaTCTAAAAGCAATCTGTGCAATGAAAAACAGGAGCCATCACGCTGCTGTTTGCGTTTCTGCCGACGGGGTTTATGCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGCGTAGTAAGCCGTCGGCGCGGCGCCTCGCTTTAAGcc                                         >  4:7329/1‑149 (MQ=255)
                                                                                                tgatCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCctc                                      <  8:52000/149‑1 (MQ=255)
                                                                                                   tCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCgg                                                                             <  7:49760/107‑1 (MQ=255)
                                                                                                   tCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCgg                                                                             >  8:49760/1‑107 (MQ=255)
                                                                                                       aaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTTAATATGTTtgctg                                                                                    <  4:146645/96‑1 (MQ=255)
                                                                                                            aaaCTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCCGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGAc                          <  4:78256/149‑1 (MQ=255)
                                                                                                                      aTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTc                >  4:138751/1‑149 (MQ=255)
                                                                                                                       tGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGATGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTcg               >  1:58661/1‑149 (MQ=255)
                                                                                                                         gcaaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgta                                                                          <  7:112588/86‑1 (MQ=38)
                                                                                                                         gcaaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgta                                                                          >  8:112588/3‑88 (MQ=38)
                                                                                                                            aaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtt          <  3:138751/149‑1 (MQ=255)
                                                                                                                            aaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtt          >  1:41036/1‑149 (MQ=255)
                                                                                                                             aCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                 <  2:147485/109‑1 (MQ=255)
                                                                                                                             aCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGc                                                 >  1:147485/1‑109 (MQ=255)
                                                                                                                              cAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtttt        <  2:41036/149‑1 (MQ=255)
                                                                                                                                    ccATCAGGCTCCTGTTTGCATTTCTGACGACGATGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCa  <  2:58661/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
GGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCA  >  NC_000913/2496889‑2497169

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