Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT34.6% 48.2 / 15.6 26N156K (AAC→AAA) yiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (16/1);  new base T (0/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGG  >  NC_000913/4111617‑4111886
                                                                                                                                  |                                                                                                                                           
ttttATTAAACCCTTCAAGCTATTCAGATAGCCCCAGTTTAATCGGTTAAACAGCGAAGGTCACCCCTTTTTCGCCTTTGCCCGCGAAAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCcccc                                                                                                                                    <  2:304691/139‑1 (MQ=255)
               cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                     >  1:907153/1‑140 (MQ=255)
               cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                     >  2:447182/1‑140 (MQ=255)
                  gggATTCAGATGGCCCCAGTTTAAGCGTTTCAACAGCGAAGGCCGTCCTTTTTTCGCTGTGGGCCGCCACAAAATACCGCATTCCCCCTTCTTTCTGGGTGAAATAAGTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGc                                                                                                                  <  2:637358/138‑1 (MQ=255)
                    gATTCAGATACCCCCAGCTTAATCGTTACAACAGCGAGGGTCCCCCCTTTTCCGCCGTTGTCCGCGAAAACAAACCGCAGTGCCCCTTCTGTCTGGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAg                                                                                                                <  2:30257/140‑1 (MQ=255)
                      ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                              >  1:418650/1‑140 (MQ=255)
                                gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                    >  2:607579/1‑140 (MQ=255)
                                            ggtcaaaCAGCGAAGGCCCCCCCTTTTTCGCGGTTGCCCGCGACAACAAAACGCAGTCCCCCTTCTGTCTGGGTGTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                        <  2:379315/135‑1 (MQ=255)
                                            ggTACAACAGCAAGGGTCCCCCCTTTTCCGCTGTTGTCCGCAACAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAAGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                        <  2:534909/140‑1 (MQ=255)
                                              ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                      >  1:525729/1‑140 (MQ=255)
                                                  aCAGCAAAGGCCCCCCCTTTTTCCCCTTTGCCCGCAAAAACAAACCGCAGTCCCCCTTCTGCCTCGTTTAAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTCTGTGTTTTGGTCATTaa                                                                                  <  2:418650/140‑1 (MQ=255)
                                                      cGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                                                                              >  1:699675/1‑140 (MQ=255)
                                                        agggccAGCCCTTTTGCGCCGTTGTCCGCGAAAACATACCGCGGTGCCCTTTCTGTCTGGGTGTAATATCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            <  2:525729/135‑1 (MQ=255)
                                                        aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            >  1:725010/1‑140 (MQ=255)
                                                        aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            >  2:199382/1‑140 (MQ=255)
                                                                            tttGCCCGCAAAAACATAACGCAGTGCCCCTTCTGTCTGGTTGTAATAAGTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                        <  1:665549/139‑1 (MQ=255)
                                                                              tgcccGCGACAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAATGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                      <  1:447182/137‑1 (MQ=255)
                                                                              tgcccGCCACAACATACCCCAGTGCACTTTTTGTCTTGGTGTAATAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                      <  1:893574/137‑1 (MQ=255)
                                                                                                 cAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                   >  2:394688/1‑140 (MQ=255)
                                                                                                    tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                >  2:357140/1‑140 (MQ=255)
                                                                                                    tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                >  1:761696/1‑140 (MQ=255)
                                                                                                         cTTTTGTCTGGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                           <  2:725010/140‑1 (MQ=255)
                                                                                                         cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                           >  2:765008/1‑140 (MQ=255)
                                                                                                           tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                         >  2:591308/1‑140 (MQ=255)
                                                                                                              gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                      >  1:479203/1‑140 (MQ=255)
                                                                                                                ctcGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                    >  1:500992/1‑140 (MQ=255)
                                                                                                                     tgtAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct               <  2:698901/140‑1 (MQ=255)
                                                                                                                     tgtAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct               <  1:394688/140‑1 (MQ=255)
                                                                                                                         ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc           >  1:245712/1‑140 (MQ=255)
                                                                                                                              tttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTgg      <  1:199382/139‑1 (MQ=255)
                                                                                                                               ttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGTTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg       <  2:25232/138‑1 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGTTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg       >  1:25232/1‑138 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  1:88855/1‑140 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  1:398746/1‑140 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  1:158490/1‑140 (MQ=255)
                                                                                                                               tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGGGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:28658/1‑140 (MQ=255)
                                                                                                                                 tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg   >  2:869268/1‑140 (MQ=255)
                                                                                                                                  gtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg  >  1:239637/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                           
ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGG  >  NC_000913/4111617‑4111886

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