Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT47.1% 41.3 / 24.0 51S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (6/21);  new base T (24/0);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  NC_000913/1604091‑1604386
                                                                                                                                                  |                                                                                                                                                      
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGttt                                                                                                                                                    <  1:641790/151‑1 (MQ=255)
      tGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGttttggt                                                                                                                                                >  2:319584/1‑145 (MQ=255)
       gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGttt                                                                                                                                                    >  2:541197/1‑144 (MQ=255)
         gaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGAcc                                                                                                                                           <  2:849646/151‑1 (MQ=255)
               cGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTTGACCATCgg                                                                                                                                      >  2:655406/1‑150 (MQ=255)
                gCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATCGGcc                                                                                                                                    >  2:392345/1‑151 (MQ=255)
                  aTTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTGGTTCGGCCCTCGGCct                                                                                                                                   >  1:514085/1‑150 (MQ=255)
                  aTTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCctc                                                                                                                                  <  2:341205/151‑1 (MQ=255)
                         gCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCca                                                                                                                           <  1:543125/151‑1 (MQ=255)
                            aTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTTGGCCAACGGCCTCTGCGCcatt                                                                                                                         >  1:791910/1‑150 (MQ=255)
                              tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGAACATCGGCCTCTTCGCcatt                                                                                                                         >  2:173955/1‑148 (MQ=255)
                               aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTGGTTCGACCAACGGCCTCTGCGCcctttcat                                                                                                                    >  1:813679/1‑151 (MQ=255)
                                aCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcatt                                                                                                                   <  2:433355/151‑1 (MQ=255)
                                 cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcattg                                                                                                                  <  2:52085/151‑1 (MQ=255)
                                   cgcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGGCTCTGCGCCATTGCATTggg                                                                                                                >  2:24627/1‑151 (MQ=255)
                                       aaTGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTa                                                                                                            <  2:528402/151‑1 (MQ=255)
                                          gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTg                                                                                                         <  2:122326/151‑1 (MQ=255)
                                                 ctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTCGACCAACGGCCTCCGCGGCCTTGCATTTGGGGGACTTTTTTaa                                                                                                   >  2:755942/1‑150 (MQ=255)
                                                     atatGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtgg                                                                                              <  1:96724/151‑1 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTTGACCATCGGGCCCTGCGCCATTTCATTGGGGGTACTGTTTTTAAgtggt                                                                                             >  2:589552/1‑150 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtggtg                                                                                            <  2:813735/151‑1 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtggtg                                                                                            <  2:746858/151‑1 (MQ=255)
                                                          ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGAAATTTTGTTTTGTTTGCCCATCCGCCTCTGCGCCATTGCATTGGGCGGACTGTTTTCAAgtgg                                                                                              >  2:73770/1‑146 (MQ=255)
                                                            gCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTCGTTTTGGTCGACCCTCCGCCCCTGCGCCATTGCCTTTGGGGGGCTGttttcaaatgggg                                                                                            >  2:571113/1‑136 (MQ=255)
                                                              tGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGCGTAGAAATTTTGTTTGGTTTGACCGTTCCCCCCTGCGCCAATTCTTTGGGCGCGCTGTTTGTAAGAGGTGTTGCGAt                                                                                     >  2:527116/1‑151 (MQ=255)
                                                               gTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATg                                                                                    <  1:773933/151‑1 (MQ=255)
                                                                   aGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCCCTGCGCCATTGCATTTGGCGTACTGTTTTTAAGTGGGGTGGCGATGcc                                                                                  >  1:399500/1‑149 (MQ=255)
                                                                         aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTTTTTGGTTCCACCATCCGCCTCTGCGCCCCTGCATTTGGCGGACTGGTTTTAAGTGGGGGGGCGGTGcccccgggg                                                                           >  1:253957/1‑143 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGAGTCCTTGCATTGGGCGGACTGTTTCAAAGGGGGGTGCCGATTCCGCTGGCg                                                                           >  2:360525/1‑149 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGAt                                                                         <  2:612384/151‑1 (MQ=255)
                                                                              cTTTATCTGCATTGTCATTTTCTCCCTACCTCTAACGATGGTGATTGTCAGTGGCGGGTTCGATATTTTGTTTTGTTTGACCCTCGGCCTCTGCGCCCTTGCCTTTGGCGTACTGGTTCAAAATGGTGTGTCGGTGCCGCTGGTGGTACt                                                                      >  2:848038/1‑150 (MQ=255)
                                                                                 tATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCTTCCATCGGCCCCTGCGCCCTTGCATTTGGGGGACCTGTTTCAAGTGGGGGGCCGGGGCCGCCGGCgggcctcc                                                                    >  1:89837/1‑142 (MQ=255)
                                                                                   tCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTTCACCCTCCGCCCCTGCGCCATTGCATTGGGGGTGCTGGTTTAAAGAGGTGTGCCGCTGCCGCTGGGGAtcatcccg                                                                  >  1:697958/1‑142 (MQ=255)
                                                                                       cattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGCACGGTGTTTCGTTTTGTTCGACCATCAGCCCCCCCGCCATTTCATTGGGGGTACTGTTTTCAAGTGGGGGGCCCATGCCGCTGGCGGTACCCCTGagc                                                               >  2:82935/1‑146 (MQ=255)
                                                                                        attggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGAGAACTTTTCTTTTGTTCGGCCCACCGCCCCTGCGCCCTTGGCTTTGGCGCTCTGTTTTTAAAGGGGGTGCCCGAGCCGCCGGCgg                                                                          >  1:12397/1‑135 (MQ=255)
                                                                                                  tCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCTTTTGGTTTGGCCATCGGCCTCTGGGCCCTTTCATTTGGCGTACTGTTTTTAAGGGGGGGGCCGAAGCCGCTGGCGGTACTCCTGGCCTTActcctccg                                                    >  1:31091/1‑146 (MQ=255)
                                                                                                   cGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCa                                                <  1:630138/151‑1 (MQ=255)
                                                                                                        tACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCCGCGCCACTTCCTTTGGCGTTCTGTTTTAAAGGGGGGTGGCGGTGCCGCCGGCGGAACTCCTGACCTTAATtctccg                                                    >  1:811678/1‑140 (MQ=255)
                                                                                                         aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                           <  1:45799/150‑1 (MQ=255)
                                                                                                               aaCGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGGCCTTGCATTTGGGGTACTGTTTCAAAGGGGGATGCCGATGCCGCTGGCGGTACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGGTg                                    >  1:405929/1‑151 (MQ=255)
                                                                                                                  gATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCCTTTGGCGTACTTTTTTTAAGTGGGGTGCCGCTGCCGCCGGCGGTACTCCTGACCTTTCTTCTCGGCGCCTTGTGCggg                                       >  1:425603/1‑145 (MQ=255)
                                                                                                                  gATGGTGATTGTCAGTGGCGGGATCGATAATTTGTTTGGTTTGCCCATCGGCCTCTGCGCCATTTCATTGGGGGTACTGTTTTAAAAGGGTGGGCCGGTGCCGCCGGCGGTACTCCTGACCCTTCTGCTCCGGGCATTGTGGGGGGTggt                                  >  1:830801/1‑148 (MQ=255)
                                                                                                                        gATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGcc                           <  2:498660/151‑1 (MQ=255)
                                                                                                                        gATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGcc                           <  1:300904/151‑1 (MQ=255)
                                                                                                                        gATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCCGCCTCTGCGCCATTGCATTGGGCGGACTGTTTTAAAGTGGGGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCCACGcc                           >  1:418557/1‑151 (MQ=255)
                                                                                                                          ttGTCAGTGGCGGGATCGATATTTTGTTTGGTTCCACCAACGGCCTCTTCGCCCTTGCATTTGGGGGACTGTTTTAAAGTGGGGTGCCGATGCCGCCGGGGGTACCCCCGACCTTACTGCTCGGGGCCATTTGCGGGCTGGTCAAACCCgg                         >  1:696200/1‑151 (MQ=255)
                                                                                                                            gTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGat                       <  2:54678/151‑1 (MQ=255)
                                                                                                                                    cGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATc               <  1:51217/151‑1 (MQ=255)
                                                                                                                                       gATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCtat            <  2:836028/151‑1 (MQ=255)
                                                                                                                                             tatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa       >  1:779827/1‑150 (MQ=255)
                                                                                                                                                 tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt   <  1:653431/150‑1 (MQ=255)
                                                                                                                                                  cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTa  <  1:361264/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  NC_000913/1604091‑1604386

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