Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT40.0% 57.7 / 11.8 45S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (8/19);  new base T (18/0);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.10e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

TTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  NC_000913/1604093‑1604386
                                                                                                                                                |                                                                                                                                                      
ttATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGATATTTTGttt                                                                                                                                                    >  2:212762/1‑149 (MQ=255)
     gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCg                                                                                                                                              <  1:391762/150‑1 (MQ=255)
             cGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCgg                                                                                                                                      <  1:311101/150‑1 (MQ=255)
              gCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGcc                                                                                                                                    <  2:433029/151‑1 (MQ=255)
                  ttGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTGGTTCGTCCATCGGcccc                                                                                                                                  >  2:674562/1‑147 (MQ=255)
                       gCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGcc                                                                                                                            >  1:197776/1‑150 (MQ=255)
                        cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATGGGTATTTTGTTTGGGTCGACCATCGGCCCCTGCGcc                                                                                                                            >  2:178577/1‑149 (MQ=255)
                          aTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATCGGCCTCCGCGCcatt                                                                                                                         >  1:269701/1‑150 (MQ=255)
                                          tAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCCTTGCCATTGGGGTACtttt                                                                                                       >  2:808707/1‑151 (MQ=255)
                                            gATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGCCTCTGCGCCCTTGCATTTGGCGGACTGttt                                                                                                      >  1:322505/1‑150 (MQ=255)
                                                  aaTATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGTTATTTTCTTTTGTTCGACCATCGGGCTCTGCGCCCTTGCATTTGGGGTACTGtttttcagt                                                                                                >  2:662723/1‑144 (MQ=255)
                                                        ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGCCCCTCGGCCTCTGCGCCATTGCATTGGGGGGGCTGTTTTTAAGTGGggtg                                                                                          >  2:668319/1‑150 (MQ=255)
                                                             gTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTTGACCATCGGCCTCTGCGCCCTTTCATTGGGGGGACTGTTTCAAAGTGGGGTGCCGATg                                                                                    >  1:400794/1‑151 (MQ=255)
                                                                  gCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCt                                                                               <  1:716199/151‑1 (MQ=255)
                                                                     ccAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGc                                                                            <  2:76297/151‑1 (MQ=255)
                                                                            cTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCCGCCTCTGCGCCCTTGCATTGGGCGGACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTc                                                                     >  1:842475/1‑151 (MQ=255)
                                                                                 tCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCCTTGCATTGGGCGGACTGTTTTAAAATGGTGTGCCGATGCCGCTGGGGATACTCCTGAc                                                                >  1:119143/1‑151 (MQ=255)
                                                                                      attggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTTGTTCGACCATCGGCCCCTGGGCCATTGCATTTGGCGGACTGTTTTTAAGAGGGGTGCCGCTGCCCCCGGCGGTACACCTGGccct                                                             >  2:347751/1‑147 (MQ=255)
                                                                                             ttgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGTTATTTTGTTTGGGTTGACCAACGGGCCCTGCGGCCTTGCATTGGGGGTGCTTTTTCAAAGTGGGGTCCCGCTGCCGCTGGCGATACTCCCGGCCTTACTGCCCgg                                                    >  2:209853/1‑151 (MQ=255)
                                                                                               gTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGCGATATTTTGTTTTGTTTGACCCTCGGCCCCTGCGCCCTTGCATTGGGGGGACTGTTTTCAAGTGGGGTGGCGGGGCCGCTGGCGGTACTCCTGACCCTACTGCTCGgcg                                                  >  1:135998/1‑151 (MQ=255)
                                                                                                  gCCCTTTCGCTAACGATGGTGATTGTCCGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGATGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGGCCTTACTGCTCGgcgc                                                 <  2:367516/149‑1 (MQ=255)
                                                                                                       aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCAACGGCCTCTGCGCCATTGCATTTGGCGTACTGTTTTAAAAAGGGGTGCCGATGCCGCTGGCGGTACTCCTTACCCTACTGCTCGGCGCCTtgt                                            >  1:295565/1‑149 (MQ=255)
                                                                                                        ccGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCg                                         <  1:341796/151‑1 (MQ=255)
                                                                                                         cGCTAACGATGGTGATTGTCAGTGGCGGGGTCCATATTTTCTTTTGTTCGACCCCCCGCCCCTTGGGCCTTGCCTTGGGGGGACTGTTTTTAAGGGGGGGGGCCGTGCCGCTTGCGGTACTCCTGACCCTTTTGCTCGGCGCCTTTTTCgg                                        >  2:77646/1‑151 (MQ=255)
                                                                                                         cGCTAACGATGGTGATTGTCAGTGGCGGGAGCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGGCCTTACTGCGCGGCGCATTGTGCgg                                        <  2:593213/151‑1 (MQ=255)
                                                                                                         cGCTAACGATGGTGATTGTCAGTGGCGGCGTTCATTTTTCGTTTTTTTCGCCCAGCCGCCCCCCGGCCACTGCATTGGGCGTGCTCTTTTAAACTGGGGTGGCGGTGCCGCTTGCGGTACTCCTGACCTTCATGCTccg                                                    >  2:412072/1‑137 (MQ=255)
                                                                                                           cTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTTCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCggg                                       <  2:456993/150‑1 (MQ=255)
                                                                                                             aaCGATGGTGATTGTCAGTGGCGGGGTCGGTATTTCCTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCCTTGGGCGTACTGTTTTAAAATGGTGTGGCGATGCCGCTGGCGGTACTCCTGTCCTTTCTGCTCGGGGCATTGTGgggg                                       >  2:415273/1‑148 (MQ=255)
                                                                                                               cGATGGTGATTGTCAGTGGCTGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGCCCTTACTGCTCGGCGCATTGTGCGCGCTGa                                   <  2:447437/150‑1 (MQ=255)
                                                                                                               cGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCCTTGGGCGTACTGTTTTAAAGTGGTGTGCCGATGCCGCTGGCGGTACTCCTGACCTTACTGCTCGGCGCCTTGTGCggg                                       >  1:483181/1‑146 (MQ=255)
                                                                                                                gATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGGTCGACCATCGGCCTCTGGGCCCTTGCATTGGGGGGAATGTTTTTAAATGGTGGGGCGATGCCGCCGGCGAGACTCCTCACCTTACTTCTCGGCGCCTTGTGCGGGCTggt                                  >  2:166252/1‑148 (MQ=255)
                                                                                                                 aTGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCa                                <  1:678950/151‑1 (MQ=255)
                                                                                                                    gTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACg                             <  2:462366/151‑1 (MQ=255)
                                                                                                                       aTTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCg                          <  1:468994/151‑1 (MQ=255)
                                                                                                                         tGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCgg                         <  2:20815/150‑1 (MQ=255)
                                                                                                                          gTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGa                        <  1:668481/150‑1 (MQ=255)
                                                                                                                          gTCAGTGGCGGGATCGATATTTCGCTTGGGTCGACCATCCGCCTCTCCGCCATTTCATTTGGCGCAATTGTTTAAAGTGGGGTGCCCATGCCGCCGGCGATAACCCCGGCCTTACCTCCCGCGGCATCGTGCGGGcgcgttcacacccg                         >  2:62822/1‑136 (MQ=255)
                                                                                                                            cAGTGTCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatta                     <  2:851701/151‑1 (MQ=255)
                                                                                                                            cAGTGGCGGGATCGGTATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTTCATTGGGGGTACTGTTTTAAAGGGGGGTGCCGATGCCGCTGGGGGTACCCCTTACCTTACTTCTCGGCGCCCTGTGCGGGGGGGTCAACCCCCGatt                      >  2:613941/1‑150 (MQ=255)
                                                                                                                            cAGTGGCGGGATCGATATTTCGTTTGGTTCGAGCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatta                     <  2:413660/151‑1 (MQ=255)
                                                                                                                            cAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatta                     <  1:58971/151‑1 (MQ=255)
                                                                                                                                        cGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAcc         >  2:482169/1‑151 (MQ=255)
                                                                                                                                        cGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAcc         >  2:458548/1‑151 (MQ=255)
                                                                                                                                         gATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAcc         >  2:839849/1‑150 (MQ=255)
                                                                                                                                               tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt   <  2:520623/150‑1 (MQ=255)
                                                                                                                                                cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTa  <  1:742067/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTA  >  NC_000913/1604093‑1604386

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