Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT44.2% 43.2 / 18.5 43S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (6/18);  new base T (19/0);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 2.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387
                                                                                                                                                  |                                                                                                                                                       
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGttt                                                                                                                                                     >  2:127083/1‑151 (MQ=255)
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGtt                                                                                                                                                      <  2:807668/150‑1 (MQ=255)
    aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGtt                                                                                                                                                 <  2:579176/151‑1 (MQ=255)
       gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGtt                                                                                                                                                 >  2:558514/1‑148 (MQ=255)
            tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGTTATTTCTTTTTGtt                                                                                                                                                 >  1:4178/1‑143 (MQ=255)
             gTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGACGATATTTTTTTTGGTTCGGACATc                                                                                                                                         >  2:90811/1‑150 (MQ=255)
                     tGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAACTTTTTTTTGTTTTGTTTGGCCATCCGcc                                                                                                                                     >  2:667171/1‑146 (MQ=255)
                              tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgc                                                                                                                       <  2:22801/150‑1 (MQ=255)
                               aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATTTTTTGTTTTGTTTGGCCATCGGGCTCCGCGCccttgcct                                                                                                                     >  1:414147/1‑149 (MQ=255)
                                  ccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTTGACCATCCGCTCTTGCGCcagtgcctt                                                                                                                    >  1:790804/1‑149 (MQ=255)
                                          gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACt                                                                                                           <  1:731167/150‑1 (MQ=255)
                                          gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTg                                                                                                          <  2:79131/151‑1 (MQ=255)
                                            tataTCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGCGGTATTTTGTTTTGTTTGACCACCCGGCTCTTCGCCCTTTGCTTGGGCGTAATGtt                                                                                                        >  1:86368/4‑151 (MQ=255)
                                                   cAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACTATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgt                                                                                                 <  2:781303/151‑1 (MQ=255)
                                                    aaTATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGCGGTATTTTTTTTTGTTTGCCCATCGGCCCCTGCGCCCTTGGATTTGGCGGACTGtttttaa                                                                                                   >  1:527306/1‑144 (MQ=255)
                                                      tatGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTCGACCATCGGCCTCTTCGCCCCTTCATTGGGGGTACTTTTTCaaaaggg                                                                                               >  2:539762/1‑146 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtggt                                                                                              <  2:110478/150‑1 (MQ=255)
                                                        tGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCCTTTTGTTCGACCATCGGCCTCTGGCCCCTTGCCTTTGGCGGGCTGTTTCTAAGTGgggt                                                                                            >  1:380539/1‑151 (MQ=255)
                                                            gCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGTTATTTTTTTTTGTTTGGCCATCCGCCCCTTCGCCACTTGATTTGGGGGTCTGttttttaatggggt                                                                                            >  1:165899/1‑136 (MQ=255)
                                                            gCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCg                                                                                        <  2:10746/151‑1 (MQ=255)
                                                                ttCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCATTTGGGGTACTGTTTTAAAGGGGGGTGCCGATg                                                                                     >  2:529634/1‑150 (MQ=255)
                                                                ttCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTCGACCATCCGCCTCTGCGCCATTGCCTTGGGGGGACTGTTTTTAAGTGGTGTGCCGATg                                                                                     >  2:49236/1‑150 (MQ=255)
                                                                    gCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTTTACCATCCGCCCCTTCGCCCTTGGACTTGGGGGACTTTTTTCAAGTGGGGtgccgctgccgc                                                                                 >  1:540143/1‑150 (MQ=255)
                                                                       ccAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTCGACCATCTGCCCCTGCGCCCTTGCATTGGGGGTACTGTTTTTAAGAGGGGTgccgctgccgct                                                                                >  2:378252/1‑148 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTTGACCATCGGCCTCTGCGCCCTTGCATTTGGGGGACTGTTTCAAAAAGGGGGGGCGATGCCGCCggggg                                                                           >  2:84210/1‑147 (MQ=255)
                                                                           tGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTTGACCATCGGCCTCTGCGCCCTTGCCTTGGGCGGACTGTTTTAAAAGGGGGTGCCGATGCCGCTgggg                                                                            >  2:541117/1‑146 (MQ=255)
                                                                                ttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTcc                                                                     <  2:276814/150‑1 (MQ=255)
                                                                                 tATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCt                                                                    <  2:95725/150‑1 (MQ=255)
                                                                                     tGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTTACCCTCCGCCCCTTCGCCATTGCATTTGGGGTACTGTTTTTAAGTGGGGGGCCGATGCCGCCGGCGATACTCCTGAcc                                                                >  1:743875/1‑150 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCt                                                               <  1:211154/150‑1 (MQ=255)
                                                                                          tggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCCTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACt                                                           <  1:540600/150‑1 (MQ=255)
                                                                                                tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTGGTTCTTCCCTCGGCCTCTGCGGCCTTGGATTTGGGGTGCTGTTTTTAAGGGGTGTGCCGGTGCCGCTGGCGGTACCCCCGACCtt                                                              >  1:297874/1‑141 (MQ=255)
                                                                                                            gCTAACGATGGTGATTGTCAGTGGCGGGAGCGATATTTTTTTTTGTTCGACCATCGGCCTCTGCGCCAGTGCATTTGGCGTGCTGTTTCAAAGGGGTGTGCCGAGGGCGCTGGCGGTACTCCCGGCCCTTCTTCTCCGCGCa                                                 >  1:40505/1‑142 (MQ=255)
                                                                                                                aCGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGa                                    <  1:70951/151‑1 (MQ=255)
                                                                                                                   aTGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCACATTGTGCGGGCTGATCa                                 <  1:654235/151‑1 (MQ=255)
                                                                                                                   aTGGTGATTGTCAGTGGCGGGATCGATAATTCGTTTGGTTCTAACCACGGCCCCTGCGCCCTTGCATTTGGCGGACTGTTTTAAAGGGGTGTGGCGATGCCCCCCGCGGTACTCCTGGACTTACTGCTCTGCGCCTTGGTCGGGccggt                                   >  1:524697/1‑145 (MQ=255)
                                                                                                                         aTTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCg                           <  2:121352/151‑1 (MQ=255)
                                                                                                                               aGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaa                     <  1:265548/151‑1 (MQ=255)
                                                                                                                                  ggCGGGATCGATATTTCGCTTGGTTCGGCCATCGGCCTCTGGGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGGTGGCGATACTCCTGACCTTACTGCTCGGGGCATTGTGGGGGCTGATCAACGCCGGattaatt                   >  1:774294/1‑150 (MQ=255)
                                                                                                                                          cGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAcc          >  1:69168/1‑151 (MQ=255)
                                                                                                                                             tatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa        >  1:763897/1‑150 (MQ=255)
                                                                                                                                                ttCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGt     <  2:570222/150‑1 (MQ=255)
                                                                                                                                                 tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt    <  2:774379/150‑1 (MQ=255)
                                                                                                                                                  cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTaa  <  2:775939/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 < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: