Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,515,7360TG57.1% 8.3 / 13.8 14S125S (TCA→TCChofMDNA utilization protein HofM
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base G (8/0);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 3.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGT  >  NZ_CP009273/3515647‑3515816
                                                                                         |                                                                                
cTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgagt                                                                                  <  2:571775/90‑1 (MQ=255)
  cacaTGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgagtga                                                                                <  2:187022/90‑1 (MQ=255)
   acaTGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGCtgagtgag                                                                               >  1:262385/1‑90 (MQ=255)
            cGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGCTGAGGGAGTCTTccccg                                                                      >  2:147452/1‑86 (MQ=255)
                     tCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAg                                                             <  2:687107/90‑1 (MQ=255)
                     tCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAg                                                             <  2:539962/90‑1 (MQ=255)
                       tGCCCGACTAAGCCGCGTTGCCAGCACTTTCCTTTGCGAGGCAGCCACGCTATAGGCGGGACTGAGTTAGTCTTCGCTGTAATCGAAgcg                                                           <  2:605157/90‑1 (MQ=255)
                           agagTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGCTGAGGGAGGCTTTTCTGTAATCGAAGCGCAgg                                                       >  2:262304/1‑90 (MQ=255)
                                       gTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGGGGGTCCTTGCTGTTATTGAAGCGCAGGGAGGCCgggtac                                           >  2:527131/1‑86 (MQ=255)
                                                   ttgctttgcGCGGCAGTCACGTTATAGGCGGGACTGGGGGAGTCTTCTCTGTTATCGAAACGCAGGGGGGCCGGGTCCATATCCAGCTcg                               >  2:621435/1‑90 (MQ=255)
                                                   ttgctttgcGCGGCAGTCACGTTATAGGCGGGACTGAGGGAGTCTTCGCTGTAATCGAAACGCAGGGAGGCCGGGTCCATATCCAACTcc                               >  2:231093/1‑89 (MQ=255)
                                                      ctttgcGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCgg                            <  1:348395/90‑1 (MQ=255)
                                                               gCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTc                   <  1:262304/90‑1 (MQ=255)
                                                               gCAGTCACGTTATAGGCGGGACTGAGGGAGTCTTCCCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCCGCTCGCGGGCCATCGTc                   >  1:225492/1‑90 (MQ=255)
                                                                  gTCACGTTATAGGCGGGGCTGAGGGGGGCTTCGCTGTAATCGAAGCGCAGGGGGGCCGGGTCCATTTTCCGCTCGCGGGCCATCGGCCCt                >  1:596498/1‑90 (MQ=255)
                                                                         tataGGCGGGACTGAGGGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGcc         >  1:518630/1‑90 (MQ=255)
                                                                                gggACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGGCCGGATCCATATCCAGCTCCCGGGCCCACGTCCCCGACAACCAGGCCgg  >  2:645676/1‑89 (MQ=255)
                                                                                         |                                                                                
CTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATATCCAGCTCGCGGGCCATCGTCCCTGACAGCCAGGCCGT  >  NZ_CP009273/3515647‑3515816

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