Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,832,7180TG53.3% 12.3 / 13.7 15intergenic (‑95/‑165)ascG/ascFDNA‑binding transcriptional regulator AscG/PTS cellobiose/arbutin/salicin transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (8/0);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTC  >  NZ_CP009273/2832635‑2832804
                                                                                   |                                                                                      
cTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGTGACCCGGTTccca                                                                                  >  2:322531/1‑90 (MQ=255)
         cAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTTAGGTGACCCGTCTGACAAAAATaaaa                                                                         >  1:339226/1‑90 (MQ=255)
                cGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                  <  1:164154/90‑1 (MQ=255)
                     tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAAAGAACAATTc                                                             >  2:12420/1‑90 (MQ=255)
                       tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAGTATAAAAAATGAACAATTCAc                                                           <  1:231858/90‑1 (MQ=255)
                       tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                           <  2:310055/90‑1 (MQ=255)
                         cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTTAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTTActc                                                         >  2:349331/1‑90 (MQ=255)
                           gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAAAGAAAAATTTActctc                                                       >  2:170377/1‑90 (MQ=255)
                                gTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTCAGGGGACCGGGTTCACAAATATAAAAAAAGAAAAAATCACTCTCTTGCt                                                  >  2:56938/1‑90 (MQ=255)
                                gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                  >  1:62175/1‑90 (MQ=255)
                                gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAAAAATTTACTCTCTTGCt                                                  >  1:290395/1‑90 (MQ=255)
                                          tGTCCTTGCACAATCGGCGGGGAAAAAATTCAGGGGACCGGGTTCACAAATATAAAAAATGAAAAATTTACTCTCTTGGTTTTTTTGTGa                                        >  2:113173/1‑90 (MQ=255)
                                                     aaTCGGCGGGGAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCAt                             >  2:320017/1‑90 (MQ=255)
                                                                aaaaTATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGaa                  <  2:79460/90‑1 (MQ=255)
                                                                                cGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTc  <  2:145006/90‑1 (MQ=255)
                                                                                   |                                                                                      
CTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTC  >  NZ_CP009273/2832635‑2832804

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