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

CAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATAT  >  NZ_CP009273/500850‑501022
                                                                                           |                                                                                   
cAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGt                                                                                     >  1:115637/1‑90 (MQ=255)
cAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGt                                                                                     >  2:117307/1‑90 (MQ=255)
cAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGt                                                                                     >  1:6627/1‑90 (MQ=255)
                   ttATTGGGCTGACACCCGATGCCACAGAAAAAATTGGTACCCCGAAATCCTTCACTGATATCGATTTTCGTAAGCCCGCCGATgaagcga                                                                  <  1:117307/90‑1 (MQ=255)
                   ttATTGGGCTGACAACCGATGACACAGAAAAAATTGGTAACCCGAAATACTTCACTGATATCAAATTTCGTAAGCCCGCCGATgaagcga                                                                  <  1:634900/90‑1 (MQ=255)
                   ttATTGGCCTAAAACCCAATAACAAAAAAAAAATTGGTACCCCGAAATCCTTCACTAATATCAATTTTCGAAAGCCCGCCGATgaagcga                                                                  <  2:115637/90‑1 (MQ=255)
                     ttggggCTGACACCCGATGACACACAAAAATTTGGAACCCGGGAATCCTTCACTAATTTGGATTTTGGAAGCCCCGCCGATgaagcgaag                                                                <  2:178505/87‑1 (MQ=255)
                             gACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAACTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGAt                                                        >  1:121463/1‑90 (MQ=255)
                                    gATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAg                                                 >  2:5622/1‑90 (MQ=255)
                                           caGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAac                                          >  1:465402/1‑90 (MQ=255)
                                           caGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAac                                          >  2:15697/1‑90 (MQ=255)
                                                       gTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGc                              >  2:260381/1‑90 (MQ=255)
                                                          ccccGAAATCCTTCACTAATATCAATTTTCGAAACCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAg                           <  1:509116/89‑1 (MQ=255)
                                                           cccGGAAATACTCCAATAATATAAAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAg                           <  1:380331/90‑1 (MQ=255)
                                                                       ttCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATc                >  2:186567/1‑90 (MQ=255)
                                                                       tccACTAATATCGATTTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATc                <  1:549064/88‑1 (MQ=255)
                                                                           ttaatatAGATTTTCGAACCCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcgg            <  2:233923/87‑1 (MQ=255)
                                                                           atGATATCAATTTTCGAACCCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcgg            <  1:609639/89‑1 (MQ=255)
                                                                             tatagCGATTTTCGAAACCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcggcg          <  2:258680/89‑1 (MQ=255)
                                                                              atatCAATTTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGa         <  1:108179/90‑1 (MQ=255)
                                                                                ctCGAATTAGGAACGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGAcc       <  1:427990/89‑1 (MQ=255)
                                                                                atCGAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGAcc       <  1:65046/90‑1 (MQ=255)
                                                                                atCGAATTTCGAAACCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGAcc       <  1:193187/90‑1 (MQ=255)
                                                                                     aTTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCatat  <  2:182183/90‑1 (MQ=255)
                                                                                           |                                                                                   
CAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATAT  >  NZ_CP009273/500850‑501022

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: