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

TTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAA  >  NZ_CP009273/2832636‑2832801
                                                                                  |                                                                                   
tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAGAAACATTCTCGTGACCGGGTTcccga                                                                              >  1:352668/1‑86 (MQ=255)
tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGTGACCGGGTTCACaa                                                                              >  2:742113/1‑90 (MQ=255)
        cAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATaaaa                                                                      >  2:667281/1‑90 (MQ=255)
        cAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGACCGGGTTCACAAATATaaaa                                                                      >  1:224857/1‑90 (MQ=255)
              aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGaa                                                                <  2:410980/90‑1 (MQ=255)
                  aGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTTAGGTGACCGGGTTTACAAATATAAAAAATGAACAAt                                                            >  1:61534/1‑90 (MQ=255)
                    tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTTCGGGGACCCGTTTCACAAAAATAAAAAATGAACAAttt                                                          >  1:328766/1‑89 (MQ=255)
                    tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTCAGGAGACCGGGGGCCCCAAAATAAAAAAAGAACAATTc                                                          >  1:393326/1‑90 (MQ=255)
                    tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGCGCCCGGATTCACAAATATAAAAAAGGAACaaata                                                          >  1:410273/1‑87 (MQ=255)
                      tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                        <  2:172029/90‑1 (MQ=255)
                        cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActc                                                      <  2:328766/90‑1 (MQ=255)
                        cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGACCGGGTTCACAAATATAAAAAATGAAAAATTCActc                                                      >  1:671451/1‑90 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                               <  1:837413/79‑1 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                               >  2:837413/1‑79 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCCCAAATATAAAAAATGAACAATTCActctc                                                    >  2:108556/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAAAAATTCACTCTCTTGCt                                               >  2:415459/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCCGGTTCACAAAAATAAAAAAAGAAAAAATCACTCTCTTGCt                                               >  1:120645/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAATGAAAAATTTACTCTCTTTCt                                               >  2:406574/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAATGAAAAATTCACTCTCTTGCt                                               >  2:180747/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTAAGGTGACCCGGTTCACAAATATAAAAAATGAACAACTCACTCTCTTGCt                                               >  2:427911/1‑90 (MQ=255)
                                   gCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTAtt                                           <  1:831069/90‑1 (MQ=255)
                                   gCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTAtt                                           <  2:741374/90‑1 (MQ=255)
                                        ttGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTg                                      <  2:61534/90‑1 (MQ=255)
                                             cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                 <  2:498930/90‑1 (MQ=255)
                                             cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAAAGAACAAATTACTCTCTTTCTTTTTTTGTGAAAAc                                 >  1:283645/1‑90 (MQ=255)
                                                    aaTCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCAt                          >  2:697216/1‑90 (MQ=255)
                                                                 aaTATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAc             <  1:568230/90‑1 (MQ=255)
                                                                 aaTATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAc             <  1:236483/90‑1 (MQ=255)
                                                                            gACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTaa  >  2:246836/1‑90 (MQ=255)
                                                                                  |                                                                                   
TTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAA  >  NZ_CP009273/2832636‑2832801

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