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

CAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACG  >  NC_000913/1824925‑1825193
                                                                                                                                         |                                                                                                                                    
cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                   >  1:720885/1‑140 (MQ=255)
   aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTACCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                >  1:1169329/1‑140 (MQ=255)
     ttttCCATCCTCCCTCCCGGCCTTTTTTTTTCTTTTTTATCTTCAGGGTACGCGGGGGTGGTTAGGGGCATCTGGGGTGTTTTACTTTAGGGGTGGTCAGTTTTTTTTTACACTTGTTTAAAAAGAGAAAAGAAAAccgc                                                                                                                              <  2:582679/138‑1 (MQ=255)
      aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                             >  1:1180036/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:993072/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:1129465/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:1322972/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:1340919/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  2:115828/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:1520562/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  2:1334250/1‑140 (MQ=255)
        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGCGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa                                                                                                                           >  1:109117/1‑140 (MQ=255)
             aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                      >  1:1116902/1‑140 (MQ=255)
               tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                    >  1:1081247/1‑140 (MQ=255)
                 gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTccc                                                                                                                  >  1:192581/1‑139 (MQ=255)
                  cccccGGCCCTTTTTTTTCTTTGTTCTCGTGCAGGAACCCGGGGATGGTTATGGGCAACTGGGGTGGTTTTGGTTGGGGGTGATCAGCTTTTTTTAACAAGTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                                                                                                 <  2:109117/138‑1 (MQ=255)
                       ggCCCTTTTTTTTCATTTTTATCGTTCAGGTACCGGGGGATGTTTATCGGCAACTGGGGTATTTTACTTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                            <  2:1116902/140‑1 (MQ=255)
                       ggCCCGTTTTTTTTATTTTTATCGTTCAGGTAACCGTGGATGGTTATCGGCACCTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                            <  2:1322972/140‑1 (MQ=255)
                             ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                      >  1:522745/1‑140 (MQ=255)
                             ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                      >  2:5609/1‑140 (MQ=255)
                              ttttttCTTTTTTATCGTTCAGGTACCGGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGGGGTGATCGGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                     <  2:662184/138‑1 (MQ=255)
                              ttttttCTTTGTTATCTTTCAGGTAACGGTGGTTGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGTTTTTTTTTACACTTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                     <  2:1129465/138‑1 (MQ=255)
                              ttttttCATTGTTATCGTTCAGGTACCGGTGGATGTTTACCGGCAACTGGGGTGTTTTACGTTGGTGGTGATCAGGTTTTTTTTACAACTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                     <  2:1340919/138‑1 (MQ=255)
                                   tCATTTTTATCGTTCAGGTACCGGTGGTTGTTTCTCGGCAACTCGGGTGTTTTACGTTGGTGGTGATCAGGCTTTTTTTACACCTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                <  2:1081247/140‑1 (MQ=255)
                                    taTTGTTATCGTTCAGGTAGCCGTGGATGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCGGGTTTTTTTAACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                               <  2:106059/139‑1 (MQ=255)
                                       tGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                            >  1:662184/1‑140 (MQ=255)
                                                    ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                >  2:1471135/1‑140 (MQ=255)
                                                          ggTGGATGGTTATCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                          <  2:1483858/139‑1 (MQ=255)
                                                                  tttATCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                  <  2:1520562/139‑1 (MQ=255)
                                                                   ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                                                 >  2:5950/1‑140 (MQ=255)
                                                                    tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                >  1:335568/1‑140 (MQ=255)
                                                                     aTCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCGGTTTTTTTTTACAACTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                               <  2:1127253/140‑1 (MQ=255)
                                                                                            ttAGTGTTGACCGGGTTTTTTTTACATCTGTTAGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                        <  2:1054863/140‑1 (MQ=255)
                                                                                             tAGTGGTGATGAGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                       <  1:1334250/140‑1 (MQ=255)
                                                                                              ggtggtgATCGGGTTTTTTTTAAAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                      <  1:5609/139‑1 (MQ=255)
                                                                                               gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCAGGc                                     >  2:838289/1‑140 (MQ=255)
                                                                                               gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                     >  1:514439/1‑140 (MQ=255)
                                                                                                tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa                                    >  1:1054863/1‑140 (MQ=255)
                                                                                                     aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                               >  1:242473/1‑140 (MQ=255)
                                                                                                       cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc                             >  1:629310/1‑140 (MQ=255)
                                                                                                                    cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa                >  2:450034/1‑140 (MQ=255)
                                                                                                                        tGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt            <  2:415457/140‑1 (MQ=255)
                                                                                                                          tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt          >  2:1147436/1‑140 (MQ=255)
                                                                                                                          tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTTCCGCTGGACGTTCt          >  2:823719/1‑140 (MQ=255)
                                                                                                                          tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTGCCGCTGGACGTTCt          >  1:1234631/1‑140 (MQ=255)
                                                                                                                          tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt          >  2:237484/1‑140 (MQ=255)
                                                                                                                          tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt          >  2:50058/1‑140 (MQ=255)
                                                                                                                                 agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc   <  1:1471135/140‑1 (MQ=255)
                                                                                                                                  aaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg  <  2:1169329/137‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACG  >  NC_000913/1824925‑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: