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

CAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGG  >  NC_000913/4111624‑4111887
                                                                                                                                  |                                                                                                                                       
cAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                                                                                              >  1:1346431/1‑140 (MQ=255)
  aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCTGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagt                                                                                                                            >  1:851448/1‑140 (MQ=255)
   aCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagtc                                                                                                                           >  1:1460697/1‑139 (MQ=255)
        cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                      >  1:1447699/1‑140 (MQ=255)
        cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                      >  1:518571/1‑140 (MQ=255)
         aaGCGATTCAGATAGCCCCAGCTTAATCGGTTCAACAGCGAAGGTCCCCCCTTTT‑C‑CGCCGTTGCCCGCAAAAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACt                                                                                                                     <  2:1346431/140‑1 (MQ=255)
               ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                                                                                               >  1:210824/1‑140 (MQ=255)
                         cccAGTTTAATGGGTCAAACAGGGAGGGCCCCCCCTTTT‑C‑CGCGGTTGCCCGCGACAACATACCGCGGTCCCCCTTCTGTCTGGGTTTAATACTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                     <  2:518571/139‑1 (MQ=255)
                         cccAGTTTAATCGGTCAAACAGCGAGGGTCACCCCTTTT‑C‑CGCCGTTGCCCGCAAAAACAAACCGCAGTCCCCCTTCTGTCTGGGTATAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                     <  2:188983/139‑1 (MQ=255)
                          ccAGTTTAATCGTTTCAACAGCGAGGGTCGCCCCTTTT‑T‑CGCCTTTGTCCGCGAAAACATACCGCGGTGCCCCTTCTGTCTGGGTGAAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                    <  2:230020/140‑1 (MQ=255)
                          ccAGCTTAATGGGTCCAACAGCGAAGGCCACCCCTTTT‑C‑CGCCTTTGCCCGCAAAAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCtt                                                                                                    <  2:210824/140‑1 (MQ=255)
                           cAGTTTAATCGGTAAAACAGC‑AAGGGCCCCCCCTTTCT‑CGCCTTTGCCCGCAAAAAAATAACGCAGTCCCCCTTTTTTCTCGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                                                                                   <  2:1541336/140‑1 (MQ=255)
                              cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                                >  1:931631/1‑140 (MQ=255)
                                   tgggTTCAACAGGGAAGGGCACCCCTTTT‑T‑CCCGGTTGTCCGCGACAACATAGGGGAGTGCTCCTTTTTTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                           <  2:1459451/138‑1 (MQ=255)
                                   cagtttCAACAGCGAGGGTCACCCCTTTT‑C‑CGCGTTTGCCCGCAACAACATACCGCAGTGCCCCTTCTGTCTGGGTTTAATATGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                           <  2:1105018/136‑1 (MQ=255)
                                   aaggTACAACAGCGAAGGTCAGCCCCTTT‑C‑CGCCGTTGCCCGCGACAAAATAACGCAGTCCACCTTCGGTCTGGGTAAAAAAACTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                                                                                           <  2:1068095/138‑1 (MQ=255)
                                     gtttaaaCGGCGAAGTTCCCCCCCTTT‑T‑CGCCGTTGCCCGCGACAACATACCGCAGTGCCCCTTCTGTCTCGTTGAAATACCTTTTTTTTCTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                         <  2:16656/135‑1 (MQ=255)
                                     gttaaaaCAGCGAAGGTCACCCCTTTT‑C‑CGCTGTTGCCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                         <  2:1014281/135‑1 (MQ=255)
                                       ttCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                       >  1:324940/1‑140 (MQ=255)
                                        tCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt                                                                                      >  1:639187/1‑140 (MQ=255)
                                         cAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                     >  1:958007/1‑140 (MQ=255)
                                               cGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                                                                               >  2:581984/1‑140 (MQ=255)
                                                 aaGGTCCCCCCTTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                             <  2:407414/140‑1 (MQ=255)
                                                 aaGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggggg                                                                             >  1:318373/1‑139 (MQ=255)
                                                 aaGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggctg                                                                             >  1:1105018/1‑137 (MQ=255)
                                                 aaGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggctg                                                                             >  1:1068095/1‑137 (MQ=255)
                                                 aaGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                             >  2:41387/1‑140 (MQ=255)
                                                 aaGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                             >  1:1176254/1‑140 (MQ=255)
                                                  aaggcccccccTTT‑TCCGCCGTTGTCCGCAAAAACATAACGCAGTCCCCTTTCTTTCTGGGTTAAATAGCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                             <  2:218178/133‑1 (MQ=255)
                                                       aGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgtgt                                                                       >  2:650796/1‑140 (MQ=255)
                                                                   cgcgGTTGTCCGCGACAACATAACGCGGTCCCCCTTCTGTCTGGGTTAAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAcc                                                             <  2:334260/140‑1 (MQ=255)
                                                                       tttGTCCGCGACAACATACCGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                         <  1:224623/139‑1 (MQ=255)
                                                                            ccGAGAAAAAATACCGCAGTCCCCTTTCTGTCTGGGTTTAATATCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg                                                    <  2:860262/140‑1 (MQ=255)
                                                                                          cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                      >  2:422519/1‑140 (MQ=255)
                                                                                             aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                                   >  1:410312/1‑140 (MQ=255)
                                                                                             aGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgc                                   >  1:845753/1‑140 (MQ=255)
                                                                                               tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                 >  2:157326/1‑140 (MQ=255)
                                                                                                gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                                >  1:632318/1‑140 (MQ=255)
                                                                                                  aCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                              >  2:1117320/1‑140 (MQ=255)
                                                                                                  aCCTTCTGTCGGGGTATAACAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGctttcgctttcg                              <  1:137782/140‑1 (MQ=255)
                                                                                                    tttCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                            <  2:713687/139‑1 (MQ=255)
                                                                                                         gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                       >  1:735302/1‑140 (MQ=255)
                                                                                                         gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                       >  2:430432/1‑140 (MQ=255)
                                                                                                              ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                  >  1:1289070/1‑140 (MQ=255)
                                                                                                              ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                  >  2:189983/1‑140 (MQ=255)
                                                                                                               gtgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc                 >  1:713687/1‑140 (MQ=255)
                                                                                                                tgtATTACCGTTTTTTTTACCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                <  2:623065/140‑1 (MQ=255)
                                                                                                                tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                <  1:41387/140‑1 (MQ=255)
                                                                                                                  taaaaaCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg              <  2:410312/140‑1 (MQ=255)
                                                                                                                   aataatTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc             <  2:1405709/140‑1 (MQ=255)
                                                                                                                     taaCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGccc           <  2:31271/140‑1 (MQ=255)
                                                                                                                        cGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg        <  1:650796/140‑1 (MQ=255)
                                                                                                                          tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa      >  1:632046/1‑140 (MQ=255)
                                                                                                                          tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa      >  2:1167379/1‑140 (MQ=255)
                                                                                                                             ttttttCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg   <  1:189983/139‑1 (MQ=255)
                                                                                                                             gtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg   <  1:1117320/140‑1 (MQ=255)
                                                                                                                              tttttCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATggg  >  2:1483027/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTT‑T‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGG  >  NC_000913/4111624‑4111887

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