Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT47.8% 21.4 / 19.7 23N156K (AAC→AAA) yiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (11/1);  new base T (0/11);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 9.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111622‑4111885
                                                                                                                              |                                                                                                                                          
ttCAACCTTTCAAGCGATTCAGATAGCGCCAGTTTAATGGTTTCAACAGCGAAGGTCACCCCCTTTC‑CGCCTTTGCCCGCAAAAACATAACGCGGTCCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGt                                                                                                                              <  2:630928/140‑1 (MQ=255)
    acccTTAAAGCGATTAAGATAGCCCCAGCTTAATCGGTTCAACAGGGAAGGCCCCCCCTTTTC‑CGCCGTTGCCCGCAAAAACATACCGCAGTCCCCCTTCTGTCTCGGTGTAATACCTTTTTTTTTCCCCCGCCGTtagt                                                                                                                          <  2:924614/138‑1 (MQ=255)
       cTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGc                                                                                                                       >  2:835367/1‑140 (MQ=255)
                           cccAGCTTAATCGTTAAAACAGCGAGGGTCACCCCTTTTC‑CGCCTTTGTCCGCAAAAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   <  2:17047/139‑1 (MQ=255)
                                cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                              >  1:445932/1‑140 (MQ=255)
                                cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                              >  1:685044/1‑140 (MQ=255)
                                cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGCTAGTAGCGACTGCAGTTTTTTCTGGTTTtgtg                                                                                              >  1:62915/1‑140 (MQ=255)
                                  aaaTCGGTACAACGGGGAAGGTCACCCCTTTTT‑CGCCGTTGCCCCCGACAACAAACCGCAGTGCCCCTTCTGTCTGGGTGTAATACCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                            <  2:518365/139‑1 (MQ=255)
                                       gtttCAACGGCAAGGGCCACCCCTTTTT‑CGCCGTTGCCCGCAAAAACATACCGCAGTGCCCCTTCTGTCTGGGTATAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:340426/138‑1 (MQ=255)
                                        gTTAAACAGCAAAGGTCACCCCTTTTT‑CGCCGTTGCCCGCGAAAACATAACGCAGTGCCCTTTCTGTCTGGGTGAAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTc                                                                                      <  2:112303/140‑1 (MQ=255)
                                                   aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggggc                                                                           >  1:340426/1‑139 (MQ=255)
                                                    agggcccccccTTTTCCCCCGTTGTCCGCAAAAAAATACCGCGGTCCCCCTTCTTTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:746041/133‑1 (MQ=255)
                                                       ccacccctTTTT‑CGCCTTTGCCCGCAACAACATACCGCGGTGCCCTTTTTGTCTGGGTTTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt                                                                       <  2:392020/136‑1 (MQ=255)
                                                           ccccTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                                                                   >  1:133040/1‑140 (MQ=255)
                                                                        tttGCCCGCAACAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  1:472624/139‑1 (MQ=255)
                                                                                             cAGTGCACCTTCTGTCTCGGTGTAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcg                                  <  1:591981/140‑1 (MQ=255)
                                                                                                tGCACCTTCTGTCTCGGAGGAATAACGTTTGCTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGATTTGGTCATTAATAGAGTCAGCAAGCGGATCACCGCTGCCATGTATTCCTGAgcttgcgctt                               >  2:391709/1‑140 (MQ=255)
                                                                                                 gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  2:498689/1‑140 (MQ=255)
                                                                                                                   taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg            >  1:339872/1‑140 (MQ=255)
                                                                                                                   taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg            >  2:265775/1‑140 (MQ=255)
                                                                                                                    aatatCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc           <  1:712757/140‑1 (MQ=255)
                                                                                                                     ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:148646/1‑140 (MQ=255)
                                                                                                                     ataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:631997/1‑140 (MQ=255)
                                                                                                                         cGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg      <  1:265775/140‑1 (MQ=255)
                                                                                                                           tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:964542/1‑140 (MQ=255)
                                                                                                                           tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    <  1:784919/140‑1 (MQ=255)
                                                                                                                           tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:859165/1‑140 (MQ=255)
                                                                                                                           tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTGTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:784919/1‑140 (MQ=255)
                                                                                                                             tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  1:201073/1‑140 (MQ=255)
                                                                                                                              |                                                                                                                                          
TTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111622‑4111885

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: ‑