Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0310CT21.9% 66.8 / 10.6 32intergenic (‑94/+109)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base C (20/5);  new base T (0/7);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 2.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTC‑CGGCGTCAGGATGTTGTAAATTT  >  NC_000913/1824813‑1825252
                                                                                                                                                                                                                           |                                                                                                                                                                                                                               
gatgCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAc                                                                                                                                                                                                                        <  1:471320/228‑1 (MQ=255)
            gCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgag                                                                                                                                                                                                             <  2:571404/227‑1 (MQ=255)
                 ccAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TCTCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTga                                                                                                                                                                                                              >  2:640235/1‑221 (MQ=255)
                  cAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgag                                                                                                                                                                                                             <  1:641257/221‑1 (MQ=255)
                                         tGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACCCGGGTGATTTAAGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAATGAAAACCGCCGATCCTg                                                                                                                                                                                   >  1:570728/1‑224 (MQ=255)
                                                                                           tCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGGTTTTTgg                                                                                                                            >  1:373249/1‑228 (MQ=255)
                                                                                                                        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCCTTAAGTTTTTT‑GCTGGAAAACTATTCATCAGtttgagggc                                                                                                >  2:165510/1‑220 (MQ=255)
                                                                                                                               tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                                                                  >  1:137437/1‑155 (MQ=255)
                                                                                                                                cGCTCCCGCCCTGTTCTTTTCATTTTTATCTTTC‑TGTATCCTGGGGTGGTGTGTCGCCATCGAGAGTGTTTTCGGTAGGTGTTGACCGTTCTTTTTTAACAGCTGTTAAATAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                                                                                                                                 <  2:137236/155‑1 (MQ=255)
                                                                                                                                    cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTTGCTGaa                                                                                                                       >  1:550509/1‑192 (MQ=255)
                                                                                                                                             ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCCTAAGGTTTTTGGGTGAAAGAATATTTAGCAGTTGTAAGCATTTTAACTTTTGTCGCTGGGCg                                                                          >  2:184187/1‑228 (MQ=255)
                                                                                                                                               gtttgCATTGTTATCTTTGAGGTAGCGGTGGATGTT‑TATCGGCAACTGGGGTGATTTACGTTAGTGGTACTCAGGTTTTTTTTACAACTGTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                                                                        <  1:184465/227‑1 (MQ=255)
                                                                                                                                                          tATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                                                                                 >  1:366832/1‑97 (MQ=255)
                                                                                                                                                           aTCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTTGCTGACAGATTATTCAGCAGTTGCAGGCATTTTACCCTTTGCCGCTTGACGTTCTTTCAGACGCt                                                            >  1:692710/1‑228 (MQ=255)
                                                                                                                                                            tCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                   >  1:373844/1‑125 (MQ=255)
                                                                                                                                                            tCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTTGCTGAAAGATAATTCAGCCGTTGTCGGCGTTTTACCTTTTTCCCCTTGACGATCTGTCAGaatcat                                                           >  2:117860/1‑223 (MQ=255)
                                                                                                                                                                     tAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCa                                                                                                                                                                          >  1:77187/1‑109 (MQ=255)
                                                                                                                                                                         cGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                   >  2:354356/1‑112 (MQ=255)
                                                                                                                                                                         cGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                   >  1:647024/1‑112 (MQ=255)
                                                                                                                                                                          gTGGATGGT‑TAGCGGCAACTGGGGGGATTTACGTTGGTGGTGACCAGGTTTTTTTTAAACCTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                                                                  <  1:354947/112‑1 (MQ=255)
                                                                                                                                                                          gTGGATGGT‑TACCGGCAACTGGGGTGATTTACTTTGGTGGTGATCGGTTTTTTTTTAAAACTGTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                                                                  <  2:645985/112‑1 (MQ=255)
                                                                                                                                                                             gATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTCTTCAGCAGTTGCAGGCAATTTAc                                                                                          >  1:253070/1‑180 (MQ=255)
                                                                                                                                                                              ttgtt‑tATCGACACCTGGGGTGATTTACGTTAGTGGTAATCAGGTTTTTTTCACACTTTTTTGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc                                                                                         <  2:252679/176‑1 (MQ=255)
                                                                                                                                                                              aTGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                                                                               >  1:539109/1‑127 (MQ=255)
                                                                                                                                                                               tgtt‑tACCGGAAACTCGGGTAATTTACGTTAGGGGTGATCAGGTTTTTTTTACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                                                              <  2:538240/124‑1 (MQ=255)
                                                                                                                                                                                       cGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTa                                                                                                                                                                                                                         >  1:174831/1‑45 (MQ=255)
                                                                                                                                                                                        ggCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAc                                                                                                                                                                                                                        <  2:174565/45‑1 (MQ=255)
                                                                                                                                                                                                            ttAGTGGTGACCAGTTTTTTTTAACACCTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                                                                         <  2:77074/72‑1 (MQ=255)
                                                                                                                                                                                                               gtggttATCAGGTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTC‑CGGCGTCAGGAtgttgt        <  1:165761/228‑1 (MQ=255)
                                                                                                                                                                                                                     aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTC‑CGGCGTCAGGATGTTGTAAAttt  >  1:513421/1‑228 (MQ=255)
                                                                                                                                                                                                                       cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGCAAGATTATTTAGCAGTTGCAGGCCTTTTTCCTTTTTCCGCCGGACGGTCTGTTAGAACGTTTTCTAAAATAGACTTTAATTTCTTTTTTTTTTCTGGGCGCCGGGAGGTGGTAAAttt  >  2:540837/1‑227 (MQ=255)
                                                                                                                                                                                                                       cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCCGCAAGTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAAACGCTTTTCAAAAttcacgt                                             >  1:674061/1‑179 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                               
GATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGT‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGC‑GGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTC‑CGGCGTCAGGATGTTGTAAATTT  >  NC_000913/1824813‑1825252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑