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

GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGT  >  NZ_CP009273/1329865‑1330025
                                                                                   |                                                                             
ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAt                                                                         >  1:122161/1‑90 (MQ=255)
      ttttCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctt                                                                   <  1:345343/90‑1 (MQ=255)
                        ttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCTTTTGCTTTTCCACAAGGTTTctccc                                                 >  1:609101/1‑89 (MQ=255)
                        ttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTctcct                                                 >  2:340850/1‑90 (MQ=255)
                         tctGGGCCGATCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTctcctc                                                <  1:842829/90‑1 (MQ=255)
                                 gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGTTTTTCCTCTTTTTTCa                                        >  2:573586/1‑90 (MQ=255)
                                 gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTTTCa                                        >  1:7041/1‑90 (MQ=255)
                                      cctcTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTg                                   >  2:169819/1‑90 (MQ=255)
                                         ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgt                                <  1:573586/90‑1 (MQ=255)
                                          tGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgtt                               <  1:22541/90‑1 (MQ=255)
                                                      cACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTGGGGTGTTTTCAAATCg                   >  1:317213/1‑90 (MQ=255)
                                                      cACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTCTCCCCTTTTTTCAATTTTGGGTGTTTTCAAATCg                   >  1:179090/1‑90 (MQ=255)
                                                        cTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTTTCAATTTGGGGTGTTTTCAAATCGtt                 >  2:99528/1‑90 (MQ=255)
                                                            tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgc             <  1:340850/90‑1 (MQ=255)
                                                                  aCCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCCCAAGGGTTTTCCTCTTTTTTCCATTTGGGGTGTTTTTAAATCGTTACGCGATGtt       >  1:438463/1‑90 (MQ=255)
                                                                       gCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt  <  2:19938/90‑1 (MQ=255)
                                                                                   |                                                                             
GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGT  >  NZ_CP009273/1329865‑1330025

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: