Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,856,0900AC53.8% 3.1 / 12.4 13W48G (TGG→GGG) yeaCYeaC family protein
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base C (7/0);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGTTTGCTCTTCATGACCATCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGCTTGCCACAGCATCACCAGTTGCAGGCAGTTCTCTTTTTGTTCCTCGGTCAACGCAACGCCATCAGGCCATTTCCCC  >  NZ_CP009273/1856012‑1856160
                                                                              |                                                                      
gctgTTTGCTCTTCATGACCATCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGGGTGGCCCAGCATCACCa                                                             >  2:402517/1‑90 (MQ=255)
gctgTTTGCTCTTCATGACCATCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGCTTGCCCCAACATCAccc                                                             >  2:17967/1‑89 (MQ=255)
               tGACCATCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGGTTTGCCCACCCTCACCACTTgccgggaggttt                                              >  1:134246/1‑86 (MQ=255)
                    aTCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGCTTGCCACAg                                                                     >  1:275020/1‑62 (MQ=255)
                    aTCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGCTTGCCACAg                                                                     <  2:275020/62‑1 (MQ=255)
                           cATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGCTTGCCACAGCATCACCAGTTGCAGGCAGTTCTCTTTTTGTTCCt                                  <  1:330061/90‑1 (MQ=255)
                               tgtgTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGGTTTCCCCAGCATCACCCATTGCAGGCAGGTTTCTTTTTTTTTCTCCGt                              >  1:160627/1‑90 (MQ=255)
                                          tCATATGCTGTGCTTCGATATTATGGCGGGGTTTGCCCACCATCACCAGGTTTAGGGAGTTTTTTTTTTTTTTCTTGGgccacgccacac                   >  1:334387/1‑88 (MQ=255)
                                                   gtgCTTCGATATTATGGCGGGCTTTCCCCAACATCACCAGTGGCAGGCCGTTCTCTTTTTTTTCCCCGGGCAACGCCAAGCCATCAGGcc          >  2:92218/1‑90 (MQ=255)
                                                   gtgCTTCGATATTATGGCGGGCTTGCCCCAGCATCACCAGTTGCAGGGAGTTCTCTTTTTTTTTCTCGGGCAACGCCACGCCATCAGGcc          >  2:116908/1‑90 (MQ=255)
                                                   gtgCTTCGATATTATGGCGGGCTTGCCACAGCATCACCACTTGCAGGGAGGTCTCTTTTTTTTTCTCGGGCCACGCAACGCCCACAgggc          >  1:8840/1‑88 (MQ=255)
                                                      cTTCGATATTATGGCGGGCTTGCCACAGCATCACCAGTTGCAGGCAGTTCTCTTTTTGTTCCTCGGTCAACGCAACGCCATCAGGCCAtt       <  1:157043/90‑1 (MQ=255)
                                                           atatTATGGCGGGCTTGCCACAGCATCACCAGTTGCAGGCAGTTCTCTTTTTGTTCCTCGGTCAACGCAACGCCATCAGGCCATTTcccc  <  2:8840/90‑1 (MQ=255)
                                                                              |                                                                      
GCTGTTTGCTCTTCATGACCATCTGACCATTTGTGTCAATTGTCATATGCTGTGCTTCGATATTATGGCGGGCTTGCCACAGCATCACCAGTTGCAGGCAGTTCTCTTTTTGTTCCTCGGTCAACGCAACGCCATCAGGCCATTTCCCC  >  NZ_CP009273/1856012‑1856160

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