Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA33.3% 62.8 / 15.1 39intergenic (‑124/+79)ves/spyHutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (24/2);  new base A (0/13);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 1.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACG  >  NC_000913/1824922‑1825193
                                                                                                                                            |                                                                                                                                    
tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                                                      >  1:1068161/1‑140 (MQ=255)
           tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATTGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:657965/1‑140 (MQ=255)
           tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:99309/1‑140 (MQ=255)
                aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                      >  1:1206294/1‑140 (MQ=255)
                  tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                    >  1:1282751/1‑140 (MQ=255)
                     cccccGGCCCTTTTTTTTTTTTGTTATCGTTCAGGTAACGGTGGTTTTTTATCGGCAACTGGGGTGATTTACGTTGGTGGTGATCAGTTTTTTTTTAAACTTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                                                                                                 <  2:657965/138‑1 (MQ=255)
                     cccccGGCCCGTTTTTTTCATTGTTATCGTTCAGGTACCGGTGGTTGTTTACCGGCAACTGGGGTTTTTTAGGTTAGTGGTGATCAGTTTTTTTTTACAATTTTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCAc                                                                                                                 <  2:1206294/138‑1 (MQ=255)
                     cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                 >  2:443796/1‑140 (MQ=255)
                          ggCCCTTTTTTTTCTTTTTTATCGTTCAGGTACCGGTGGATGTTTATCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                            <  2:1282751/140‑1 (MQ=255)
                          ggCCCTTTTTTTTCTTTGTTTTCGTTCAGGTACCGGTGGATGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTTTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                            <  2:1053572/140‑1 (MQ=255)
                          ggCCCGTTCTTTTCTTTGTTATCGTTCAGGTACCGGTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                            <  2:99309/140‑1 (MQ=255)
                          gccccTTTTTTTTTTTTGTTTTTGTCGAGGTAACCGTGGATGTTTCTGGGAAACTGGGTTTATTTACGTTGGTGGTAATCAGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                            <  2:1271871/138‑1 (MQ=255)
                                ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                      >  1:274555/1‑140 (MQ=255)
                                ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                      >  1:1271871/1‑140 (MQ=255)
                                 ttttttCTTTTTTATCGTTCAGGTACCGGTGGTTGGTTATGGGCAACTGGGGTGTTTTATGTTGGGGGTGTTCAGGTTTTTTTTACAATTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                     <  2:274555/138‑1 (MQ=255)
                                 ttttttCTTTGTTATGGTCCAGGTACCCGTGGGTGTTTAGCGGCACCTGGGTTGTTTTGGGTTAGTGGGGATCCGGCTTTTTTAACAGGTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                     <  2:1406792/138‑1 (MQ=255)
                                                    cAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                  >  1:537075/1‑140 (MQ=255)
                                                      ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                >  1:1156480/1‑140 (MQ=255)
                                                            cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCttt                                                                          >  2:843690/1‑139 (MQ=255)
                                                                 aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt                                                                     >  1:1480455/1‑140 (MQ=255)
                                                                 aTGGTTATCGGAAACTGGGGTGTTTTGCTTTGGTGGTGATCGGTTTTTTTTTACAACTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt                                                                     <  2:1480455/140‑1 (MQ=255)
                                                                    tttATCGGCAACTGGGGTGTTTTACGTTGGTGGTGATCAGGTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                  <  2:330850/139‑1 (MQ=255)
                                                                    gTTATCGGAACCTGGGGTGATTTCGTTTAGTGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                  <  2:587071/140‑1 (MQ=255)
                                                                     ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTg                                                                 >  1:376145/1‑140 (MQ=255)
                                                                       ttCGGCAACTGGGGTGATTTACGTTGGGGGTGATCAGGTTTTTTTTAAAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                               <  1:1290540/139‑1 (MQ=255)
                                                                        tCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                                                              >  2:37875/1‑140 (MQ=255)
                                                                               cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                       >  1:63836/1‑140 (MQ=255)
                                                                               cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                       >  2:1416632/1‑140 (MQ=255)
                                                                                 ggggttttttAGGTTAGGGGTGAGGGGTTTTTTTTTAAAATTGTTGGAGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                     <  2:883849/133‑1 (MQ=255)
                                                                                         ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGGAAGATTTTTcagcag                                             >  1:883849/1‑140 (MQ=255)
                                                                                              ttAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                        >  1:673922/1‑140 (MQ=255)
                                                                                               tAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                       <  1:37875/140‑1 (MQ=255)
                                                                                                 gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                     >  2:1642274/1‑140 (MQ=255)
                                                                                                  tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCa                                    >  2:1097409/1‑140 (MQ=255)
                                                                                                      gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCAtttt                                >  1:1023264/1‑140 (MQ=255)
                                                                                                      gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt                                >  1:905285/1‑140 (MQ=255)
                                                                                                       aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                               >  2:161995/1‑140 (MQ=255)
                                                                                                            gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                          >  1:1254783/1‑140 (MQ=255)
                                                                                                                      cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa                >  1:179757/1‑140 (MQ=255)
                                                                                                                        tcTGTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg              <  2:1023264/139‑1 (MQ=255)
                                                                                                                         ttGTTAGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGt             <  2:179757/139‑1 (MQ=255)
                                                                                                                          tGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt            <  2:673922/140‑1 (MQ=255)
                                                                                                                                   aaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa    <  2:1156480/135‑1 (MQ=255)
                                                                                                                                     aaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg  <  2:376145/137‑1 (MQ=255)
                                                                                                                                       gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg              >  2:1279872/1‑126 (MQ=255)
                                                                                                                                       gaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg              <  1:1279872/126‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                    
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACG  >  NC_000913/1824922‑1825193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: