Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,329,9480TG53.8% 12.6 / 17.1 13intergenic (+24/‑140)BW25113_RS25560/acnAhypothetical protein/aconitate hydratase AcnA
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base G (7/0);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTA  >  NZ_CP009273/1329864‑1330033
                                                                                    |                                                                                     
gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCa                                                                                  >  1:65808/1‑90 (MQ=255)
  gcgcGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtt                                                                                <  2:65808/90‑1 (MQ=255)
               gTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACaa                                                                   >  2:484069/1‑90 (MQ=255)
                                                 aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGTTTTTCCTCTTTTTTCAATTTTGGTTGTTTTc                                 >  1:237128/1‑90 (MQ=255)
                                                 aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGTTTTTCCTCTTTTTTCAATTTTGGGTGTTATc                                 >  1:24046/1‑90 (MQ=255)
                                                          tGATCCTTGACCCCGCTGCGGCGGGGGTGTTCTTTGCTTTTCCACAAAGTTTTTCCTCTTTTTTCAATTTGGGGTTTTTTTAAAACGTTa                        >  1:295371/1‑90 (MQ=255)
                                                           gATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGTTTTTCCTCTTTTTTTAATTTGGGGTGTTTTCAAATCGtttc                       >  1:486376/1‑88 (MQ=255)
                                                                   aCCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGtt               >  2:408491/1‑90 (MQ=255)
                                                                   aCCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGtt               <  2:424159/90‑1 (MQ=255)
                                                                   aCCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCACAAGGTTTTTCCTCTTTTTTTAATTTGGGGTTTTTTCAAATCGTTACGCGATGtt               >  1:495378/1‑90 (MQ=255)
                                                                        gCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt          <  1:408491/90‑1 (MQ=255)
                                                                          tGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa        <  2:425256/90‑1 (MQ=255)
                                                                                ggggTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTa  <  2:55203/90‑1 (MQ=255)
                                                                                    |                                                                                     
GGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTA  >  NZ_CP009273/1329864‑1330033

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