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

TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111612‑4111883
                                                                                                                                         |                                                                                                                                        
ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CTCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                      >  2:834608/1‑140 (MQ=255)
     ttttATTCAAACCTTCAAGCGATTCAGATAGCCCCAGTTTAATCGGTTCAACAGCGAGGGTCACCCCTTTTT‑CGCCTTGGCCCGCGAAAAAATAC‑CGCGGTGCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCcccc                                                                                                                                 <  2:3318/139‑1 (MQ=255)
      tttATCCAACCCTACAAGCGATTCAGATAGCCCCAGTTTAATCGGTTCAACAGCGAGGGCCCGCCCCTTTT‑CGCCTTTGCCCGCGAAAACATAC‑CGCAGTGCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTTCCccg                                                                                                                                <  2:283447/140‑1 (MQ=255)
              aaCCTTCAGGCGATTCAGATAGCGCCAGCTTAATCGGTCCAACAGCGAGGGTCACCCCTTTTT‑CGCCGTTGCCCGCGACAACATAC‑CGCGGTGCCCTTTCTGTCTGGGTGTAATACCGTTTTTTTTCCCCCGCCGTtagt                                                                                                                        <  2:177910/140‑1 (MQ=255)
                    cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                  >  1:752707/1‑140 (MQ=255)
                                     gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 >  2:256628/1‑140 (MQ=255)
                                     cccAGTTTAATCGGTTCAACAGCGAAGGTCACCCCTTTTC‑CGCCGTTGCCCGCAAAAACAAAC‑CGCAGTGCCCCTTCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 <  2:35271/139‑1 (MQ=255)
                                      ccAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                >  1:1254343/1‑140 (MQ=255)
                                       cactttAATGGTTACAACAGCGAAGTCCACCCCTTTTT‑CGCCGTGGCCCGCAAAAAAATAA‑CGCATTCCCCCTTCTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                                                                               <  2:31299/136‑1 (MQ=255)
                                          cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                            >  1:857979/1‑140 (MQ=255)
                                                   ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                   >  1:648646/1‑140 (MQ=255)
                                                     aaaCAGCGAAGGCCACCCCCTTTT‑CGCCGTTGTCCGCAAAAACATAA‑CGCAGTGCCCTTTCTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                 <  2:1135092/139‑1 (MQ=255)
                                                     aaaCAGCAAGGGCCACCCCTTTTC‑CGCTGTTGCCGGCGA‑AACAAAAGCGCAGCCCCCTTTCTGTCTGGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                 <  2:721234/139‑1 (MQ=255)
                                                            gAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                          >  1:831303/1‑140 (MQ=255)
                                                             aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         >  1:1056686/1‑140 (MQ=255)
                                                             aaGGTCACCCCCTTTC‑CGCCGTTGCCCGCGACAACATAA‑CGCAGTGCCCCTTCTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:1056686/140‑1 (MQ=255)
                                                              agggcccccccTTTTCCGCCGTGGCCCGCAAAAACATAC‑CGCGGTGCCCTTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                         <  2:563116/133‑1 (MQ=255)
                                                               ggTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                       <  2:556214/140‑1 (MQ=255)
                                                                          ttt‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                            >  1:708284/1‑140 (MQ=255)
                                                                          ttt‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                            >  1:1223344/1‑140 (MQ=255)
                                                                               gcggTTGCCCGAGACAACAAAA‑CGCGGTCCCCCTTCTGTCTGGTTATAATATTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                        <  2:859829/137‑1 (MQ=255)
                                                                                  tttGCCCGCAACAAAATAC‑CGCGGTGCCCTTTTTGTCTGGGTGTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                     <  2:204161/139‑1 (MQ=255)
                                                                                           gACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTAtt                                            >  1:608534/1‑140 (MQ=255)
                                                                                                   aa‑CGCAGTGCCCCTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgct                                    <  2:752707/140‑1 (MQ=255)
                                                                                                         aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                               >  1:715603/1‑140 (MQ=255)
                                                                                                         aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                               >  2:241740/1‑140 (MQ=255)
                                                                                                           tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  1:31073/1‑140 (MQ=255)
                                                                                                            gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                            >  2:408194/1‑140 (MQ=255)
                                                                                                            gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                            >  2:1057351/1‑140 (MQ=255)
                                                                                                                cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                        >  2:666379/1‑140 (MQ=255)
                                                                                                                  tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                      >  2:902747/1‑140 (MQ=255)
                                                                                                                     gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                   >  1:101841/1‑140 (MQ=255)
                                                                                                                     gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                   >  1:965571/1‑140 (MQ=255)
                                                                                                                     gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                   >  1:1009710/1‑140 (MQ=255)
                                                                                                                          ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  2:1043399/1‑140 (MQ=255)
                                                                                                                              taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg          >  1:594678/1‑140 (MQ=255)
                                                                                                                                atatttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc        <  1:834608/134‑1 (MQ=255)
                                                                                                                                ataaCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc        <  2:31073/140‑1 (MQ=255)
                                                                                                                                      ttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc         <  1:129469/133‑1 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc         >  2:129469/1‑133 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:681562/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:265977/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:176067/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:1131365/1‑140 (MQ=255)
                                                                                                                                      tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  2:91129/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAA‑CGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111612‑4111883

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