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

TCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111615‑4111885
                                                                                                                                     |                                                                                                                                          
tcttttATACAACCCTCCAAGCGTTTCAGATAGCGCCAGCTTAATCGGTAAAACAGCGAGGGCCCCCCCTTTTC‑CGCCGTTGCCCGCGAAAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAA‑TTTTTTTTTTccc                                                                                                                                    <  2:447093/137‑1 (MQ=255)
  ttttATTCAACCCTTCAAGCGATTCAGATAGCGCCAGCTTATTCGGTCCAACAGCGAGGGCCAGCCCCTTTC‑CGCCGTTGTCCGCGACAACATAACGCGGTGCCCTTTTTGTCTCGGTGTAATAACTTTTTTTTTTcccc                                                                                                                                   <  2:730112/139‑1 (MQ=255)
                tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                                                                                     >  1:155608/1‑140 (MQ=255)
                 cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                    >  1:1028875/1‑140 (MQ=255)
                   aGCGTTTCAGATAGCGCCATTTTAATCGGTTCAACAGCGAAGGTCCCCCCTTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAATGTTTTTTTTCCCCCGCCGTTAGTAGCGACTg                                                                                                                  <  2:1291533/140‑1 (MQ=255)
                        ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  1:1291533/1‑140 (MQ=255)
                        ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                             >  1:22444/1‑140 (MQ=255)
                            gATAGCCCCAGCTTAATCGGTCAAACAGCGAGGGCCACCCCTTTTC‑CGCCGTTGCCCGCAAAAAAATAACGCAGTCCCCCTTCTGTCGGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                         <  2:155608/140‑1 (MQ=255)
                            aaTAGCCCCAGCTAAACGGTTTCAACAGAGAGGGTCACCCCTTTTG‑CGCCGTTTCCCGAGAAAACAAACCGCAGTGCCCTTTCTGTCCGGTTTTAATAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGtttttt                                                                                                         <  2:22444/139‑1 (MQ=255)
                             aTAGCGCCAGCTTAATCGGTACAACAGCGAAGGTCACCCCTTTTC‑CGCCGTTGCCCGCGACAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                        <  2:1406858/140‑1 (MQ=255)
                              tggccccAGTTAAATGGGTAAAACAGGGAAGGCCCCCCCTTTTT‑CCCCGTTGCCCGCGACAACATACCGCGGTGCCCCTTCTGTCTGGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCt                                                                                                       <  2:1109222/135‑1 (MQ=255)
                                 cgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                                                                                    >  2:153233/1‑140 (MQ=255)
                                  gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  1:895476/1‑140 (MQ=255)
                                  gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   >  1:943270/1‑140 (MQ=255)
                                  cccAGTTTAATCGTTTCAAGGGCGAAGGTCACCCCTTTTT‑CGCGGTTTCCCGGCAAAAAATAACCCAGTGCCCCTTCTTTCTGGGTTTAATAGTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                   <  2:1144297/139‑1 (MQ=255)
                                       cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                              >  2:87605/1‑140 (MQ=255)
                                         tAATCGTTAAAACAGCGAGGTTCACCCCTTTTG‑CCCCGTTGCCCGAAAAAACAAACCGCAGTCCCCTTCCTGACGCGGTGTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtt                                                                                            <  2:1521325/140‑1 (MQ=255)
                                          aaTCGGTCCAACAGCGAAGTCCCCCCCCTTTC‑CGCCGTTGTCCGCAAAAACATACCGCAGTCCCCTTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                                                                           <  2:506397/140‑1 (MQ=255)
                                           aTCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                                                                                          >  2:237278/1‑140 (MQ=255)
                                              ggTTCAACAGATACGGTCCCCCCCTTTG‑CGCGGTTTGCCGCGACAACAAACCGCGGGCCTCTTTCTGTTGCGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                       <  2:943270/140‑1 (MQ=255)
                                                   aaCAGCGAAGGTCACCCCTTTTT‑CGCCGTTGTCCGCGACAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                                                                                  <  2:942138/140‑1 (MQ=255)
                                                         gAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg                                                                            >  1:942138/1‑140 (MQ=255)
                                                           agggccaccccTTTTCCGCCTTTGCCCGCAAAAAAATACCGCGGTCCCCTTTCTGTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                           <  2:528694/133‑1 (MQ=255)
                                                            ggccccccctTTTT‑CGCGGTTGTCCGCAACAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATATCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                         <  2:525452/134‑1 (MQ=255)
                                                                       ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                              >  2:321088/1‑140 (MQ=255)
                                                                       ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                              >  2:251785/1‑140 (MQ=255)
                                                                               tttGCCCGCGACAACATACCGCAGTCCCCCTTCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                       <  2:895476/139‑1 (MQ=255)
                                                                                 tGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                     >  2:1017450/1‑140 (MQ=255)
                                                                                  gTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCa                                                    >  2:35953/1‑140 (MQ=255)
                                                                                    ccGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                  <  1:251785/140‑1 (MQ=255)
                                                                                          caacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTcc                                            >  1:214718/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGTGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  1:766418/1‑140 (MQ=255)
                                                                                                       tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                               >  2:438866/1‑140 (MQ=255)
                                                                                                              tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                        >  2:1331860/1‑140 (MQ=255)
                                                                                                                      ggTGTAATATCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                <  1:680692/140‑1 (MQ=255)
                                                                                                                        tgtAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:662523/140‑1 (MQ=255)
                                                                                                                        tgtAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGTATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct              <  1:87605/140‑1 (MQ=255)
                                                                                                                          taataaCGTTTGTTTTTCCCCGCCGTTAGTAGTGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg            >  1:172207/1‑140 (MQ=255)
                                                                                                                          taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg            >  2:408031/1‑140 (MQ=255)
                                                                                                                           aataaCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc           <  1:153233/140‑1 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:1240284/1‑140 (MQ=255)
                                                                                                                                  tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa    >  1:1006073/1‑140 (MQ=255)
                                                                                                                                    tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt   <  1:477949/139‑1 (MQ=255)
                                                                                                                                    tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGAt   >  2:477949/1‑139 (MQ=255)
                                                                                                                                    tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg  >  2:215057/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                          
TCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG  >  NC_000913/4111615‑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: ‑