Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,695,5600TC55.0% 20.4 / 27.3 20I116T (ATC→ACC) malYbifunctional maltose regulon transcriptional repressor/cystathionine beta‑lyase MalY
Reads supporting (aligned to +/- strand):  ref base T (9/0);  new base C (0/11);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 5.95e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAA  >  NZ_CP009273/1695480‑1695638
                                                                                |                                                                              
ggtggtGTATGGCCCTTCTGTCATCTTTATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACAcc                                                                       >  1:727898/1‑90 (MQ=255)
 gtggtgTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACAccc                                                                      >  2:815153/1‑90 (MQ=255)
            cccTTCTCTCATCTATTTGGTTTAAGAACTGATTCGTCGGTGGTCAAAAACGGGTGAAGGGGTGGTGACCCACACACCCGCCTATGACGc                                                           <  2:620775/90‑1 (MQ=255)
             ccTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCa                                                          >  2:331915/1‑90 (MQ=255)
                       tttatgTGTTTTAAGAACTGTTTCGTGAGGGGTCTAAAACGGGGGAGGGGGTGGTGACCCCCACCCCCGCCTATGACGCATTTTACAAgg                                                <  1:802661/88‑1 (MQ=255)
                        ttataTGTTTTAAGAACTGTTTCGTCAGGGGTCTGAAACAGGGGAAGGCGTGGTGACCCACACACCCGCCTATGACGCATTTTACAAGGc                                               <  2:797217/89‑1 (MQ=255)
                          atatGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCa                                             >  1:628806/1‑90 (MQ=255)
                               gTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGaa                                        >  2:545895/1‑90 (MQ=255)
                               gTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGaa                                        >  2:461259/1‑90 (MQ=255)
                                 ttCAGAACTGTTTCGTCGGTGGTCTAAAAAGGGGGAGGGGGTGGTGACCCACCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAgg                                      <  1:461259/90‑1 (MQ=255)
                                      aaCTGATTCGTCAGTGGTCTGAAACGGGTGAAGGGGGGGTGACCCACACCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAAcc                                 <  2:376709/90‑1 (MQ=255)
                                          gATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAgcg                             >  1:738468/1‑90 (MQ=255)
                                           tttcttcAGTGGTCTAAAACAGGGGAAGGCGGGGTGACCCACACCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAgcgc                            <  2:602050/89‑1 (MQ=255)
                                             tcgtcgGTGGTCTGAAACGGGTGAAGGCGTGGTGACCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGcac                          <  2:284081/90‑1 (MQ=255)
                                             ccGTCAGTGGTGTGAAACGGGGGAGGGGGTGGGGACCCACCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGcac                          <  1:59836/89‑1 (MQ=255)
                                                        cTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCg               >  1:51216/1‑90 (MQ=255)
                                                         taaaaCAGGGGAGGGGGTGGTGACCCCCCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGt              <  1:277447/88‑1 (MQ=255)
                                                         tGAAAGAGGGGAGGGGGGGGTGACCCACCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGt              <  1:304808/90‑1 (MQ=255)
                                                          aaaaCGGGAGAGGGGGTGGTGACCCCCCCCCCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGtt             <  1:815153/89‑1 (MQ=255)
                                                                     aGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGaa  >  2:501579/1‑90 (MQ=255)
                                                                                |                                                                              
GGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAACAGGTGAAGGCGTGGTGATCCACACACCCGCCTATGACGCATTTTACAAGGCCATTGAAGGTAACCAGCGCACAGTAATGCCCGTTGCTTTAGAGAA  >  NZ_CP009273/1695480‑1695638

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