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

TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  NC_000913/4111612‑4111883
                                                                                                                                        |                                                                                                                                         
ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                       >  1:260287/1‑140 (MQ=255)
                                     gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                 >  1:193417/1‑140 (MQ=255)
                                      ccAGTTTAATCGTTTCAACAGCGAAGGTCCCCCCCTTTT‑CGCCTTTCCCCGCAACAACATACCGCATTCCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                <  2:50361/140‑1 (MQ=255)
                                          cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                            >  1:22402/1‑140 (MQ=255)
                                             aaccGGTTCAACAGCGAGGGTCACCCCTTTTT‑CGCCGTTGCCCGCAACAACAAAACGCATTCCCCTTTCTCCCTGGGTTTAATAACTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                                                                         <  2:260287/137‑1 (MQ=255)
                                                 ggTCCAACAGCAAGGGCCACCCCTTTTC‑CGCCGTTGTCCGCAACAACAAACCGCGGTCCCCCTTCTGTCGCGGTATAATAGTGTTTTTTTTC‑CACCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                     <  2:366075/140‑1 (MQ=255)
                                                     aaaCAGCGAAGGCCCCCCCTTTTC‑CGCTTTTGCCCGCGAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGAAATACTTTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                 <  2:924258/139‑1 (MQ=255)
                                                                  ccaccccTTTTCCGCCTTTCCCCGCGACAACATACCGCGGTCCCCTTTCTGTCTGGGTATAATACTTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt                                                                     <  2:404613/137‑1 (MQ=255)
                                                                    gCCCCTTTTCCGCGGTTGTCCGCGACAACATACCGCAGTGCCCTTTCTGTCTGGGTTTAATATTTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtgt                                                                   <  2:951555/140‑1 (MQ=255)
                                                                          ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                            >  1:201062/1‑140 (MQ=255)
                                                                          ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                            >  1:404613/1‑140 (MQ=255)
                                                                             ccGCCGTTGTCCGCGAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAc                                                          <  1:1024924/139‑1 (MQ=255)
                                                                               gCCTTTGCCCGAAACAACAGACCGCATTCCCCCTTTTGTCTC‑GTTTAAAGACATTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                        <  2:587575/140‑1 (MQ=255)
                                                                                ccGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGc                                                       >  1:951555/1‑140 (MQ=255)
                                                                                    tGTCCGCGACAACAAAACGCAGTCCCCCTTCTGTCTCGGTGTAATACCGTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                   <  1:509694/140‑1 (MQ=255)
                                                                                                caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa                                       >  2:738579/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  2:877102/1‑140 (MQ=255)
                                                                                                          tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                             >  2:301937/1‑140 (MQ=255)
                                                                                                                         ggTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  1:62070/1‑140 (MQ=255)
                                                                                                                         ggTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  1:519955/1‑140 (MQ=255)
                                                                                                                         ggTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  2:292139/1‑140 (MQ=255)
                                                                                                                         ggTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct              >  2:341697/1‑140 (MQ=255)
                                                                                                                          gtgtAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc             >  1:44315/1‑140 (MQ=255)
                                                                                                                           tgtAATACTTTTTTTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct            <  1:183399/140‑1 (MQ=255)
                                                                                                                              aataatTTTTGTTTTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc         <  1:877102/140‑1 (MQ=255)
                                                                                                                               atatttttttttttC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc        <  2:775014/134‑1 (MQ=255)
                                                                                                                               ataaCGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc        >  1:829519/1‑140 (MQ=255)
                                                                                                                                   ttttttttttT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  2:22402/138‑1 (MQ=255)
                                                                                                                                   ttttttttttT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  2:62070/138‑1 (MQ=255)
                                                                                                                                   gttttgtttTC‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  2:1076889/138‑1 (MQ=255)
                                                                                                                                   cGTTTTTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg    <  1:301937/140‑1 (MQ=255)
                                                                                                                                     tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  2:309606/1‑140 (MQ=255)
                                                                                                                                     tttgtttTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa  >  1:658659/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTT‑CCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGA  >  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: ‑