Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,4100TG41.8% 92.6 / 11.7 56intergenic (+112/‑35)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base T (4/28);  new base G (23/0);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 9.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT  >  NC_000913/15261‑15536
                                                                                                                                                        |                                                                                                                               
aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt                                                                                                                                <  1:759228/151‑1 (MQ=255)
  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg                                                                                                                               <  2:180421/150‑1 (MQ=255)
   gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgtg                                                                                                                             <  2:754998/151‑1 (MQ=255)
       gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCTAAACGCCTACTTTAGGGGTGTgg                                                                                                                            >  2:171165/1‑148 (MQ=255)
        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                                                                                                        <  2:1077734/151‑1 (MQ=255)
           aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                                                                                     <  2:580354/151‑1 (MQ=255)
            agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                                                    <  1:490325/151‑1 (MQ=255)
             gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                                                                                                   <  1:663325/151‑1 (MQ=255)
                  ttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGCTTACCCCCATAAACGCCGACCTAAAGGGTGTGGTATTACGcc                                                                                                                   >  1:1115647/1‑146 (MQ=255)
                   tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCCAAGCGCTAAATTAAAGGGTGTGGGATTATGCCTggt                                                                                                               >  1:566365/1‑147 (MQ=255)
                    gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                                                                                                            <  2:794083/151‑1 (MQ=255)
                      cgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTAATCCCCTAAGCGCCCACCTTAGGGGTGGGGGAGTTCGACCGATATg                                                                                                            >  2:898420/1‑149 (MQ=255)
                        acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTAACTTAAAGGGTGTGGGATTTCGCCTGATATGAttt                                                                                                        >  1:1118499/1‑151 (MQ=255)
                         ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                       <  1:384032/151‑1 (MQ=255)
                             aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAGAGCTCAATTAAAGGGTGGGGGGTTTCGCCTGCTATGAatt                                                                                                        >  1:182743/1‑146 (MQ=255)
                             aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATAAGATTTAAc                                                                                                     <  2:699596/149‑1 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCCTAAGCGCGAAATTTAGGGGTGTGGGTTTTCGCCCGGTATTAttt                                                                                                        >  2:990812/1‑145 (MQ=255)
                              cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                                                                                                   <  1:990839/150‑1 (MQ=255)
                                cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGca                                                                                                <  2:1000454/151‑2 (MQ=255)
                                cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGca                                                                                                <  2:1000456/151‑2 (MQ=255)
                                  cTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCCAAACGCTGACTTTAAGGGTGTGGGGTTATCCCTGGTATGtttttaaatgcc                                                                                                >  2:1026181/1‑141 (MQ=255)
                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATAAATCCCATAAACGCTAACTTAAAGGGTGTGGTATTACGCCCGAAATGATTTTAAGTGcc                                                                                                >  2:561837/1‑147 (MQ=255)
                                     aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa                                                                                           <  1:373740/151‑1 (MQ=255)
                                      ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                          <  1:807768/151‑1 (MQ=255)
                                      ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCGAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                          <  2:654213/151‑1 (MQ=255)
                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                       <  1:752302/151‑1 (MQ=255)
                                          cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                                      <  1:840106/151‑1 (MQ=255)
                                              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAAAGCTGAAATAAAGGGTGTGGTGTTAAGCCTGATATGGTTTTAAACGCCGATGaa                                                                                          >  1:506212/1‑143 (MQ=255)
                                              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACCCCCCCAAACGCTCAATTAAGGGGTGTGGGATTACGCCTGATATTTTTTAAACTGcc                                                                                                >  1:1072006/1‑137 (MQ=255)
                                              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                                  <  2:408385/151‑1 (MQ=255)
                                                     gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                                                            <  2:266980/150‑1 (MQ=255)
                                                        cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTCACTTAAGGGTTGTGGTGTTACGCCTGATATGATTTAACGTGCCGATGAATTAATCTCACGATaa                                                                           >  1:837599/1‑148 (MQ=25)
                                                         gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCAAAAAGCCTAACTTTAAGGGTGTGGGATTAAGCCTGCTATGATTTAACAAGTCGAAGAATTTCTCTCCCCATAACTg                                                                        >  1:1073297/1‑150 (MQ=255)
                                                           gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                                     <  1:614548/151‑1 (MQ=255)
                                                              cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAAAGCTAAATTTAAGGGTGTGGGATTTAGCCTGGTATGATTTTAAGTGCCGATGAATTAATCTCCCGATAACTgg                                                                       >  2:463900/1‑146 (MQ=255)
                                                              cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGACTCCCCTAAAGGCTAACTTCAAGGGTTTGGTATTTCCCCCGCTTTGATTTTACGTGCCGGAGGATTAATCTCAccaaaaatagggcg                                                                   >  1:14890/1‑137 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAACGCCAACATAAAGGGTGTGGGGGTATGCCCGGTTTTATTTAACCTGCCGGAGAAATAATCTCACGATAACAGGGCa                                                                    >  1:617618/1‑147 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                                <  1:540027/151‑1 (MQ=255)
                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATAAACCCCCTAAACGCTAACTTTAAGGGTGTGGTGTTTCACCCGGTTTTATTTAATCTGCCGATGGATTACTCTCCCGCAAAATGGTCAGCa                                                                 >  1:565225/1‑150 (MQ=255)
                                                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCATAAAAGCTCAATTAAGGGGTGTGGGGTTTCGCCTGATTTGATTTTAAGGGCCGCTGAATTAATTTCACGaaaaataggccacaaattt                                                            >  2:579739/1‑135 (MQ=255)
                                                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCCAAAAGCTAACTTAAAGGGTGTGGGGTTACGCCTGTTTTTATTTAACGTGCCGATGAAATACTCTCACCATAAATTGTCAGCAATTCt                                                            >  2:707816/1‑150 (MQ=255)
                                                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                                          <  2:672461/151‑1 (MQ=255)
                                                                          aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCTAACTTAAGGGGTGTGGGATTACGCCTGATATGATTTTACGTGCCGATGAATTACTCTCACGATAAATGGTCAGCAAATCTGGccc                                                       >  1:1174788/1‑150 (MQ=21)
                                                                            aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAAAGCCCACTTAAGGGGTGGGGGATTACGCCCGCTATGATTTAACGTGCCGGTGAATTACTCTCACCATAACTGGTCAGGAATTCTGGccc                                                       >  2:416263/1‑148 (MQ=255)
                                                                               aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTAACTCCCATAAACGCTAAATTTAGGGGTGTGGGGTTACGCCTGATATGATTTAAAGTGCCGACGGATTAATTTCACAAAAACATGGCAGACATTCTTGGCCACAtt                                                  >  2:835729/1‑150 (MQ=255)
                                                                                     cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCCCTCCCCTAAACGCCAACTTTAGGGGTGTGGGATTTCGCCCGTTATGATTTTAAATGCCGATGAATTACTCTCACCATAAATGGTCAGCAAATCTTGGCCCTAtt                                                  >  1:677184/1‑144 (MQ=255)
                                                                                        tGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGAGACTCCCCTAAGGGGCAACTTAAGGGGTGTGGGGTTACGCCTGATATGATTTTACATTGCGGTGAAATACTCTCCAGAGAAAATGTCAGCAAATCTTGGCCACATTGGTAAGc                                          >  2:913716/1‑148 (MQ=255)
                                                                                            gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaaga                                   <  1:480022/151‑1 (MQ=255)
                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCCCAAACGCCTACCTAAAGGGTTTGGGATTTCTCCCGGTATGATTTTAAGTTCCGATGAATTTCTCTCCCCATAAATGGTCCGCAATTCTTGCCCCCATTTGTAAGccc                                        >  2:621493/1‑145 (MQ=255)
                                                                                               atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                                <  1:368365/151‑1 (MQ=255)
                                                                                                 gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                              <  1:407607/151‑1 (MQ=255)
                                                                                                     gCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAACGCCAAAATAAGGGGTGGGGGGTTAAGCCTGATATGATTTAACGTGCCGCTGGAATAATCTCACGATAACTGGTTCGCAATTTTGGCCCCCATTTGTAAAAACGAAGAAAAGGATa                           >  2:627366/1‑150 (MQ=9)
                                                                                                       aCCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACtt                        <  2:1161389/151‑1 (MQ=255)
                                                                                                         cTGTAAAGTGGCGGGGGTCACCCCCCTAAACCCCCAATTTAAGGGTGTGGGATTACGCCTGATATGATTTAAAGTGCCCCTGAATTAATCTTCCGATAAATGGGCAGCCAATTTGGGCCACATTTGTAAAACCCGAAAAACgggtaatt                         >  2:386378/1‑144 (MQ=2)
                                                                                                                   gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg            <  2:138699/151‑1 (MQ=34)
                                                                                                                         aTCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg      <  2:206196/151‑1 (MQ=33)
                                                                                                                             cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctct  >  1:1064235/1‑151 (MQ=11)
                                                                                                                                                        |                                                                                                                               
ATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT  >  NC_000913/15261‑15536

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