Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,832,7180TG56.2% 10.5 / 17.5 16intergenic (‑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 (9/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 4.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAA  >  NZ_CP009273/2832650‑2832801
                                                                    |                                                                                   
aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGac                                                                <  2:446675/90‑2 (MQ=255)
aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGaa                                                                <  2:251475/90‑1 (MQ=255)
aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGaa                                                                >  1:404235/1‑90 (MQ=255)
aCGCAGAGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGACCACAAAAATAAAAAATGaa                                                                >  1:269190/1‑90 (MQ=255)
        tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCAc                                                        >  1:193192/1‑90 (MQ=255)
            gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATTTAAAAAATGAACAATTCActctc                                                    >  1:151209/1‑90 (MQ=255)
               ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAAAGAACAATTCACTCTCTTg                                                 >  1:128521/1‑90 (MQ=255)
                 gTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTTAGGGGACCGGGTTCACAAATATAAAAAAAGAACAAATCACTCTCTTTCt                                               >  2:355908/1‑90 (MQ=255)
                 gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATTTAAAAAATGAACAATTCACTCTCTTGCt                                               <  2:151209/90‑1 (MQ=255)
                 gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                               >  1:135269/1‑90 (MQ=255)
                    tGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGGCCCGGTTCCCAAATATAAAAAAAGAAGAAATAACACtctttctttt                                            >  2:460411/1‑88 (MQ=255)
                            gTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGaa                                    >  1:388173/1‑89 (MQ=255)
                             tCCTTGCACAATCGGCGGGAAAAATATTCAGGTGCCCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACa                                   <  2:340425/90‑1 (MQ=255)
                               cctGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                 <  2:118146/88‑1 (MQ=255)
                                ttGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACt                                >  1:397574/1‑90 (MQ=255)
                                                              gACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTaa  >  2:258749/1‑90 (MQ=255)
                                                                    |                                                                                   
ACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAA  >  NZ_CP009273/2832650‑2832801

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