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

AAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGT  >  NC_000913/4536064‑4536341
                                                                                                                                     |                                                                                                                                                
aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATg                                                                                                                                  <  1:497496/150‑1 (MQ=255)
    ggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagc                                                                                                                             <  1:240957/151‑1 (MQ=255)
       atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAacaagaa                                                                                                                            >  1:469574/1‑147 (MQ=255)
        tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATGTAaaaaaaaaga                                                                                                                                >  1:604659/1‑140 (MQ=255)
         atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa                                                                                                                        <  1:308183/151‑1 (MQ=255)
         atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa                                                                                                                        <  1:434540/151‑1 (MQ=255)
            ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAAATAAAAAAAAAGAAAACaaaaa                                                                                                                        >  2:246925/1‑148 (MQ=255)
                cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAAAAAGCAaaaaaataa                                                                                                                   >  1:98407/1‑149 (MQ=255)
                  aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGATTTTTTTACTGAAAAAAAATACCAGCAAAAAAAtaaaaa                                                                                                                >  2:102839/1‑150 (MQ=255)
                        aGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATGAAAAAAAAAAAAAGCAAAAAAAtaataaatacc                                                                                                           >  2:587133/1‑144 (MQ=255)
                             aCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATTAAAAAAAAAGCAAGAAAAAATATAAAAAATGGCTAAttt                                                                                                     >  2:535280/1‑150 (MQ=255)
                                  ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTg                                                                                                <  2:279802/150‑1 (MQ=255)
                                       ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcat                                                                                           <  1:329436/150‑1 (MQ=255)
                                          ttGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATTAAAAAAAAAGCAAAGAAAAATAAAATAAAGGCCTAATTTTTTTTAacaaca                                                                                         >  2:432729/1‑149 (MQ=255)
                                           tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCAtt                                                                                       <  2:477594/150‑1 (MQ=255)
                                            gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGTATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTa                                                                                      <  2:141455/150‑1 (MQ=255)
                                              gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTa                                                                                      <  2:361592/148‑1 (MQ=255)
                                              gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGc                                                                                   <  2:610314/151‑1 (MQ=255)
                                                         ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcatc                                                                        <  1:527377/151‑1 (MQ=255)
                                                            gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTc                                                                     <  2:584647/151‑1 (MQ=255)
                                                                  aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTg                                                               <  2:466493/151‑1 (MQ=255)
                                                                       ttATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAGAAAAAACAAACAAAAATATAATAAATGCCTAATTTTTTTTTAAATCATTTAACAAATCCTCATCGTCTttttttc                                                             >  1:469192/1‑145 (MQ=255)
                                                                            gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAGAAAAAGCAAACAAAAATATAATAATAGGCTAATTTTTTTTTAAATCATTAAACATATCATCATCGTCTAtttttccgaa                                                         >  1:85251/1‑140 (MQ=255)
                                                                              ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGCTTTTTTGATTGAAAGAAAAAGCAAGCAATAATATAATAAATGACTAATTTTTTTGTAAATTATTTAGCATATCATCATCGTCTTTTTTCCCGAATcc                                                      >  2:31094/1‑148 (MQ=255)
                                                                                   tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTAGGATTAAAGAAAAAACCAACAAAAAAATAATAATTGCCCAATTTTTTTTTAAATCATTAAGAATATCATCATCTTCTATTTTCCAGAAAcc                                                      >  1:550055/1‑143 (MQ=255)
                                                                                      tCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAGAAAATGCAAACAAAAATATAATAAAAGGCTAATTTTTTTTTAAAATATTAAAAATATCATCATCTTCTATTTTTCAGAATCCGTCGCCtt                                              >  1:394314/1‑148 (MQ=255)
                                                                                       cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                          <  1:76281/151‑1 (MQ=255)
                                                                                       cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                          <  1:121040/151‑1 (MQ=255)
                                                                                        aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGATTTTTTTTTTTATAAAAAAAGCCAGAAAAAAAATAAAAAATGCGTAATTTTTTTTAACATCATTAAAAATATCATAATCTACTATTTTTCAGAAAACGCCGACTttta                                           >  2:131885/1‑149 (MQ=255)
                                                                                             ccTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATTAAAAAAAAATAGAGAAAAAATATAATAAATGACTATTTTTTTTTTACATCATTAAACAAATAACCATACTCTATTTTTCAGAATTCCTCCAAGTATAAtttcc                                     >  2:344087/1‑146 (MQ=255)
                                                                                                ttGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGa                                  <  2:2746/150‑1 (MQ=255)
                                                                                                   cATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAGCAAGAAAAAAAATAAAAAAAGCCTAATTTTTTTTTAAAACATTAAACATATCATCATCGTCTATTTTTCAGAAACCGTCGACTTATAATTGCCAGaa                                 >  1:399291/1‑148 (MQ=255)
                                                                                                               gggATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACg                  <  2:174801/151‑1 (MQ=255)
                                                                                                                     aaTGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgct             <  1:54371/150‑1 (MQ=255)
                                                                                                                           gggTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACg      >  1:194749/1‑151 (MQ=255)
                                                                                                                               gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgt  >  1:281025/1‑151 (MQ=255)
                                                                                                                                     |                                                                                                                                                
AAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGT  >  NC_000913/4536064‑4536341

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