Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7520TC37.1% 41.6 / 11.6 35K155E (AAA→GAA) yiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base T (18/4);  new base C (0/13);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111617‑4111885
                                                                                                                                        |                                                                                                                                     
ttttATTAAACCCTTAAAGCGATTCAGATAGCGCCAGTTTAATGGGTTCAACAGCGAAGGCCACCCCCTTTC‑CGCCGTTGTCCGCAACAACATAACGCGGTGCCCCTTCTGTCTCGGTTTAATAACTTTTTTTTTCcccc                                                                                                                                   <  2:162434/139‑1 (MQ=255)
                      ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  2:262506/1‑140 (MQ=255)
                       tCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttt                                                                                                            >  1:502048/1‑140 (MQ=255)
                               cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                    >  1:245798/1‑140 (MQ=255)
                                   aGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGCTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                                                                                >  2:145690/1‑140 (MQ=255)
                                                  aCAGCGAGGGTCCCCCCTTTTC‑CGCCGTTGTCCGCGAAAACAAACCGCGGTCCCCTTTCTGTCTGGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                 <  2:265044/140‑1 (MQ=255)
                                                  aCAGCAAAGGCCACCCCCTTTT‑CGCCGTTGCCCGCGAAAACATACCGCAGTGCCCTTTCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa                                                                                 <  2:265045/140‑1 (MQ=255)
                                                     gCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTCCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgg                                                                              >  1:252497/1‑140 (MQ=255)
                                                        aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCGGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  1:193809/140‑1 (MQ=255)
                                                         gagacggCCCCTTTTCCGCGGTTGCCCGCAACAACATAACGCAGTCCCCTTTCTGTCTGGGTGAAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:855379/134‑1 (MQ=255)
                                                           ggcccccccTTTTCCGCTTTTGCCCGCAAAAACATACCGCAGTGCCCCTTTTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                         <  2:234107/135‑1 (MQ=255)
                                                           gTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg                                                                        >  1:297145/1‑140 (MQ=255)
                                                                     ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                              >  1:717782/1‑140 (MQ=255)
                                                                      tt‑CGCCGTTGCCCGCAACAACATAACGCGGTCCCCCTTCTGTCTGGGTGTAATAATTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCa                                                             <  2:245798/140‑1 (MQ=255)
                                                                         cGCTGTTGCCCGCGAAAACATACCCCAGTGCACCTTCTGTCTCGTTTTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAcc                                                           <  2:297145/140‑1 (MQ=255)
                                                                          cccGTTGCCCGCGAAAACATAACGCAGTCCCCTTTCTGTCTGGGTTTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                          <  1:25437/139‑1 (MQ=255)
                                                                           ccGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGactgtctc                                                         >  1:855379/1‑133 (MQ=255)
                                                                             tttGTCCGCAACAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  1:262506/139‑1 (MQ=255)
                                                                                  ccGCGACAACATACCGCGGTCCCCCTTCTGTCTGGTTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTCCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                  <  2:252497/140‑1 (MQ=255)
                                                                                           caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                         >  2:446027/1‑140 (MQ=255)
                                                                                                  cAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCAGCCATGTATTCCTGAgctttcg                                  >  2:699851/1‑140 (MQ=255)
                                                                                                     tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:277204/1‑140 (MQ=255)
                                                                                                      gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                              >  2:408140/1‑140 (MQ=255)
                                                                                                      gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctct                              >  1:616655/1‑138 (MQ=255)
                                                                                                          tttCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          <  2:538227/139‑1 (MQ=255)
                                                                                                          tttCTGTCTGGGTGAAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          <  2:704184/139‑1 (MQ=255)
                                                                                                          cTTCTGTCGCGTTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                          <  2:277204/140‑1 (MQ=255)
                                                                                                               ttttcGGTTAAAAAACGTTTTTTTTCCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                     <  1:145690/137‑1 (MQ=255)
                                                                                                                 ctcGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttttt                   >  2:560644/1‑140 (MQ=255)
                                                                                                                    ggTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                <  2:503299/140‑1 (MQ=255)
                                                                                                                      tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              >  1:503299/1‑140 (MQ=255)
                                                                                                                       ttaatacCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTg             <  2:717782/139‑1 (MQ=255)
                                                                                                                         aatactgtttttttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc           <  1:855686/134‑1 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          <  1:868780/134‑1 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc          >  2:868780/1‑134 (MQ=255)
                                                                                                                                tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  2:573943/1‑140 (MQ=255)
                                                                                                                                  tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:724380/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                     
ATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111617‑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: ‑