Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273500,9390TA53.3% 9.4 / 20.0 15R190R (CGT→CGAushAbifunctional UDP‑sugar hydrolase/5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (6/1);  new base A (0/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGC  >  NZ_CP009273/500850‑501027
                                                                                         |                                                                                        
cAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGt                                                                                          >  2:469910/1‑90 (MQ=255)
cAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGt                                                                                          >  2:10956/1‑90 (MQ=255)
                  gTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATgaagcg                                                                        >  2:88191/1‑90 (MQ=255)
                             gACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGAt                                                             >  2:232483/1‑90 (MQ=255)
                                   cGATGACACAGAAAAAATTGGTACCCCGGAATCCTTCACTGATATCGAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGa                                                       <  1:196236/90‑1 (MQ=255)
                                    gATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAg                                                      >  2:153976/1‑90 (MQ=255)
                                     aTGACAAAGCAAAATGTGGCACCCCGGAATCCTAAACGAATTTCTTTTTCCGAACCCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGc                                                     <  2:320024/90‑1 (MQ=255)
                                       gACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTg                                                   >  1:221001/1‑90 (MQ=255)
                                          acaGAAAAAATTGGAACCCCGAAATCCTTCACTGATACCAAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCaa                                                <  1:65794/90‑1 (MQ=255)
                                                   aTTGGAACCCCGAAATACTCCACTGATATCAAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGaa                                       <  1:208483/90‑1 (MQ=255)
                                                    ttGGAACCCCGAAATACTTCACTAATATAGATTTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGaaa                                      <  1:133765/90‑1 (MQ=255)
                                                            ccGGAATTCTTCATTGTTACCATTTTTCAAAACCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGAc                              <  2:408405/90‑1 (MQ=255)
                                                                    cTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATtattat                      >  1:341340/1‑90 (MQ=255)
                                                                         ttaatatCAATTTTCGAAACCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcgg                 <  1:415745/87‑1 (MQ=255)
                                                                         cTGAAATCAATTTTCGAACGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcgg                 <  1:415746/90‑1 (MQ=255)
                                                                                     ttCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATgg     >  1:17986/1‑90 (MQ=255)
                                                                                        gTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGc  <  1:232483/90‑1 (MQ=255)
                                                                                         |                                                                                        
CAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGC  >  NZ_CP009273/500850‑501027

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