Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273339,1550AG53.8% 5.5 / 15.6 13G187G (GGA→GGG) yahKNADPH‑dependent aldehyde reductase YahK
Reads supporting (aligned to +/- strand):  ref base A (0/6);  new base G (7/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGT  >  NZ_CP009273/339072‑339199
                                                                                   |                                            
gggATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATg                                        <  2:149225/90‑1 (MQ=255)
gggATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCCGGGGCataag                                        >  2:390480/1‑88 (MQ=255)
     caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATgggggt                                   >  2:19992/1‑88 (MQ=255)
     caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGGATATtggggg                                   >  1:423901/1‑88 (MQ=255)
     caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGGAAATgggggt                                   >  1:35808/1‑88 (MQ=255)
     caccacGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGGCTGGGGCATATgggggt                                   >  2:243918/1‑88 (MQ=255)
            tATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATg                                        >  1:133157/1‑78 (MQ=255)
            tATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATg                                        <  2:133157/78‑1 (MQ=255)
                        cGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATg                <  1:168294/90‑1 (MQ=255)
                              tGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGca          <  2:280355/90‑1 (MQ=255)
                                    gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGGCATATGGGGATTAAGCTGGccc                        >  1:130079/1‑70 (MQ=255)
                                    gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGccc                        <  2:130079/70‑1 (MQ=255)
                                      cgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATgtggt  <  1:1443/90‑1 (MQ=255)
                                                                                   |                                            
GGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGT  >  NZ_CP009273/339072‑339199

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