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

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  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)
                               aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTGGTTCGACCAACGGCCTCTGCGCcctttca‑t                                                                                                                    >  1:813679/1‑151 (MQ=255)
                                aCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgca‑tt                                                                                                                   <  2:433355/151‑1 (MQ=255)
                                 cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgca‑ttg                                                                                                                  <  2:52085/151‑1 (MQ=255)
                                   cgcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGGCTCTGCGCCATTGCA‑TTggg                                                                                                                >  2:24627/1‑151 (MQ=255)
                                       aaTGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTa                                                                                                            <  2:528402/151‑1 (MQ=255)
                                          gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTg                                                                                                         <  2:122326/151‑1 (MQ=255)
                                                 ctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTCGACCAACGGCCTCCGCGGCCTTGCA‑TTTGGGGGACTTTTTTaa                                                                                                   >  2:755942/1‑150 (MQ=255)
                                                     atatGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAgtgg                                                                                              <  1:96724/151‑1 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTTGACCATCGGGCCCTGCGCCATTTCA‑TTGGGGGTACTGTTTTTAAgtggt                                                                                             >  2:589552/1‑150 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAgtggtg                                                                                            <  2:813735/151‑1 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAgtggtg                                                                                            <  2:746858/151‑1 (MQ=255)
                                                          ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGAAATTTTGTTTTGTTTGCCCATCCGCCTCTGCGCCATTGCA‑TTGGGCGGACTGTTTTCAAgtgg                                                                                              >  2:73770/1‑146 (MQ=255)
                                                            gCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTCGTTTTGGTCGACCCTCCGCCCCTGCGCCATTGCC‑TTTGGGGGGCTGttttcaaatgggg                                                                                            >  2:571113/1‑136 (MQ=255)
                                                              tGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGCGTAGAAATTTTGTTTGGTTTGACCGTTCCCCCCTGCGCCAATTCT‑TTGGGCGCGCTGTTTGTAAGAGGTGTTGCGAt                                                                                     >  2:527116/1‑151 (MQ=255)
                                                               gTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATg                                                                                    <  1:773933/151‑1 (MQ=255)
                                                                   aGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCCCTGCGCCATTGCA‑TTTGGCGTACTGTTTTTAAGTGGGGTGGCGATGcc                                                                                  >  1:399500/1‑149 (MQ=255)
                                                                         aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTTTTTGGTTCCACCATCCGCCTCTGCGCCCCTGCA‑TTTGGCGGACTGGTTTTAAGTGGGGGGGCGGTGcccccgggg                                                                           >  1:253957/1‑143 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGAGTCCTTGCA‑TTGGGCGGACTGTTTCAAAGGGGGGTGCCGATTCCGCTGGCg                                                                           >  2:360525/1‑149 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGAt                                                                         <  2:612384/151‑1 (MQ=255)
                                                                              cTTTATCTGCATTGTCATTTTCTCCCTACCTCTAACGATGGTGATTGTCAGTGGCGGGTTCGATATTTTGTTTTGTTTGACCCTCGGCCTCTGCGCCCTTGCC‑TTTGGCGTACTGGTTCAAAATGGTGTGTCGGTGCCGCTGGTGGTACt                                                                      >  2:848038/1‑150 (MQ=255)
                                                                                 tATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCTTCCATCGGCCCCTGCGCCCTTGCATTTGGGGGACCTGTTTC‑AAGTGGGGGGCCGGGGCCGCCGGCgggcctcc                                                                    >  1:89837/1‑142 (MQ=255)
                                                                                   tCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTTCACCCTCCGCCCCTGCGCCATTGCA‑TTGGGGGTGCTGGTTTAAAGAGGTGTGCCGCTGCCGCTGGGGAtcatcccg                                                                  >  1:697958/1‑142 (MQ=255)
                                                                                       cattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGCACGGTGTTTCGTTTTGTTCGACCATCAGCCCCCCCGCCATTTCA‑TTGGGGGTACTGTTTTCAAGTGGGGGGCCCATGCCGCTGGCGGTACCCCTGagc                                                               >  2:82935/1‑146 (MQ=255)
                                                                                        attggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGAGAACTTTTCTTTTGTTCGGCCCACCGCCCCTGCGCCCTTGGC‑TTTGGCGCTCTGTTTTTAAAGGGGGTGCCCGAGCCGCCGGCgg                                                                          >  1:12397/1‑135 (MQ=255)
                                                                                                  tCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCTTTTGGTTTGGCCATCGGCCTCTGGGCCCTTTCA‑TTTGGCGTACTGTTTTTAAGGGGGGGGCCGAAGCCGCTGGCGGTACTCCTGGCCTTActcctccg                                                    >  1:31091/1‑146 (MQ=255)
                                                                                                   cGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCa                                                <  1:630138/151‑1 (MQ=255)
                                                                                                        tACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCCGCGCCACTTCC‑TTTGGCGTTCTGTTTTAAAGGGGGGTGGCGGTGCCGCCGGCGGAACTCCTGACCTTAATtctccg                                                    >  1:811678/1‑140 (MQ=255)
                                                                                                         aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                           <  1:45799/150‑1 (MQ=255)
                                                                                                               aaCGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGGCCTTGCA‑TTTGGGGTACTGTTTCAAAGGGGGATGCCGATGCCGCTGGCGGTACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGGTg                                    >  1:405929/1‑151 (MQ=255)
                                                                                                                  gATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCC‑TTTGGCGTACTTTTTTTAAGTGGGGTGCCGCTGCCGCCGGCGGTACTCCTGACCTTTCTTCTCGGCGCCTTGTGCggg                                       >  1:425603/1‑145 (MQ=255)
                                                                                                                  gATGGTGATTGTCAGTGGCGGGATCGATAATTTGTTTGGTTTGCCCATCGGCCTCTGCGCCATTTCA‑TTGGGGGTACTGTTTTAAAAGGGTGGGCCGGTGCCGCCGGCGGTACTCCTGACCCTTCTGCTCCGGGCATTGTGGGGGGTggt                                  >  1:830801/1‑148 (MQ=255)
                                                                                                                        gATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGcc                           <  2:498660/151‑1 (MQ=255)
                                                                                                                        gATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGcc                           <  1:300904/151‑1 (MQ=255)
                                                                                                                        gATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCCGCCTCTGCGCCATTGCA‑TTGGGCGGACTGTTTTAAAGTGGGGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCCACGcc                           >  1:418557/1‑151 (MQ=255)
                                                                                                                          ttGTCAGTGGCGGGATCGATATTTTGTTTGGTTCCACCAACGGCCTCTTCGCCCTTGCA‑TTTGGGGGACTGTTTTAAAGTGGGGTGCCGATGCCGCCGGGGGTACCCCCGACCTTACTGCTCGGGGCCATTTGCGGGCTGGTCAAACCCgg                         >  1:696200/1‑151 (MQ=255)
                                                                                                                            gTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGat                       <  2:54678/151‑1 (MQ=255)
                                                                                                                                    cGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATc               <  1:51217/151‑1 (MQ=255)
                                                                                                                                       gATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCtat            <  2:836028/151‑1 (MQ=255)
                                                                                                                                             tatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa       >  1:779827/1‑150 (MQ=255)
                                                                                                                                                 tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt   <  1:653431/150‑1 (MQ=255)
                                                                                                                                                  cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTa  <  1:361264/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCA‑TTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  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: ‑