Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,344,2120AC54.5% 17.2 / 29.8 22S725A (TCA→GCA) arcBaerobic respiration two‑component sensor histidine kinase ArcB
Reads supporting (aligned to +/- strand):  ref base A (10/0);  new base C (0/12);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.55e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACA  >  NZ_CP009273/3344129‑3344297
                                                                                   |                                                                                     
ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgc                                                                                 >  1:129814/1‑90 (MQ=255)
ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgc                                                                                 >  2:358545/1‑90 (MQ=255)
   cccGACTTTTTCTCCCCGGGCCGGAAGGTAAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCcgccgc                                                                              <  1:342420/89‑1 (MQ=255)
          ttATCTCCCCGGGCGGGAGGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTa                                                                       <  1:280503/90‑1 (MQ=255)
          ttATCTCCCCGGGCGGGAAGGTCAGGAGACTAATTTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTa                                                                       <  1:383323/90‑1 (MQ=255)
               tccccAGGCCGAAGGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                  <  1:358545/88‑1 (MQ=255)
                 cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                >  2:93312/1‑90 (MQ=255)
                   cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGt                                                              >  2:189762/1‑90 (MQ=255)
                    aGGCGGAAAGGTCGGGAAACTAATTTTCCTCCCCCAGTTTTTGCAGAGGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGtc                                                             <  2:19406/90‑1 (MQ=255)
                          gAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct                                                       >  2:313148/1‑90 (MQ=255)
                          gAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct                                                       >  2:79167/1‑90 (MQ=255)
                                cgggAGACTGATTTTGCTGCCCCAGTTGTTGCAGGGGGCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                 <  1:499730/88‑1 (MQ=255)
                                   gagaCTAATTTTCCTGCCCCAGTTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTCTAATTTTATGTCCTTCCTCAACAATGc                                              <  2:129814/90‑1 (MQ=255)
                                    agaTTGATTTTCCTGCCCCAGTTTTTGCAGTGGGCCAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                             <  2:475761/90‑1 (MQ=255)
                                    agaCTGATTTTGCTGCCCCAGTTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                             <  1:115904/90‑1 (MQ=255)
                                       cTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                          >  2:459937/1‑90 (MQ=255)
                                              tGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATcc                                   >  2:249369/1‑90 (MQ=255)
                                                   aCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc                              <  1:313148/90‑1 (MQ=255)
                                                             tGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                    >  2:281697/1‑90 (MQ=255)
                                                              gaaGATGGCGAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                   <  2:219210/88‑1 (MQ=255)
                                                                      cGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc           <  1:284930/90‑1 (MQ=255)
                                                                               aCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACa  >  1:492408/1‑90 (MQ=255)
                                                                                   |                                                                                     
TTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACA  >  NZ_CP009273/3344129‑3344297

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