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

TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT  >  NC_000913/1824922‑1825196
                                                                                                                                             |                                                                                                                                       
tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                                                         >  1:1460960/1‑140 (MQ=255)
      aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                   >  1:1496855/1‑140 (MQ=255)
           tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                              >  1:1456943/1‑140 (MQ=255)
           tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                              >  1:880655/1‑140 (MQ=255)
             cATCCTCGCTCCCGGCCTTTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAGAAAACCGCCGATc                                                                                                                            <  1:1359582/140‑1 (MQ=255)
               tccTCGCTCCCGGCCCGTTCTTTTCATTTTTATCGTTCAGTTACCCGTGGATGGTTATCGCCACCTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAACCCGCCGATCCt                                                                                                                          <  2:42500/138‑1 (MQ=255)
                aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                         >  2:924154/1‑140 (MQ=255)
                aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                         >  1:275082/1‑140 (MQ=255)
                 cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt                                                                                                                        >  2:1538259/1‑140 (MQ=255)
                  tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                       >  1:893618/1‑140 (MQ=255)
                     cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTGGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                    >  2:1141734/1‑140 (MQ=255)
                       cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg                                                                                                                  >  1:888145/1‑140 (MQ=255)
                       cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                 >  1:658671/1‑140 (MQ=255)
                         cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCa                                                                                                                >  1:1115176/1‑140 (MQ=255)
                          ggccttttCTTTTTATTGTTATTGTACAGATAATGGTGTATGTTTATCGTCATTGGGGTTATTTTACGTTATTGTTAAGCAGGTTTTTTTAACAAGTGTTAGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                               <  2:1115176/134‑1 (MQ=255)
                          ggCCCGTTTTTTTTATTTTTTATCTTTAAGGTACCGGTGGATGTTTATCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCAGTTTTTTTTCAAACCTTTTAGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCa                                                                                                                <  2:315523/140‑1 (MQ=255)
                          ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGGTCTTGTCCActacct                                                                                                               >  1:864895/1‑135 (MQ=255)
                            cccTTTTTTTTCATTGTTATCGTTCAGGTAACGGTGGATGTTTATCGGAAACTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGGTCTTGTCCActaccttg                                                                                                             <  2:864895/140‑8 (MQ=255)
                            cccTTTCTTTTCATTGTTATCTTTCGGGTACCCGTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                             <  2:1333003/140‑1 (MQ=255)
                             cgttttttttCTTTTTTATCGTTCAGGTAACGGGGGATGTTTATCGGCAACTGGGGTTATTTACTTTAGTGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTAc                                                                                                            <  2:977065/134‑1 (MQ=255)
                                       tcttttttATCGTTCAGGAAACGGTGGTTGGTTATCGGCACCTGGGGTGTTTTACTTTAGGGGTGATCGGTTTTTTTTTACATTTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                   <  2:1460960/134‑1 (MQ=255)
                                       tatttGTTTTCGTCCAGGTAACGGGGGATGTTTATCGGAATCGCGGGTGTTTTAGTTTGGTGGTGATCGGGTTTTTTTAACATTTTTTGAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                   <  2:880655/137‑1 (MQ=255)
                                              tATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                            >  1:309363/1‑140 (MQ=255)
                                                   ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                       >  1:1273716/1‑140 (MQ=255)
                                                       ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCATCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                   >  2:1528205/1‑140 (MQ=255)
                                                                   gggTTACCGGCAACTGGGGTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACATCTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt                                                                       <  2:717706/139‑1 (MQ=255)
                                                                        aTCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                  >  1:1250327/1‑140 (MQ=255)
                                                                                cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAgct                                                          >  1:717706/1‑138 (MQ=255)
                                                                                          tttACGTTGGTGGTGATCGGTTTTTTTTAACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                 <  2:275082/140‑1 (MQ=255)
                                                                                                tttgtggtgCTCGGTTTTTTTTCACAGTTCTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGTATTTCTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                           <  2:1456943/137‑1 (MQ=255)
                                                                                                ttAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                           <  2:1250327/140‑1 (MQ=255)
                                                                                                 tggtggtcACCAGTTTTTTTTAACAATTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                          <  2:1496855/138‑1 (MQ=255)
                                                                                                 tAGGGGTGATCAGTTTTTTTTAAAAATTGTTGAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                          <  2:309363/140‑1 (MQ=255)
                                                                                                   gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                        >  2:1112609/1‑140 (MQ=255)
                                                                                                      gtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCCGGCAtt                                     >  2:88138/1‑140 (MQ=255)
                                                                                                        gATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt                                   >  1:1334472/1‑140 (MQ=255)
                                                                                                         aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                  >  1:117666/1‑140 (MQ=255)
                                                                                                         aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCATAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTTCAGGCAttttt                                  >  2:815011/1‑139 (MQ=255)
                                                                                                                tttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                           >  1:1248298/1‑140 (MQ=255)
                                                                                                                   ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg                        >  2:588567/1‑140 (MQ=255)
                                                                                                                   ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg                        >  1:1551825/1‑140 (MQ=255)
                                                                                                                      aaaatcTGTTAGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                     <  2:1551825/135‑1 (MQ=255)
                                                                                                                              tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt             >  2:158454/1‑140 (MQ=255)
                                                                                                                                        gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt               >  1:1060416/1‑128 (MQ=255)
                                                                                                                                        gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt               <  2:1060416/128‑1 (MQ=255)
                                                                                                                                        gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAAATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGACGCt   >  2:1540594/1‑140 (MQ=255)
                                                                                                                                         aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt  >  1:628199/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                       
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT  >  NC_000913/1824922‑1825196

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