Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA42.1% 43.0 / 15.7 38intergenic (‑124/+79)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (15/7);  new base A (0/16);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 1.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824913‑1825208
                                                                                                                                                     |                                                                                                                                                   
ttttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                    >  1:592483/1‑150 (MQ=255)
 tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                                   >  2:777332/1‑150 (MQ=255)
       cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                            >  2:325955/1‑151 (MQ=255)
         tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                           >  1:432835/1‑150 (MQ=255)
              aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                                     >  1:843893/1‑151 (MQ=255)
                 tATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                                  >  2:652966/1‑151 (MQ=255)
                  aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                                  >  1:476216/1‑150 (MQ=255)
                   ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                                >  2:80399/1‑151 (MQ=255)
                     ccATACTCGCTCCCGGCCCGTTTTTTTCATTGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCg                                                                                                                              <  1:855471/151‑1 (MQ=255)
                          cTCGCCCCCGGCCCGTTTTTTTCATTTTTATCGTTCAGTTACCCGTGGATGTTTATCGGCAACTCGGGTGTTTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCa                                                                                                                            <  1:624888/148‑1 (MQ=255)
                              cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                     >  1:656186/1‑151 (MQ=255)
                                cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                    >  1:57992/1‑150 (MQ=255)
                                         tttttttCATTTTTATCTTTCGGGTAACCGTGGTTGTTTATCGCCAACTCGGGTGTTTTCCTTTAGTGGTGATCGGTTTTTTTTCACAACTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                               <  2:720683/143‑1 (MQ=255)
                                         ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                           >  2:11097/1‑150 (MQ=255)
                                             tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                        >  1:340074/1‑149 (MQ=255)
                                                     ttATCGTTCAGGTACCGGTGGTTGGTTACCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGTTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                <  1:383133/149‑1 (MQ=255)
                                                      tATCGTTCGGGTAACGGGGGATGTTTATCGGCACCTCGGGTGATTTACGTTAGGGGTGATCGGGTTTTTTTTACAATTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                              <  2:675441/150‑1 (MQ=255)
                                                           ttaaGGTAACCGTGGATGTTTATCGGCAACTGGGGTGTTTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                             <  1:812212/143‑1 (MQ=255)
                                                           ttCAGTTACCCTGGGGTGTTTATCGGCAACTGGGTTGATTTAGTTTAGTGGTGATCAGGTTTTTTTAACAACTTTTTGAAAAAAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                            <  2:101706/147‑1 (MQ=255)
                                                           ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt                                                                                           <  2:274255/148‑1 (MQ=255)
                                                                    ccGTGGATGGTTATCGGCAACTGGGGTGTTTTCCGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                               <  1:160744/151‑1 (MQ=255)
                                                                      gTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                                                                             <  1:580677/151‑1 (MQ=255)
                                                                      gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                               >  1:547729/1‑149 (MQ=255)
                                                                               tATCGGAACTTCGGGGGTTTTACTTTAGTGGTGACCGGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                     <  2:609190/150‑1 (MQ=255)
                                                                                 acagacaactgggTTATTTACGGTTGTTGGTGATCGGGTTTTTTTAACATCTGTTGGGAAAGAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                     <  2:337915/139‑1 (MQ=255)
                                                                                      aaCTGGGGTATTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                               <  2:174638/149‑1 (MQ=255)
                                                                                      aaCTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                                 <  1:313578/147‑1 (MQ=255)
                                                                                      aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                              >  2:630623/1‑150 (MQ=255)
                                                                                         tggggTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                              <  1:487056/145‑1 (MQ=255)
                                                                                           gggTGATTTACGTTAGTGGTGATCGGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt                                                           <  2:27532/148‑1 (MQ=255)
                                                                                           gggTGATTTACGTTAGGGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag                                                            <  1:204408/147‑1 (MQ=255)
                                                                                            ggTGTTTTCCGTTAGGGGTGATCATGTTTTTTTCACACCTGTTTGAGAAAAAAAATGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                                       <  1:162472/151‑1 (MQ=255)
                                                                                                     cGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt                                               >  1:871212/1‑150 (MQ=255)
                                                                                                           gtgttgTTCGGGTTTTTTTTACAACTGTTTGAAAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                         <  2:325547/151‑1 (MQ=255)
                                                                                                           gtggtgATCAGTTTTTTTTTACACTTTTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc                                            <  2:790945/148‑1 (MQ=255)
                                                                                                             ggtgAGCAGTTTTTTTTTACACCTGTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                              <  2:886304/144‑1 (MQ=255)
                                                                                                                     ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                >  1:663283/1‑150 (MQ=255)
                                                                                                                           ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                         >  1:138521/1‑151 (MQ=255)
                                                                                                                            tttACAACTTTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                         <  2:277419/150‑1 (MQ=255)
                                                                                                                                      ttggaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                  <  2:651339/141‑1 (MQ=255)
                                                                                                                                               agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     >  1:752616/1‑151 (MQ=255)
                                                                                                                                               agaaaaaaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAt      <  2:705018/150‑1 (MQ=255)
                                                                                                                                               aaaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  1:700220/147‑1 (MQ=255)
                                                                                                                                                 aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTTCCGCTGGACGGTCTGTCCGACGCTTCTCAAAAtt     >  2:852460/1‑149 (MQ=255)
                                                                                                                                                  aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg   >  1:755604/1‑150 (MQ=255)
                                                                                                                                                  aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  <  2:824911/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824913‑1825208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: