Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,1970AT43.6% 35.4 / 41.0 39intergenic (+166/+417)yjhR/nanSpseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (5/17);  new base T (17/0);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 5.23e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCC  >  NC_000913/4536063‑4536345
                                                                                                                                      |                                                                                                                                                      
taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGTAATTGAaaaaaa                                                                                                                                          >  1:518380/1‑149 (MQ=255)
       cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagcaa                                                                                                                                 <  2:130571/151‑1 (MQ=255)
            tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATTAAAAAAAATgcaaaaaaaaa                                                                                                                              >  1:319066/1‑142 (MQ=255)
                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATg                                                                                                                   <  1:730660/150‑1 (MQ=255)
                             gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGTATTAAAAAAAAAAAAAACAATGATATAATAAAAGTCTAAttt                                                                                                         >  1:56524/1‑151 (MQ=255)
                              aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTc                                                                                                        <  2:468075/151‑1 (MQ=255)
                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTTCTTTTTGGAATGAAAGAAAAATGAAGCCAAAATAAAATAAAAGCCTCAATTTCt                                                                                                       >  2:773992/1‑151 (MQ=255)
                                tttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAGCCAACAATAATATAATAAATGGCTAAtttttt                                                                                                      >  2:86952/1‑151 (MQ=255)
                                    cGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                  <  1:906677/151‑1 (MQ=255)
                                          aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatcat                                                                                            <  2:607894/151‑1 (MQ=255)
                                          aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatcat                                                                                            <  2:10277/151‑1 (MQ=255)
                                              cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATTAAAAAAAAAGCAAACAATAAAATAAAAAAAGGCTAATTTTTTTTAACATCATTaa                                                                                         >  1:155199/1‑150 (MQ=255)
                                              cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCCTTTTTGTATGTAAAAAAAAACAAGAAAAAATATAAAAAAAGGCCAACTTCtttttaca                                                                                                >  2:109908/1‑138 (MQ=255)
                                                      cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAAAAAAAGCAAGCAATAAAATAAAAAAAGGCTAATTTTTTTTTAAATCATTAAGAata                                                                                    >  1:255657/1‑147 (MQ=255)
                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATTAAAAAAAAAACAAGCAATAATATAAAAAAAGCCTAATTTTTTTTTAAATAATTAAACata                                                                                    >  1:822387/1‑143 (MQ=255)
                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                            <  1:905191/151‑1 (MQ=255)
                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                            <  1:13180/151‑1 (MQ=255)
                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcat                                                                             <  1:840260/150‑1 (MQ=255)
                                                                tataAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTa                                                                       <  1:557222/150‑1 (MQ=255)
                                                                      tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGTATGAAAAAAAAAAAAAACAATAAAATAATAATTGCCTAATTTTTTTTTAAATCATTAAGCATATTATCATCGTCTAttttt                                                                  >  1:407012/1‑147 (MQ=255)
                                                                         tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCa                                                                <  1:623179/148‑1 (MQ=255)
                                                                            tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATc                                                           <  1:596818/150‑1 (MQ=255)
                                                                               ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                        <  1:759615/150‑1 (MQ=255)
                                                                                 aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAGCAAGCAATAATATAATAATAGCCTAATTTTTTTTAAAATCATTAAACATAACATCATCGTCTATTTTCCAGAATCCGTc                                                       >  2:397758/1‑149 (MQ=255)
                                                                                           tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGAAAGAAAGAAAATGCAAGCAATAATAAAAAAATAGCCCAATTTTTTTTAACATAATTAAGACTATCATCATTGACTATTTGCCAGAATCCGTCGACCTATaa                                              >  2:684757/1‑148 (MQ=255)
                                                                                           tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATTGAAGAAAAAGCAAAAAATAAAATAAAAAAAGACTAATTTTTTTTTAAAACATTAAGAAAATAATCATCGTCTATTTTTCACAATCCTCCGACtttt                                                >  1:523783/1‑144 (MQ=255)
                                                                                             tCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTGAAGAAAAAACAAGCAATAAAATAATAAATGGCTTATTTTTTTTTAAAACATTTAAAAAAATATCATCGTCTTTTTTTCAAAATTCGTCGACTTTTaa                                              >  1:329407/1‑146 (MQ=255)
                                                                                                  tGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGGAAGAAAAAGGAAGCAAAAATATAATAATAGCCTAATTTTTTTTTAAAACAATAAACAAATTATCATCCTCTATTTTTCAGAATCCGTCGGCTTATAATAGcccaaa                                     >  1:438021/1‑146 (MQ=255)
                                                                                                             ccTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAGAAAAAGCAAGAAAAAAAATAATAAAAGGCTAATTTTTTTTAACAACATTAAACATATCATTATCGTCTATTTTTCAGAAACCGTCGACTTATAAAATCCCGAAAGGCATATCAgg                         >  1:181642/1‑151 (MQ=255)
                                                                                                             ccTGGGATTAATGGCGGGTGGTTTTTTGAAAGAAAAAAAAAGCAAGAAATAATATAAAAAAAGCCTAATTTCTTTTAAAATCATTAAACAAATCATCACCGTCTTTTTTCCAGAATCCCTCGACTTATAAACCCCAAAACGGCGTATc                            >  2:118851/1‑148 (MQ=255)
                                                                                                                gggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGa                        <  2:550177/149‑1 (MQ=255)
                                                                                                                     tAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgc                  <  1:109910/150‑1 (MQ=255)
                                                                                                                        tGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg              >  2:107928/1‑151 (MQ=255)
                                                                                                                            gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg          >  1:434479/1‑151 (MQ=255)
                                                                                                                             ggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGc         <  2:255660/151‑1 (MQ=255)
                                                                                                                               tGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtg       >  2:68285/1‑151 (MQ=255)
                                                                                                                                gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtg       <  2:319072/150‑1 (MQ=255)
                                                                                                                                   tttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGc   >  2:476517/1‑151 (MQ=255)
                                                                                                                                    ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc  >  2:89288/1‑151 (MQ=255)
                                                                                                                                      |                                                                                                                                                      
TAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCC  >  NC_000913/4536063‑4536345

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