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

TCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGC  >  NC_000913/4111623‑4111891
                                                                                                                                  |                                                                                                                                           
tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGtt                                                                                                                                   >  2:1077090/1‑140 (MQ=255)
    aCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                                                                                               >  2:459322/1‑140 (MQ=255)
        tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                                                                                           >  1:784775/1‑140 (MQ=255)
          aaGCGATTCAGATAGCGCCAGCTTAATCGTTCCAACAGCGAAGGTCCCCCCTTTTT‑CGCCGTTGCCCGCGACAACATAACGCGGTGCCCCTTCTGTCTCGGTTTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACt                                                                                                                         <  2:832267/140‑1 (MQ=255)
              aaTTCAAATGGCCCCAGTTTAATCGTTCCAACAGGGAAGGCCACCCCTTTTC‑CGCCGTTGCCCGCAACAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAg                                                                                                                     <  2:733916/139‑1 (MQ=255)
                     aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                              >  1:616589/1‑140 (MQ=255)
                       aGCGCCATTTTAATCGGTTCAACAGGGAAGGCCCCCCCTTTTC‑CGCCTTTGTCCGCAACAAAATACCGCAGTCCCCCTTCTGTCTGGGTTAAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTg                                                                                                            <  2:616589/140‑1 (MQ=255)
                        gcgcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                                                                                           >  2:1329470/1‑140 (MQ=255)
                          gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                         >  1:148317/1‑140 (MQ=255)
                          gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                         >  1:663995/1‑140 (MQ=255)
                              gCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                                                                                     >  1:749202/1‑140 (MQ=255)
                                          cAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                         >  1:137681/1‑140 (MQ=255)
                                           aaCAGCGAAGGTCACCCCCTTTT‑CGCCGTTGCCCGCGACAAAATAACGCAGTGCCCTTTGTGTCTGGGTGTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                                                                                        <  2:784775/140‑1 (MQ=255)
                                                  aagcacAGCCCCTTCG‑CGCGGTTCGCCGCACCAAAATCGCGGGGTGCACCTGCTGCGTGGGTTTAATAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                                 <  2:512306/135‑1 (MQ=255)
                                                  aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                                 >  1:447719/1‑140 (MQ=255)
                                                   agggccaccccTTTTCCGCCTTTGCCCGCAAAAAAATAACGCGGTGCACTTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                                 <  2:399253/133‑1 (MQ=255)
                                                    ggTCAGCCCCTTTT‑CGCCGTTGTCAGCAACAACATAACGGAGTGCTCCTTCTGTCTCGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa                                                                               <  2:471213/140‑1 (MQ=255)
                                                            cttttt‑CGCGGTTGCCCCCAAAAACATACCGCGGTCCCCTTTCTTTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg                                                                       <  2:954381/138‑1 (MQ=255)
                                                            cttttc‑CGCCTTTGCCCGCAAAAACATAACGCAGTGCCCTTTCTTTCTGGGTATAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg                                                                       <  2:669846/138‑1 (MQ=255)
                                                            ccTTTT‑CGCCGTTGCCCGCAACAACATAACGCGGTCCCCTTTCTGTCTGGGTGTAATATCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg                                                                       <  2:750106/140‑1 (MQ=255)
                                                               ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                                    >  1:996109/1‑140 (MQ=255)
                                                               ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                                    >  1:177578/1‑140 (MQ=255)
                                                                tt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATCACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGAGTGCGGATCa                                                                   >  1:382004/1‑140 (MQ=255)
                                                                 t‑CGCCGTTGTCCGCAACAACATAACGGAGTGCCCTTTCTGTCTGGGTGTAATATCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCAc                                                                  <  2:447719/140‑1 (MQ=255)
                                                                       gTTGCCCGCAACAAAATACCGCATTCCCCCTTCTGTCTGGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                             <  2:996109/140‑1 (MQ=255)
                                                                          gTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCa                                                          >  2:962141/1‑140 (MQ=255)
                                                                             cgcgACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGt                                                       >  2:93385/1‑140 (MQ=255)
                                                                                               tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                     >  2:698032/1‑140 (MQ=255)
                                                                                               tGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                     >  1:479815/1‑140 (MQ=255)
                                                                                                  cctttCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTTTTCCTGAGctttcgctttcg                                  <  2:1112840/137‑1 (MQ=255)
                                                                                                    cTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCAGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                                <  1:390510/140‑1 (MQ=255)
                                                                                                      tCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGc                              >  2:226336/1‑140 (MQ=255)
                                                                                                              ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTTTTCCTGAGCTTTCGCTTTCGCGGCTTTTTGcc                      >  1:1112840/1‑139 (MQ=255)
                                                                                                              ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTGTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCGGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                      >  2:938542/1‑140 (MQ=255)
                                                                                                                tgtAATAGTGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                    <  1:1131804/140‑1 (MQ=255)
                                                                                                                tgtAATAGCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                    <  1:1329470/140‑1 (MQ=255)
                                                                                                                tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct                    >  1:273075/1‑140 (MQ=255)
                                                                                                                   aataatTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGc                 <  2:439925/140‑1 (MQ=255)
                                                                                                                    atacttttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcc                <  2:479815/134‑1 (MQ=255)
                                                                                                                        ctttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg            <  1:431908/138‑1 (MQ=255)
                                                                                                                         gtttgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTgg           <  1:226336/140‑1 (MQ=255)
                                                                                                                          tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa          >  2:320031/1‑140 (MQ=255)
                                                                                                                          tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa          >  1:895777/1‑140 (MQ=255)
                                                                                                                            tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg        >  2:18494/1‑140 (MQ=255)
                                                                                                                            tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg        >  1:723724/1‑140 (MQ=255)
                                                                                                                            tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg        >  1:375849/1‑140 (MQ=255)
                                                                                                                             ttttttCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg       <  2:177578/139‑1 (MQ=255)
                                                                                                                             gtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg       <  1:320031/140‑1 (MQ=255)
                                                                                                                                  tCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGc  <  2:723724/140‑1 (MQ=255)
                                                                                                                                  tCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGc  <  2:749202/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                           
TCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGGGTAGC  >  NC_000913/4111623‑4111891

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