Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT35.7% 47.3 / 16.7 28N156K (AAC→AAA) yiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (15/3);  new base T (0/10);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 2.18e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111617‑4111885
                                                                                                                                   |                                                                                                                                            
ttttATTAAACCCTTAAAGCGTTTCAGATAGCCCCAGTTTAATCGGTTCAACAGCGAAGGTCACCCCT‑TTTTCGCCGTTGCCCGCGAAAAAAAAACGCGGTGCCCCTTCTGTCTGGGTGAAATAACGTTTTTTTTCcccc                                                                                                                                     <  2:30214/139‑1 (MQ=255)
               cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                      >  2:241411/1‑140 (MQ=255)
               cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                      >  1:334691/1‑140 (MQ=255)
               cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                      >  1:364954/1‑140 (MQ=255)
                aaGCAATTCAGATAGCGCCAGTTTATTCGTTTCAACAGCGAAGGCCACCCCC‑TTTCCGCGGTTGCCCGCAACAACATACCGCGGTGCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACt                                                                                                                     <  2:364954/140‑1 (MQ=255)
                      ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                               >  1:616362/1‑140 (MQ=255)
                                gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                     >  1:772614/1‑140 (MQ=255)
                                gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                     >  1:509016/1‑140 (MQ=255)
                                cccAGCTTAATGGGTTCAACAGGGAAGGTCCCCCCT‑TTTTCGCCGTTGTCCGCAACAACATACCGCAGTGCCCCTTCTGTCTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                     <  2:212276/139‑1 (MQ=255)
                                     cTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                                >  1:212276/1‑140 (MQ=255)
                                       tAATCGTTCAAACAGCGAAGGTCACCCCT‑TTTCCGCCTTTGCCCGCAAAAACATACCGCGGTGCCCCTTCTGTCTGGGTTTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                              <  2:738483/140‑1 (MQ=255)
                                        aaccGGTCCAACAGCGAAGGTCCCCCCC‑TTTTCGCGGTTGCCCGCAACAACATACCGCAGTGCCCTTTCTGTCTGGGTGTAATATGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                                                                             <  2:113646/137‑1 (MQ=255)
                                            gtttCAACAG‑GAAGGGCCCCCCCTTTTTCGCCTTTGCCCGCAAAAACATAACGCAGTGCCCTTTCTGTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                         <  2:584862/138‑1 (MQ=255)
                                            gttacAACAGCGAAGGCCCGCCCT‑TTTTCGCCGTTGTCCGCAACAACATACCGCAGTGCCCCTTTTTTCGGGGTGTATTATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                         <  2:719975/136‑1 (MQ=255)
                                            gatccAACGGCGAGGGACCCCCCT‑TTTCCCCCTTTGTCCCCAAAAACAAAACGCAGTGCACCTTCTGCCTTGGTGTAATACTGTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                         <  2:772614/136‑1 (MQ=255)
                                              tacAACAGCGAAGGTCCCCCCT‑TTTCCGCCTTTGCCCGCAAAAACAAAACGCGGTCCCCCTTCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                       <  2:777620/138‑1 (MQ=255)
                                               tCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                      >  1:795974/1‑140 (MQ=255)
                                                  aCAGCGAGGGTCCCCCCT‑TTTCCCCCTTTGTCCGCAAAAACATAACGCAGTGCCCCTTCTTTCTGGGTGTAATACTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                   <  2:132661/140‑1 (MQ=255)
                                                   cAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                                                                                  >  1:765959/1‑140 (MQ=255)
                                                        aaGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                             >  1:777620/1‑140 (MQ=255)
                                                          ggccgcccct‑TTTTCGCCGTTGTCCGCAACAACATACCGCAGTGCCCCTTTTGTCTGGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                           <  2:509016/134‑1 (MQ=255)
                                                                            ggTTGCCCGCAACAACATAACGCGGTGCTCCTTCTGTCTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCt                                                        <  2:596325/139‑1 (MQ=255)
                                                                            ctttGTCCGCAACAACATAACGCAGTCCCCTTTCTGTCTCGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGTGCGGATCACCg                                                          <  1:210227/138‑1 (MQ=255)
                                                                             gtggcccGCAAAAACATAACGCGGTCCCCCTTCTGTCTGGGTTAAATAACTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTg                                                       <  2:795974/135‑1 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:490092/1‑140 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  2:232442/1‑140 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCCGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:79966/1‑140 (MQ=255)
                                                                                                        cctttCTGTCTGGGTGTAATAATTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                            <  2:618981/137‑1 (MQ=255)
                                                                                                          cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTCGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          >  1:624193/1‑140 (MQ=255)
                                                                                                                      tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:713442/140‑1 (MQ=255)
                                                                                                                      tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:232442/140‑1 (MQ=255)
                                                                                                                      tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              >  2:713442/1‑140 (MQ=255)
                                                                                                                              ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  1:623214/138‑1 (MQ=255)
                                                                                                                              tgtttttttCTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTCGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  2:624193/139‑1 (MQ=255)
                                                                                                                               gtttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTgg     >  1:428335/1‑140 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGTGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:234554/1‑140 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:590421/1‑140 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:326605/1‑140 (MQ=255)
                                                                                                                                  tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAtc  >  1:113047/1‑139 (MQ=255)
                                                                                                                                  tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:1700/1‑140 (MQ=255)
                                                                                                                                  tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:523871/1‑140 (MQ=255)
                                                                                                                                  tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  1:41852/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                            
ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCC‑TTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGA‑‑GTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111617‑4111885

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