Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,404,2150TG28.1% 51.2 / 10.3 32G244G (GGT→GGG) yhdHacrylyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (2/21);  new base G (9/0);  total (11/21)
Fisher's exact test for biased strand distribution p-value = 1.96e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTC  >  NC_000913/3404084‑3404353
                                                                                                                                   |                                                                                                                                           
tttGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTAc                                                                                                                                     >  1:271899/1‑140 (MQ=255)
    ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTg                                                                                                                                 <  2:451917/140‑1 (MQ=255)
        aTCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGccca                                                                                                                             >  2:191359/1‑138 (MQ=255)
          cccGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAAcc                                                                                                                           <  1:191359/140‑1 (MQ=255)
              tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTg                                                                                                                                 <  2:734778/130‑1 (MQ=255)
              tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTg                                                                                                                                 >  1:734778/1‑130 (MQ=255)
              tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATCTGCCCACCACgg                                                                                                                       >  2:608414/1‑140 (MQ=255)
              tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACgg                                                                                                                       >  1:773897/1‑140 (MQ=255)
               ccTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGt                                                                                                                      <  2:8827/140‑1 (MQ=255)
                  ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                   <  2:739963/140‑1 (MQ=255)
                  ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                   <  2:238494/140‑1 (MQ=255)
                  ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                   <  1:702021/140‑1 (MQ=255)
                   tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCCACCACGGTCATg                                                                                                                  >  1:381175/1‑140 (MQ=255)
                                   ggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTCCAACCACCGTCATGCCCTTTATTCTGCGTa                                                                                                  >  2:770739/1‑140 (MQ=255)
                                    gctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTaa                                                                                                 <  1:191006/140‑1 (MQ=255)
                                        gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTTCTCTTCCCACCACGGGCCTGCCATTTTTTTTGCGTAAttgc                                                                                             >  2:387043/1‑137 (MQ=255)
                                         ggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCCTTTTTTTTTCGTAATATcc                                                                                            >  2:88951/1‑140 (MQ=255)
                                             aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCCCTTCCCACCACGGGCATGCCCTTTTTTTTGCGTGATGTTCGttt                                                                                        >  2:203327/1‑140 (MQ=255)
                                             aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATCTGCCCACCCCGGTCATGCCATTTTTTTTGCGTAATGTCCGttt                                                                                        >  2:181550/1‑140 (MQ=255)
                                              ttGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCCCGGGCATGCCATTTTTTTTGCGTAATGTCCGtttt                                                                                       >  2:512283/1‑139 (MQ=255)
                                               gacagccGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑c                                                                                     <  1:642937/135‑1 (MQ=255)
                                                gACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑Ca                                                                                    <  2:584625/140‑1 (MQ=255)
                                                gACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑Ca                                                                                    <  1:598889/140‑1 (MQ=255)
                                                 acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑Caa                                                                                   <  2:647654/140‑1 (MQ=255)
                                                  cacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTTCTCTGCCCACCACGGTCATGCCATTTATTTTGGGTAAAGTCCGTTTG‑CAAg                                                                                  >  1:163306/1‑140 (MQ=255)
                                                    ccGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCTGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAggg                                                                                <  1:578229/140‑1 (MQ=255)
                                                    ccGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAggg                                                                                <  1:512283/140‑1 (MQ=255)
                                                      gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGt                                                                              <  1:239302/140‑1 (MQ=255)
                                                        tGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTTCTCTGCCCACCCCCGGGCTGCCATTTATTCTGCGTAATGTCCGTTTGCCAAgggggg                                                                             >  2:239902/1‑138 (MQ=255)
                                                               aaaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTATTCCGGCCACCCCGGTCCCGCCCATTTTTCTGCGGGATGTTCGGTTG‑CCCGGGGTgggtctggt                                                                     >  2:642937/1‑133 (MQ=255)
                                                                 aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCCAACACGGGTATGGCCTTTTTTTTTCGTAAATTTCCTTTT‑CAAAGGGGGGGTTTAGTaa                                                                   >  2:572422/1‑140 (MQ=255)
                                                                 aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCATTTATTCTGCGGAATGTCCGTTTG‑CAAGGGGGGGATTCAGTaa                                                                   >  1:306322/1‑140 (MQ=255)
                                                                             aGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCCACCACGGGCATGCCATTTTTTCTGCGGAATGTCCGGTTG‑CAAAGGGGGGGTTCAGTAATGACGCCACCAg                                                       >  2:84629/1‑140 (MQ=255)
                                                                                  tGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCCACCCCGGCCATAGCACTTTTTCTGCGGAATTTACCTTTT‑CAAGGGGTGGGTTCAGTAATGACGCCACCAGAAccc                                                  >  2:697000/1‑138 (MQ=255)
                                                                                   ggCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTGTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCGCCAGACcgca                                                 <  2:192696/140‑2 (MQ=255)
                                                                                                    gcggcTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAg                                >  1:93243/1‑140 (MQ=255)
                                                                                                                   ccTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTGTTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAt                 <  1:181550/140‑1 (MQ=255)
                                                                                                                   ccTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGAt                 <  2:474358/140‑1 (MQ=255)
                                                                                                                       tgGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTAc             <  1:496827/140‑1 (MQ=255)
                                                                                                                         gTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCg           <  2:381175/140‑1 (MQ=255)
                                                                                                                            tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa        <  2:773897/140‑1 (MQ=255)
                                                                                                                             ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt       <  2:196771/140‑1 (MQ=255)
                                                                                                                                 ggtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCAtt   <  2:306322/140‑1 (MQ=255)
                                                                                                                                  gtggtTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTc  >  2:389739/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                           
TTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTG‑CAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTC  >  NC_000913/3404084‑3404353

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