Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273333,7850GT52.6% 20.7 / 36.0 19V433L (GTG→TTG) fdrAacyl‑CoA synthetase FdrA
Reads supporting (aligned to +/- strand):  ref base G (0/9);  new base T (10/0);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCT  >  NZ_CP009273/333701‑333862
                                                                                    |                                                                             
ccGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTc                                                                          <  2:34682/90‑1 (MQ=255)
         gACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGGGGTGGGTTTTTTGCTCGGGTTTTGa                                                                 >  1:634811/1‑90 (MQ=255)
         gACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGGGGTCGTGGTGGGTTTTTTGGTCGGGTTTggg                                                                 >  2:13807/1‑89 (MQ=255)
         gACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGGGGTCGGGGTGGGTTTTTTGGTCGGGTTTTGa                                                                 >  2:49080/1‑90 (MQ=255)
          aCCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGAt                                                                <  2:343347/90‑1 (MQ=255)
          aCCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGAt                                                                <  1:386632/90‑1 (MQ=255)
           ccccACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATc                                                               <  1:13807/90‑1 (MQ=255)
           ccccACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATc                                                               <  2:133844/90‑1 (MQ=255)
               accaACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTTTGCTCGGGTTTGGATCGCAt                                                           >  1:125392/1‑90 (MQ=255)
                       cATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGGTGGGTTTTTTGGTCGGGTTTGGATTGCATGAAGaatc                                                   >  1:458963/1‑87 (MQ=255)
                           aGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTc                                               <  1:595820/90‑1 (MQ=255)
                                                          aaGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGaagcga                <  1:246358/90‑1 (MQ=255)
                                                              ggCGGTGGTCGTGGTGGGTTTTTTGTTCGGGTTTGGGTCGCAAGAAGAACCGGGGGGGTCCCCCCACGAGACGATCAAAAaaacgtaaag            >  2:532004/1‑88 (MQ=255)
                                                              ggCGGTGGTCGTGATGGGTTTTTTGCTCGGGTTTGGATCGCATGAAGAACCGGGCGGCTCCCCCACCGAGACGAACAAAAaagcgaaaac            >  2:154338/1‑88 (MQ=255)
                                                              ggCGGTGGTCGGGGTGGGTTTTTTTCTCGGGTTTGGATCGCATGAAGATCCGGGCGGGGCCCCCACCGAGACGATCAAAAaaacgaaaac            >  2:685620/1‑88 (MQ=255)
                                                              ggCGGTGATCGTGGTGGATTTTTTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCCCCATCGAGACGATCAAAAaagcgaaaac            >  1:320667/1‑88 (MQ=255)
                                                              ggCGGTGATCGTGATGGATTTTTTGCTCGGGTTTTGGTCGCATGAAGAACCGGTCGGCTCCCCCCTCGAGACGAGCAAAAaagcgaaaga            >  2:480636/1‑89 (MQ=255)
                                                                       cGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATcgccgc   <  2:266231/90‑1 (MQ=255)
                                                                        gTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCt  <  2:634811/90‑1 (MQ=255)
                                                                                    |                                                                             
CCGATGATTGACCCCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCT  >  NZ_CP009273/333701‑333862

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