Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,598,4940CT55.6% ‑1.8 / 18.1 18S64L (TCG→TTG) lsrDautoinducer 2 ABC transporter permease LsrD
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base T (10/0);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 1.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACC  >  NZ_CP009273/1598409‑1598582
                                                                                     |                                                                                        
cTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGttt                                                                                      >  2:107529/1‑90 (MQ=255)
  gTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGttttg                                                                                    >  1:127004/1‑88 (MQ=255)
  gTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGttttg                                                                                    >  2:160383/1‑88 (MQ=255)
       gCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGCGGTATTTTTTTTTGtttgg                                                                               >  2:195837/1‑87 (MQ=255)
            aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTTGTTTCACCATc                                                                          >  2:212315/1‑90 (MQ=255)
             gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCg                                                                         <  1:319928/90‑1 (MQ=255)
                   ttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTGGTTTGACCATCGGCctct                                                                   >  1:241659/1‑90 (MQ=255)
                         gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGTGATATTTTGTTTTGTTCGGCCATCGCCCTCTGCGccc                                                             >  2:464722/1‑89 (MQ=255)
                               gcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcat                                                       <  1:195837/90‑1 (MQ=255)
                                   tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGCGATATTTTGTTTTGTTCGACCCCCCGCCCCTTCGCCCCTGCATTtgg                                                   >  2:38547/1‑90 (MQ=255)
                                   tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTCGACCATCGGGCTCTGCCCCCTTGCATTggg                                                   >  2:221639/1‑90 (MQ=255)
                                        cccTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTAc                                              <  1:19070/90‑1 (MQ=255)
                                             ccGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGttt                                         <  1:19021/90‑1 (MQ=255)
                                                       tGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGGGGACTGTTTTAAAgggggg                               >  1:273307/1‑88 (MQ=255)
                                                                            tCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATAc          <  1:302164/90‑1 (MQ=255)
                                                                               atatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTcc       <  1:160383/90‑1 (MQ=255)
                                                                                  tttCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGa    >  2:452165/1‑90 (MQ=255)
                                                                                    tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGAcc  >  2:337476/1‑90 (MQ=255)
                                                                                     |                                                                                        
CTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACC  >  NZ_CP009273/1598409‑1598582

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