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

GTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCT  >  NC_000913/1824929‑1825195
                                                                                                                                     |                                                                                                                                      
ttATCCCATCCTCCCCCCCGCCCCTTTTTTTTCATTTTTATCGTTCAGGTAACCGTGGATGTTTATCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAAAAGAAAAAAAAAAAccg                                                                                                                                 <  2:489097/139‑1 (MQ=255)
    tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                             >  1:1056994/1‑140 (MQ=255)
    tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                             >  2:1329375/1‑140 (MQ=255)
    tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                             >  1:517927/1‑140 (MQ=255)
    tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                             >  1:381289/1‑140 (MQ=255)
       aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc                                                                                                                          >  2:177269/1‑140 (MQ=255)
           tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                      >  1:1345725/1‑140 (MQ=255)
             gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                   >  1:574842/1‑140 (MQ=255)
              cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                   >  2:177278/1‑140 (MQ=255)
              cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                   >  1:1683960/1‑140 (MQ=255)
              cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                   >  2:710920/1‑140 (MQ=255)
                  cgccccTTTTTTTTATTTTTTTTCTTTCAGGTACCGGTGGTTGTTTATCGGCAGCTGGGTTATTTTCTGTTGGTGGTGATCGGTTTTTTTTAACAATTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCa                                                                                                               <  2:1056994/137‑1 (MQ=255)
                  cGGCCCTTTTTTTTCATTGTTATCGTTCAGGTACCGGTGGATGGTTACCGGCAACTGGGGTGATTTACGTTAGTGGTGATGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                              <  2:574842/140‑1 (MQ=255)
                  cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCa                                                                                                               >  1:365327/1‑140 (MQ=255)
                  cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCa                                                                                                               >  1:1189385/1‑140 (MQ=255)
                         ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                        >  1:762794/1‑140 (MQ=255)
                         ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                        >  1:913029/1‑140 (MQ=255)
                          ttttttCATTTTTATCGTTCGGGTACCGGTGGATGTTTATCGGCACCGGGGGTGTTTTATGTTGGTGGTGATCGGGTTTTTTTTACAATTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                       <  2:157426/138‑1 (MQ=255)
                          ttttttCATTGTTATCGTTCAGGTACCGGTGGATGTTTATCGGCAACTGGGGTGTTTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                       <  2:466494/138‑1 (MQ=255)
                               tcttttttATCGTTAAGGTAACGGGGGTTGTTTATCGGCACCTGGGGTGTTTTACGTTAGTGGTGACCGGTTTTTTTTTAAACTTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                  <  2:626094/134‑1 (MQ=255)
                                       aTCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                                                                          >  1:1611182/1‑140 (MQ=255)
                                             cAGGTACCGGTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACACCTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                    <  2:1611182/140‑1 (MQ=255)
                                             aaGGTACCGGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                    <  2:517927/139‑1 (MQ=255)
                                              aGGTACCGGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGACCAGTTTTTTTTTACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                   <  2:1141609/140‑1 (MQ=255)
                                               gttaccggTGGATGGTTACCGGCACCTCGGGTGTTTTACGTAAGTGGTGACCGGGTTTTTTTAAAAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                  <  2:1618202/133‑1 (MQ=255)
                                                   aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTGGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACt                                                                              >  2:1209866/1‑140 (MQ=255)
                                                     aagggATGGTTTTCGGCAACTCGGGTGATTTACGTTAGTGGTGATCCGGTTTTTTTTAAAACTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                            <  2:65653/137‑1 (MQ=255)
                                                                  gggCACCTGGGGTGATTTCCGTTGGTGGTGATCGGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                               <  1:772060/139‑1 (MQ=255)
                                                                   ggCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGcctg                                                              >  1:65653/1‑137 (MQ=255)
                                                                        cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGGAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                         >  1:1224192/1‑140 (MQ=255)
                                                                                       gTTAGTGGTGATCAGTTTTTTTTTAAAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                           <  1:710920/140‑1 (MQ=255)
                                                                                           gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAAATTATTTAGCAGTTGCAggg                                       >  2:1381055/1‑139 (MQ=255)
                                                                                            tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCa                                      >  2:1498510/1‑140 (MQ=255)
                                                                                            tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCa                                      >  2:1362820/1‑140 (MQ=255)
                                                                                                 aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                 >  1:459357/1‑140 (MQ=255)
                                                                                                    aGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                              >  2:625181/1‑140 (MQ=255)
                                                                                                      tttttttttACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                            <  2:762794/138‑1 (MQ=255)
                                                                                                      gttttttttACAACTGTTTGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                            <  1:177278/138‑1 (MQ=255)
                                                                                                      gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                            >  2:710408/1‑140 (MQ=255)
                                                                                                        ttttttaACAATTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                          <  1:1381055/140‑1 (MQ=255)
                                                                                                        cgttttaACACTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                          <  1:1690272/138‑1 (MQ=255)
                                                                                                                    tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt              >  1:1701195/1‑140 (MQ=255)
                                                                                                                      tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCt            >  2:498045/1‑140 (MQ=255)
                                                                                                                      tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt            >  2:580708/1‑140 (MQ=255)
                                                                                                                      tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt            >  1:285789/1‑140 (MQ=255)
                                                                                                                      tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTGGCTGGAAGATTATTTAGCAGGTGCAGGCATTTTAACTTTTGCCGCTTGGCGtctg            >  2:1690272/1‑137 (MQ=255)
                                                                                                                            aaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa      <  1:720073/135‑1 (MQ=255)
                                                                                                                              gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCCGGACGGTCTGTCAGACg    >  2:720073/1‑140 (MQ=255)
                                                                                                                              gaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg    <  2:120256/137‑1 (MQ=255)
                                                                                                                               agaaaaaaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc   <  1:976244/140‑1 (MQ=255)
                                                                                                                                gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt  >  2:1613415/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                      
GTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCT  >  NC_000913/1824929‑1825195

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: