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

TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTC  >  NC_000913/1824922‑1825197
                                                                                                                                           |                                                                                                                                        
tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                                                          >  1:263430/1‑140 (MQ=255)
  tCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                        >  1:728524/1‑140 (MQ=255)
           tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                               >  2:432668/1‑140 (MQ=255)
           tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                               >  2:759090/1‑140 (MQ=255)
                aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                          >  2:583973/1‑140 (MQ=255)
                  ccGCCCCCGGCCCTTTTTTTTCTTTTTTATCTTTCAGGTAACGGTGGATGGTTATGGGCAACTGGGGTGTTTTACGTTGGTGGTGATCAGTTTTTTTTAACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTc                                                                                                                        <  2:580127/139‑1 (MQ=255)
                    gCTCCCGGTCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCa                                                                                                                      >  2:962562/1‑140 (MQ=255)
                    gCTCCCGGCCCTTTTTTTTCTTTTTTATCGTTAAGGTACCCGTGGTTGTTTATCGGCAACGGGGGTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAATTTTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCa                                                                                                                      <  2:637757/140‑1 (MQ=255)
                           gCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                               >  2:718063/1‑140 (MQ=255)
                            cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                              >  2:819229/1‑140 (MQ=255)
                            cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAG‑AAACCGCCGATCCTGTCCACCGCATTAc                                                                                                             >  1:477223/1‑140 (MQ=255)
                            cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAG‑AAACCGCCGATCCTGTCCACCGCATTAc                                                                                                             >  1:903432/1‑140 (MQ=255)
                                 ttttttCTTTTTTATCGTTAAGGTACCCGTGGATGTTAATCGGCAACTGGGTTGATTTACGTTGGTGGTGACCAGGTTTTTTTTACACTTGTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                         <  2:17257/138‑1 (MQ=255)
                                       tttttttATCGTTCAGGTAACCGTGGTTGTTTACCGGCAACTGGGGTTTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAACTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                                   <  2:293051/135‑1 (MQ=255)
                                         ttGTTATCGTTCAGGTAACGGTGGTTGGTTATGGGCACCTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                 <  2:955303/140‑1 (MQ=255)
                                                    cAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                      >  1:774546/1‑140 (MQ=255)
                                                    agggTAACGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGGGGTAACGGGTTTTTTTTTACAGTTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                      <  2:746687/138‑1 (MQ=255)
                                                      ggTAACGGTGGTTGGTTATCGGAAACTGGGGTTTTTTAGGTTGGTGGTGATGAGTTTTTTTTTAAAATTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                    <  2:263430/140‑1 (MQ=255)
                                                            ggTGGAGGGTTACCGGCACCTGGGGTGATTTAGTTTAGGGGTGATCAGGTTTTTTTTACAACTTTTGAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                              <  2:175949/139‑1 (MQ=255)
                                                                     ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCCTAAGTTTTTTg                                                                     >  1:775015/1‑140 (MQ=255)
                                                                       ttcgcaaaCGCGGGTGTTTTACTTTAGTGGTAATCGGCTTTTTTTTAAAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                   <  1:962562/134‑1 (MQ=255)
                                                                                              ttggtggtgATCGGGTTTTTTTTCCATGTTTTTGAGAAAAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                            <  1:96451/137‑1 (MQ=255)
                                                                                                 gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                         >  1:227757/1‑140 (MQ=255)
                                                                                                    gtgATCGGGTTTTTTTTAAAACTGTTTGAGAAGAAAAAAG‑AAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                     <  2:477223/140‑1 (MQ=255)
                                                                                                    gtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAG‑AAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                     <  2:903432/140‑1 (MQ=255)
                                                                                                       aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                   >  2:12713/1‑140 (MQ=255)
                                                                                                           ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCA‑GCATTTTACCttt                              >  2:566708/1‑140 (MQ=255)
                                                                                                            gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                              >  1:475538/1‑140 (MQ=255)
                                                                                                              tttttttACATTTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                            <  1:356410/140‑1 (MQ=255)
                                                                                                                    tACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                      >  1:515703/1‑140 (MQ=255)
                                                                                                                          ttttagagaagagaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                <  1:759090/135‑1 (MQ=255)
                                                                                                                          tGTTAGAGAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                <  2:728524/140‑1 (MQ=255)
                                                                                                                                  aaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa        <  2:764820/135‑1 (MQ=255)
                                                                                                                                     aaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc     <  1:139566/138‑1 (MQ=255)
                                                                                                                                      aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCA‑GCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt   <  1:566708/139‑1 (MQ=255)
                                                                                                                                        aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGATAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc  <  2:774546/140‑1 (MQ=255)
                                                                                                                                        aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc  <  1:432668/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTC  >  NC_000913/1824922‑1825197

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