Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,434,5290AC70.6% 27.6 / 10.0 17F186V (TTT→GTT) fabBbeta‑ketoacyl‑ACP synthase I
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base C (7/5);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAG  >  NZ_CP009273/2434442‑2434614
                                                                                       |                                                                                     
gggTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAcca                                                                                     >  1:247750/1‑90 (MQ=255)
                 gtagACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGccc                                                                    <  2:97090/90‑1 (MQ=255)
                    gACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGt                                                                 <  2:178057/90‑1 (MQ=255)
                    gACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGt                                                                 >  1:178057/1‑90 (MQ=255)
                     aCAGCGCCCCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGtt                                                                <  1:123179/90‑1 (MQ=255)
                      cAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTg                                                               >  1:276040/1‑90 (MQ=255)
                      cAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTg                                                               >  1:163213/1‑90 (MQ=255)
                        gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                             >  1:207077/1‑90 (MQ=255)
                        gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                             >  2:243654/1‑90 (MQ=255)
                          gcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATc                                                           <  2:47913/90‑1 (MQ=255)
                                     gTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGc                                                >  2:177293/1‑90 (MQ=255)
                                                aaGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATAc                                     <  2:247750/90‑1 (MQ=255)
                                                          ccAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGa                           <  1:294713/90‑1 (MQ=255)
                                                                     ctTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACAcg                >  2:113210/1‑90 (MQ=255)
                                                                         gccgccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGa            <  2:40903/90‑1 (MQ=255)
                                                                          ccgccgccAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAg           <  2:275767/90‑1 (MQ=255)
                                                                            gccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTAACGATACAGTGTGCGGAAGTCGCACACGCGGAGCt         >  2:308855/1‑90 (MQ=255)
                                                                                   gCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGag  <  2:252435/90‑1 (MQ=255)
                                                                                       |                                                                                     
GGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGATGGAG  >  NZ_CP009273/2434442‑2434614

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