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

GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536050‑4536347
                                                                                                                                                    |                                                                                                                                                         
gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAgg                                                                                                                                                         <  2:317742/150‑1 (MQ=255)
     aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGaa                                                                                                                                                       >  1:719738/1‑147 (MQ=255)
       gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaag                                                                                                                                                 <  2:9338/151‑1 (MQ=255)
        agtagtAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaaga                                                                                                                                                <  2:135715/151‑1 (MQ=255)
           agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGaaaa                                                                                                                                             <  1:436145/151‑1 (MQ=255)
            gtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                                                                                                                                            <  1:55719/151‑1 (MQ=255)
                       atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa                                                                                                                                 <  2:856067/151‑1 (MQ=255)
                       atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAATACAAACaata                                                                                                                                  >  1:680119/1‑150 (MQ=255)
                          ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                               <  1:295011/150‑1 (MQ=255)
                            aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAGAAAAAGAAAGCAaaaatata                                                                                                                            >  2:203993/1‑150 (MQ=255)
                                  gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATGAAAGAAAAAGGAAGCAATAAAAtaaaaa                                                                                                                       >  2:434518/1‑148 (MQ=255)
                                          gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAGAAAAAGCAAACAATAATATAATAATGGCCTAAttt                                                                                                           >  1:911166/1‑151 (MQ=255)
                                              ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGGATGAAAAAAAAAGCAAGCAATAATATAATAATTGGCTAATtttttt                                                                                                       >  1:46860/1‑151 (MQ=255)
                                                ccGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGa                                                                                                     <  2:357289/151‑1 (MQ=255)
                                                    tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAca                                                                                                  <  2:226463/150‑1 (MQ=255)
                                                      tATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTAAAGAAAAAGCAAGCAAAAATATAATAATGGAGTTATTTTTTTTAAca                                                                                                  >  2:195318/1‑148 (MQ=255)
                                                       aTTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAcatca                                                                                               <  2:400953/150‑1 (MQ=255)
                                                           cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAg                                                                                          <  2:825846/151‑1 (MQ=255)
                                                            gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTAAAAAAAAATCAAACAAAAAAATAAAAATAGGCTTATTTTTTTTAAAAAAATTaa                                                                                           >  2:603553/1‑149 (MQ=255)
                                                               tCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCGTTTTTGAATTGAAAAAAATACACGCAATAATATAATAATAGCCCCATTTTTTTTTAAAACATTaa                                                                                           >  2:291722/1‑146 (MQ=255)
                                                                cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAATAAAACAATAATAAAATAAAAGGCTTAATTTTTTTTAAaataataagaaaa                                                                                       >  2:2608/1‑143 (MQ=255)
                                                                cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCata                                                                                      <  1:337644/150‑1 (MQ=255)
                                                                   cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTGAAGAAAAAGCAAACAATAATATAAAAAATGGCTTAtttttttttacataaataaaca                                                                                         >  2:970950/1‑131 (MQ=255)
                                                                   cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATTTAAGAAAATGCAAACAATAAAATAATAATTGCCTAATTTTTTTGAAAATCATTAAACatataa                                                                                   >  2:323920/1‑148 (MQ=255)
                                                                                 aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAAAAAAAAAAAGCAATAATATAAAAAAGGCCTAATTTTTTTTAACATCATTAAACATATCATCATCGTCTAttt                                                                      >  1:4514/1‑149 (MQ=255)
                                                                                   tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATTAAAAAAAAATTAAGGAATTATATAATAATGGGCTTATTTTTTTTTAAATCATTTAACATATACTTATCTTCTAttt                                                                      >  2:230935/1‑147 (MQ=255)
                                                                                    cTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCGTTTTTGTATTAAAAAAAAAGCAAGCAATAATATAATAATGGGCTTATTTTTTTTTAAAATCTTTAGAATATCATAATTGTATAtttttc                                                                   >  2:325450/1‑146 (MQ=255)
                                                                                      tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAg                                                                 <  1:323983/149‑1 (MQ=255)
                                                                                       aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGaa                                                               <  2:231270/150‑1 (MQ=255)
                                                                                         tGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTGGGATGAAAAAAAAAGCAAACAAAAAAAAAAAAAAAGCCTAATTTTTTTTTAAATCATTAAACAAATCATCATCGTCTATTTTCCAGAATcc                                                            >  1:328229/1‑151 (MQ=255)
                                                                                          gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTTATGGCGGGGGGTTTTTTGAATTAAAAAAAAAACAAGCAATAATATAATAATGGACTAATTTTTTTGAACATAATTAAAAATATCATCATATTCTAttttatcgaa                                                               >  1:273916/1‑140 (MQ=255)
                                                                                              aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTTTGAATTAAAAAAAAAACAAGCAAAAAAATAATAATGGCCTAATTTTTTTTAAAATCATTAAAAATATCATCATCTTCTATTTTTCAGAATCCGt                                                          >  2:548702/1‑148 (MQ=255)
                                                                                              aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGTATGAAAGAAAAAGCAAACCAAAAAATAATAATAGGCTAATTTCTTTGTAAATCATTAAACATATTATCCACGTATATTTTCCAGAAACCGt                                                          >  1:646591/1‑148 (MQ=255)
                                                                                              aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCg                                                        <  1:775482/150‑1 (MQ=255)
                                                                                                   aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtat                                                  <  1:33522/151‑1 (MQ=255)
                                                                                                     cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtata                                                 <  1:937447/150‑1 (MQ=255)
                                                                                                     cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATaa                                                <  1:631885/151‑1 (MQ=255)
                                                                                                      aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTTAATTAAAAAAAAAGCAAGCAAAAATATAATAATTGGCTAATTTTTTTTTAAATTCTTTAGCAAAACATCATCGTCTATTTTCCAGAATCCGTCGACCtatta                                                >  2:679133/1‑148 (MQ=255)
                                                                                                        tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATc                                              <  1:602852/150‑1 (MQ=255)
                                                                                                        tttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATc                                              <  2:611002/150‑1 (MQ=255)
                                                                                                         ttCCTTTGCGTAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATTGAAGAAAAAGCAAAGAATAAAATAATAAATGCCTCAATTCTTTTAAAAACATTTAACAAATACTAATCGTCTATTTTCCAGAAACCGTCGACTTATAATAGc                                            >  2:678066/1‑151 (MQ=255)
                                                                                                            cTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGGAAAAAAAAGCAAAGAATAAAATAATAAAAGGCTTATTTTTTTTTAAAAAAATAAACATATAATTATTGTCTATTTTTCCCAATTCGTCGaacttttatatcc                                           >  2:204659/1‑137 (MQ=255)
                                                                                                            cTTTGCATAATGCGCCTGGGATTAATGGCGGGGGTCTTTTGGGGTTAAAGAAAAAGCAAGCAATAAAATAAAAAAGGGCTTAATTTTTTTTAAATAATTAAACATATAATAAACGTCTATTTTTCAAAAATCGCCGACTTATAATTGCCa                                          >  2:230983/1‑150 (MQ=255)
                                                                                                             tttGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTAGGAAGAAAGAAAAAACAAGCAAAAAAATAATAATAGCCTAATTTCTTTGAACATCCTTAAACATATCATCATCCTCTATTTTCCCGAATCCGTCGACTTATTATTGcc                                           >  1:417226/1‑148 (MQ=255)
                                                                                                              ttGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGGAAGAAAATGCAAAGAAAAAAATAAAAAAGGCCTAATTTTTTTTTAAAACATTAAACAAATCATCATCCTCTATTTTTCAGAATACCTCGGCCTATAATTGcc                                           >  1:587429/1‑147 (MQ=255)
                                                                                                                        cgcCTGGGATTAATGGCGGGGGGTTTTTGGAAAGAAAGAAAAAACAAGCAATAAAATAAAAAAAGACTAAATTTTTTTTAAAACATTAAAAAAATCATTATAGTCTATTTTCCAGAAAACCTCGACTTATAAATACCAGAAAGGCata                                >  2:765121/1‑148 (MQ=255)
                                                                                                                               gATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCACCAGAACGGCATATCAGGACGGt                      <  2:680045/151‑1 (MQ=255)
                                                                                                                                        cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc             >  1:477564/1‑151 (MQ=255)
                                                                                                                                              gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgt        <  2:39417/151‑1 (MQ=255)
                                                                                                                                              gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgt        >  2:459438/1‑151 (MQ=255)
                                                                                                                                                 tttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgtg       <  2:3206/149‑1 (MQ=255)
                                                                                                                                                   tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc   >  1:451793/1‑151 (MQ=255)
                                                                                                                                                   tAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc    <  1:491726/150‑1 (MQ=255)
                                                                                                                                                    aGGAATGAAAGAAAATGCAAGCAATACTATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGc     >  2:565463/1‑148 (MQ=255)
                                                                                                                                                    aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCa  <  1:87315/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                         
GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536050‑4536347

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