Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT42.5% 49.7 / 20.8 47S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (9/18);  new base T (20/0);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 1.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387
                                                                                                                                                   |                                                                                                                                                      
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGttt                                                                                                                                                    <  2:1174539/151‑1 (MQ=255)
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGt                                                                                                                                                      <  2:1138932/149‑1 (MQ=255)
   tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGttt                                                                                                                                                    >  2:346471/1‑148 (MQ=255)
    aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGtt                                                                                                                                                <  2:207430/151‑1 (MQ=255)
         gaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGGTGGTGATTGTCAGTGGCGGGATCAGCATTTTGTTTTGTTCTAcc                                                                                                                                           >  1:871472/1‑151 (MQ=255)
            tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGACATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGAccccc                                                                                                                                        >  2:551607/1‑148 (MQ=255)
               cGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCgg                                                                                                                                      <  2:21642/150‑1 (MQ=255)
                gCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGcc                                                                                                                                    <  2:364050/151‑1 (MQ=255)
                 cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTTGACCCTCCGcc                                                                                                                                    >  1:57389/1‑150 (MQ=255)
                    ttGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTCAATTTTTGTTTTGGTCGAACATccggccct                                                                                                                                 >  2:509022/1‑143 (MQ=255)
                          cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTTTTTTGTTCGACCATCCGCCTCTGCGcc                                                                                                                            >  2:849467/1‑149 (MQ=255)
                              tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTGGTTTGTCCCTCGGCCCCTTCGCcctttc                                                                                                                       >  2:685747/1‑148 (MQ=255)
                                    gcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTTGTTTGGTCGGCCAACCGCCCCCGCGCCcttggcttgg                                                                                                                  >  2:405330/1‑149 (MQ=255)
                                       aaTGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGCTCGATATTTCGTTTTGTTTGACCACTGGCCCCTGCGCCATTGCATTggggg                                                                                                               >  1:1170963/1‑147 (MQ=255)
                                          gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACt                                                                                                           <  1:438032/150‑1 (MQ=255)
                                           ttAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAACGATATTTCGTTTGGTTCGCCCAATCCCCTCTGCGCCCTTGCATTGGGGGTGCTGt                                                                                                         >  2:1111635/1‑151 (MQ=255)
                                               aTCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCCGCCTCTGCGCCATTGCATTGGGGGGACTGttt                                                                                                       >  2:444018/1‑149 (MQ=255)
                                                tctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCCGCCTCCGCGCCCCTTCATTGGGGGGGCTGttt                                                                                                       >  2:233333/1‑148 (MQ=255)
                                                 ctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGGTATTTTGTTTTGTTTGACCAACGGCCTCTGCGCCCTTGCATTTGGCGTTCTGttttta                                                                                                    >  2:448788/1‑147 (MQ=255)
                                                 ctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCCCTGCGCCCTTGCATTGGGCGTACTGTTTTaa                                                                                                    >  1:323744/1‑150 (MQ=255)
                                                   cAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTc                                                                                                      <  1:1089902/146‑1 (MQ=255)
                                                             cTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTCTTTTGGTGGGCCACCCGCCCCTGCGCCATTGCCTTTGGCGTACTGTTTTAAAGTGGGGTGcc                                                                                         >  1:193604/1‑149 (MQ=255)
                                                               gTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGa                                                                                       <  1:212979/149‑1 (MQ=255)
                                                                 tCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAACGAAATTTCGCTTGGTTTTACCATCCGCCTCTGCGCCCTTGCATTTGGGGTACTGTTTTAAAGGGGGGTGGCGATGcc                                                                                   >  1:894575/1‑151 (MQ=255)
                                                                  cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGTGTTATTTTTTTTTGTTCGACCATCGGCGCCCGCGCCCCTGCATTGGGCGGTCTGTTTTCAAgtgggggggcgctg                                                                                     >  2:905627/1‑138 (MQ=255)
                                                                     cACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTTGTTCGACCATCGGCCTCTGCGCCATTTCATTTGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCg                                                                                  >  2:690750/1‑148 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGAt                                                                          <  1:792193/151‑1 (MQ=255)
                                                                            gACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATa                                                                         <  1:662485/150‑1 (MQ=255)
                                                                                 tATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTTGTCCCTCGGCCTCTGCGCCATTGCATTGGGGGTACTGTTTTTAAATAGGGTGCCGGTGCCGCCGGCGGTACTcc                                                                     >  2:381205/1‑149 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCtt                                                              <  2:916315/151‑1 (MQ=255)
                                                                                            gcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTTGGCCATCGGCCTCTGCGCCCTTGCATTGGGGGTACTGTTTCAAAGTGGGGTGCCGATGCCGCTGGCGATAATCCTGAACTTACTGCt                                                        >  1:810038/1‑151 (MQ=255)
                                                                                                tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCgg                                                     <  1:509036/150‑1 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTTGACCATCGGCCTCTGCGCCCTTGCATTGGGGGGACTGTTTTAAAGAGGGGGGCCGGTGCCGCTGGCGGTACTCCTGACCCTACTGCTCGgcg                                                   >  1:378680/1‑148 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCAt                                                <  1:408902/151‑1 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCAt                                                <  1:995603/151‑1 (MQ=255)
                                                                                                     cccTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGAAATTTTGTTTTGTTCGACCATCCGCCCCTGCGCCATTTTCTTGGGGGTACTGTTTTTAAGTGGGGTGCCCATGCCGCTGGCGGTCCTCCTGACCTTACTGCTCGGCGCCtt                                               >  2:979572/1‑151 (MQ=255)
                                                                                                      ccTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCCTTGGGGGGACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCAtt                                               >  1:487280/1‑150 (MQ=255)
                                                                                                         aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                            <  2:57390/150‑1 (MQ=255)
                                                                                                             cTAACGATGGTGATTGTCAGTGGCGGGGTTGGTAATTCGTTTTGTTTGACCATCGGCCTCTCCGGCCCTGCATTTGGGGTACCGTTTTCAAATGGGGGGGCGATGCCGCTGGCGGAACTCCCGCCCTTTCTGCCCGGCGGCTTGTGggggg                                       >  2:1151453/1‑150 (MQ=255)
                                                                                                             cTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTCGACCATCGGCCTCTTCGCCCATTCCTTGGGCGGACTGTTTTTAAGTGGGGTGGCGATGCCGCTGGCGGTACTCCTGACCTTTCTGCTCGGCGCGCTGTTCggg                                        >  2:658724/1‑150 (MQ=255)
                                                                                                               aaCGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGc                                       <  2:785248/149‑1 (MQ=255)
                                                                                                                         aTTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGGGGACTGTTTTAAAGTGGGGTGCCGATGCCCCTGGCGGTACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGGTGGTCAACGcc                            >  2:260105/1‑150 (MQ=255)
                                                                                                                              cAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatt                       <  1:543539/150‑1 (MQ=255)
                                                                                                                                 tGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatt                   >  2:871593/1‑151 (MQ=255)
                                                                                                                                        aTCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCtata            <  2:1065226/151‑1 (MQ=255)
                                                                                                                                         tCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAc           >  1:193942/1‑151 (MQ=255)
                                                                                                                                                  tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt    >  2:444465/1‑150 (MQ=255)
                                                                                                                                                   cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTaa  <  1:307813/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387

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: