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

ACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT  >  NC_000913/1824917‑1825205
                                                                                                                                                |                                                                                                                                                
aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                            <  2:963678/151‑1 (MQ=255)
aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                            >  1:793077/1‑151 (MQ=255)
    atatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                        >  2:805897/1‑151 (MQ=255)
      atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTTTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTGGGGTGATTAACGTTAGTGGTGATCAGGCTTTTTTTACAACTTTTTGAGAAGAAAAAAGAAAACCGCCGa                                                                                                                                       <  2:441228/150‑1 (MQ=255)
         aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                   >  2:303101/1‑151 (MQ=255)
            gTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc                                                                                                                                >  2:61181/1‑151 (MQ=255)
              aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                               >  2:553429/1‑150 (MQ=255)
                 ccatccccccccccGGCCCTTTTTTTTCTTTTTTTTCGTTCGGGTACCCGTGGTTGTTTATCGGCAACTGGGTTGATTTACTTGAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                                                                                                             <  1:144372/138‑1 (MQ=255)
                         gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg                                                                                                                   >  2:292708/1‑151 (MQ=255)
                           tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                  >  2:590493/1‑150 (MQ=255)
                           ccccGCCCCTTTCTTTTCATTTTTATCGTTCAGGTACCGGTGGTTGTTTATCGGCATCTGGGGTGATTTACGTTAGGGGTGATCGGTTTTTTTTTACACTTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACt                                                                                                                    <  2:539306/147‑1 (MQ=255)
                            cccGGCCCTTTTTTTTCATTTTTATCGTTCAGGTACCCGGGGATGGTTATCGACAACTGGGGTAATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                  <  1:99110/149‑1 (MQ=255)
                                     tttttttCTTTTTTCTCTTTCAGGTAACGGTGGTTGTTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTAAAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                         <  1:683025/146‑1 (MQ=255)
                                        ttttTATTGTTTTCGTTCAGGTACCCGGGGATGGTTATCGGCACCTGGGGTTATTTGCGTTAGTGGTGATCAGGTTTTTTTTACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                      <  1:252333/149‑1 (MQ=255)
                                                 ttATCGTTCAGGTAACCGTGGATGGTTATCGGCACCTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                            <  2:241769/150‑1 (MQ=255)
                                                    tCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt                                                                                        >  2:718713/1‑151 (MQ=255)
                                                    tCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt                                                                                        >  1:400879/1‑151 (MQ=255)
                                                       ttCAGGTACCGGTGGTTGTTTATCGGCAACTGGGTTGATTTAGGTTAGTGGTGATCGGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                       <  2:798192/149‑1 (MQ=255)
                                                                  gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                           >  2:515905/1‑150 (MQ=255)
                                                                        ggTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                           >  1:992960/1‑144 (MQ=255)
                                                                              cgacAACTGGGGTGTTTTAGTTGAGGGGTTATGAGGTTTTTTTACCAACTTTTGGAAAAGAGAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                               <  1:133798/147‑1 (MQ=255)
                                                                                   aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                           >  1:259733/1‑149 (MQ=255)
                                                                                       gggTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                     <  1:113840/151‑1 (MQ=255)
                                                                                        ggTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                     >  1:480630/1‑150 (MQ=255)
                                                                                                   ttAGTGGTGATCGGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                           <  2:1126264/149‑1 (MQ=255)
                                                                                                        ggtgATCAGGTTTTTTTTACATTTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                     <  2:445748/150‑1 (MQ=255)
                                                                                                               agtttttttttACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                              <  2:509844/146‑1 (MQ=255)
                                                                                                                   tttttttACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                                  <  1:1090071/142‑1 (MQ=255)
                                                                                                                     tttttAAAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                        <  1:553434/150‑1 (MQ=255)
                                                                                                                      ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                      >  2:390662/1‑151 (MQ=255)
                                                                                                                          aCAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                   <  1:1099591/150‑1 (MQ=255)
                                                                                                                            aaCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGTAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                   >  2:567447/1‑148 (MQ=255)
                                                                                                                                  ttGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt           <  1:901366/150‑1 (MQ=255)
                                                                                                                                   tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct         >  1:285094/1‑151 (MQ=255)
                                                                                                                                     aaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa       <  1:307016/149‑1 (MQ=255)
                                                                                                                                       aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc        <  2:159636/148‑1 (MQ=255)
                                                                                                                                            aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGACGCTTCTCAAAAtt  >  1:820374/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT  >  NC_000913/1824917‑1825205

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