Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,598,9040TC25.0% 54.7 / 13.4 36S436G (AGT→GGT) lsrKautoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base C (5/4);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTA  >  NC_000913/1598766‑1599046
                                                                                                                                          |                                                                                                                                              
ccAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTCATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCACGAAAGATCAGCGAGAAATTGCCTCCAt                                                                                                                                            >  1:114097/1‑143 (MQ=255)
 cAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATa                                                                                                                                           >  1:618293/1‑143 (MQ=255)
 cAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAAGGGGACATCCTAATGCAGTCGCTTCTTTGACCGCCGGACTGTTGCCGGGTAATCCCGAGACATCAGCGAGACTTCGACTCCATa                                                                                                                                           >  1:4175/1‑143 (MQ=255)
     tgttgtgcGGTCTCTGCCAGTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAAttt                                                                                                                                       <  2:340333/136‑1 (MQ=255)
              gTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGa                                                                                                                              >  1:500530/1‑143 (MQ=255)
                  cTGCCATTGATGAAAAACTTCCGGCACCGACGCCAGCGGCAATGGCACATCCTGATGCAGTGGCTTCTTGGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCt                                                                                                                          <  2:172210/143‑1 (MQ=255)
                          gatgaAAAAATTCCGGCACCGACGCCATCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGc                                                                                                                  <  1:217590/143‑1 (MQ=255)
                             gaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCaaa                                                                                                               >  1:163661/1‑143 (MQ=255)
                              aaaaaaTTCCGGCACCGACGCCAGCTGCAATGGCAGATCCTAATGCAGTGGCTGCTTTGAGCACCGGAAAATTGACGGGTAATCCCGAGACATAAGAGAGAATTTGCCTCCATAATTTCCCTTTTGAAACTCCGCCCGAAAAg                                                                                                              >  2:43248/1‑143 (MQ=255)
                                          cACGGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGaa                                                                                                  <  1:430927/143‑1 (MQ=255)
                                              gACGCCAGCTGCAATCGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGAt                                                                                              <  1:152396/143‑1 (MQ=255)
                                               aCGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATg                                                                                             >  1:663731/1‑143 (MQ=255)
                                                ccccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGa                                                                                            <  1:657693/141‑1 (MQ=255)
                                                  ccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGTGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCGCCGCCTGCAAAGACTAACGAGGAAGGATGAat                                                                                          <  1:669126/143‑1 (MQ=255)
                                                  ccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGCCATCAGCGCGAATTTGACCCCATAATTTCTCCTTTGACCCTCCGCCTGCAAACACTGAcactcaacgaacaca                                                                                          >  2:598425/1‑129 (MQ=255)
                                                  ccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAat                                                                                          <  2:473669/143‑1 (MQ=255)
                                                        gCACTGGCACATCCTAATGCAGTGGCTTCTTTGACCACTGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCga                                                                                    <  2:711940/143‑1 (MQ=255)
                                                        gCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCga                                                                                    <  2:569965/143‑1 (MQ=255)
                                                        gCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCga                                                                                    <  1:316611/143‑1 (MQ=255)
                                                         cAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGGCTAACGATGAAGGATGAATATTCgag                                                                                   <  1:234525/143‑1 (MQ=255)
                                                          aaTGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATCATTTCCCTTTTGAACCTCCGCCTGCACAGTCGAACGCTGACGGCTGAATGTTCgcga                                                                                  >  1:695112/1‑143 (MQ=255)
                                                          aaTGGCACATCCTAATGCAGTGGCGTCTCTGCCCAATGCACAAGTGACGTTTAATCCCGAGACAGCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCgaga                                                                                  <  2:760907/143‑1 (MQ=255)
                                                               cacaTCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGACTTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAACGCCTAAGGATGCACGATGAATCTCCGACAAATCa                                                                             >  2:388909/1‑143 (MQ=255)
                                                                           cAGTGGCTTCTTTGACCACCGGCATATTGACGGGTAATCCCGCGACATCAGCGAGCGTTTGACCCCATAATTTACCTTTTGAACCTACGCCTGCAGAGACTAACGATGAAGGATGAATCTTCGCGAACGCAGAAATTggctgc                                                                 >  2:508111/1‑143 (MQ=255)
                                                                             gTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAACATTCGAGAAATCAGCAATTTGCTGCaa                                                               <  1:478125/143‑1 (MQ=255)
                                                                                 cttcttTGACCACCGGAATATTGACGGGTAATCCCGAGACAGCAGCGAGAATTTGACTCCATAATTTCCCTCTTGCACTTGCTCTTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTa                                                           <  2:725774/143‑1 (MQ=255)
                                                                                            accGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAAGCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCCAGAAATCAGCACTTTGCTGCAACCTACACGCTCATAc                                                >  1:660672/1‑143 (MQ=255)
                                                                                                gAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTACCCTAGTGAAGCTGAGACTGGTGACACCAACGATGAAGGATGAATATTCGCGAAATCAGCACTTTGCTGTACGTTACACGCTGATAGAATc                                            >  2:277010/1‑143 (MQ=255)
                                                                                                          gggTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGATACACGCCGATACAATCACCGCAtttt                                  >  2:488152/1‑143 (MQ=255)
                                                                                                               aTCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTcc                             >  1:425346/1‑143 (MQ=255)
                                                                                                                  ccGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATGCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAgc                          >  1:213321/1‑143 (MQ=255)
                                                                                                                        cATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACgg                    >  1:761433/1‑143 (MQ=255)
                                                                                                                        cATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACgg                    >  1:628009/1‑143 (MQ=255)
                                                                                                                            aGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAAACTAACGATGAAGGATCACTCTTCGAGAAATCAGCAATTCGCTGCAAGTTACACGCTCACACAATCGCAGCCTTCTCTTCCACCGAAccggacg                >  2:731219/1‑137 (MQ=255)
                                                                                                                                 gaATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTc           <  1:566032/143‑1 (MQ=255)
                                                                                                                                  aaTTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCg          >  2:578840/1‑143 (MQ=255)
                                                                                                                                          tCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCACTTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTa  >  1:194196/1‑143 (MQ=255)
                                                                                                                                          tCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTa  >  1:609413/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CCAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTA  >  NC_000913/1598766‑1599046

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