Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT52.9% 9.9 / 35.0 34S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (5/11);  new base T (18/0);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 1.53e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-03

TATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCA  >  NC_000913/1604094‑1604380
                                                                                                                                               |                                                                                                                                                 
tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGttt                                                                                                                                               >  2:347231/1‑148 (MQ=255)
 aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGtt                                                                                                                                          <  2:233672/151‑1 (MQ=255)
  ttGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGttttg                                                                                                                                             >  1:588723/1‑146 (MQ=255)
            cGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTTGAcc                                                                                                                                     >  1:450648/1‑145 (MQ=255)
               aTTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTCGACCCTCGGcc                                                                                                                              >  1:550731/1‑149 (MQ=255)
                    gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTgcg                                                                                                                       >  2:267794/1‑150 (MQ=255)
                        aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGAGATTTTGTTTGTTTCTGCCAACgggc                                                                                                                              >  1:535818/1‑138 (MQ=255)
                        aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcatt                                                                                                                  <  1:83912/151‑1 (MQ=255)
                            aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgca                                                                                                               <  2:182666/150‑1 (MQ=255)
                                   gAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTggggg                                                                                                        >  1:524825/1‑148 (MQ=255)
                                   gAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGt                                                                                                       <  1:598787/151‑1 (MQ=255)
                                     aTGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTAc                                                                                                     <  2:524828/151‑1 (MQ=255)
                                       gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGCCTCTGCGCCCTTGCATTGGGCGGACt                                                                                                    >  1:643454/1‑150 (MQ=255)
                                        ttAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTCTTTTTGGTCGACCATCGGCCTCTGCGCCATTTCATTTGGGGGACt                                                                                                    >  2:648207/1‑149 (MQ=255)
                                         tAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGACGGTATTTTGTTTTGGTCGTCCATCGGCCTCCGCGGCCTTTCATTTGGGGTACtttt                                                                                                 >  2:291054/1‑151 (MQ=255)
                                             tctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCCTCGGGCTCTGCGCCCTTTCATTtgggggaatggtt                                                                                                >  2:52233/1‑138 (MQ=255)
                                               tcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCCACGGCCCCTGCGCCCCTTCCTTGGGCGTACTGTTTcca                                                                                             >  2:433105/1‑147 (MQ=255)
                                               tcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGGGATTGTCAGTGGCGGGGTCGGTTTTTTTTTTGGGTTGGCCAACTGGCCCTTCGCCGTTGCAGTGGGCGGACTGt                                                                                                  >  2:314220/1‑144 (MQ=255)
                                                    atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTCTTTTGTTTGACCATCGGCCTCTGCGCCCCTGCATTGGGCGGGCTGTTTCCAAgtgg                                                                                        >  2:468268/1‑149 (MQ=255)
                                                      gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGgtgt                                                                                     <  1:282977/150‑1 (MQ=255)
                                                       ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGGTTCGCCATCCGCCTCTGCGGCCATTCATTGGGCGTACTGTTTTAAAatgg                                                                                        >  1:176909/1‑146 (MQ=255)
                                                          cTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAACGATATTTTGTTTTGTTTGACCATCGGCCTCTGCGCCATTGCATTTGGCGGACTGTTTCAAAATGGGGTGcc                                                                                  >  2:127460/1‑149 (MQ=255)
                                                              tCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGGTCGGCCATCGGCCCCTGCGCCCTTGCCTTGGGGGGACTGTTTTAAAGTGGTGTGCCGATGcc                                                                            >  2:282896/1‑151 (MQ=255)
                                                                          aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCCTTTGGCGGACTGTTTTTAAGTGGTGTGGCGATGCCGCTGGCGGTAc                                                                 >  1:172127/1‑150 (MQ=255)
                                                                            tttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCCTTGCATTGGGGGTACTGTTTTAAAATGGGGTGCCGATGCCGCTGGCGGTACTcc                                                              >  1:319333/1‑151 (MQ=255)
                                                                                  tGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCCACCATCCGCCTCTGCGCCATTGCATTGGGCGTACTGTTTTAAAGAGGGGTGGCGGCGCCGCTGGGGATACTCCCGAcc                                                         >  2:314355/1‑150 (MQ=255)
                                                                                   gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCtt                                                       <  2:482395/151‑1 (MQ=255)
                                                                                     attggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGATATTTCGTTTTGGTCGACCATCCGCCTCTGCGCCACTGCATTGGGGGTGCTGTTTTTAAGGGGTGTGCCGATGCCGCTTGCGGTACTCCTGACCtt                                                       >  1:654481/1‑149 (MQ=255)
                                                                                               tCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGACATTTCGCTTGGGTTTTCCACCGGTCCCTGGACCAATGCCGTGGGCGGACTGTTTTTAAGGGGGGTGCCGGTGCCGCGGGCGTTACCCCTGATCTTAATTCTCgggg                                            >  1:559496/1‑148 (MQ=255)
                                                                                                 gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCAt                                         <  1:267794/151‑1 (MQ=255)
                                                                                                     tACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGAAATTTTGTTTTGTTTGGCCCTCGGCCTCTGCGCCAATGCATTGGGCGTACTGTTTTTAAGAGGGGTGCCGATGCCGCTGGCGATACTCCCGCCCTTACTGCTCCGCGCATtttg                                     >  2:275740/1‑151 (MQ=255)
                                                                                                                        tGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGGCCTTACTGCTCGGCTCATTGTGCGGGCTGATCAACGCCgg                   <  2:52733/150‑1 (MQ=255)
                                                                                                                         gTCAGTGGCGGGATCGATATTTCGTTTGGTTCGCCCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGa                  <  1:429620/150‑1 (MQ=255)
                                                                                                                                        gATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa  <  2:503368/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCA  >  NC_000913/1604094‑1604380

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