Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,794,4960TG54.5% 9.3 / 10.6 12V58G (GTG→GGG) nrdIclass Ib ribonucleoside‑diphosphate reductase assembly flavoprotein NrdI
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base G (6/0);  total (8/4)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACTTTGGTG  >  NZ_CP009273/2794409‑2794585
                                                                                       |                                                                                         
cGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTggggcggggc                                                                                         >  2:186312/1‑90 (MQ=255)
 gCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGcc                                                                                        <  1:389756/90‑1 (MQ=255)
     aaTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGAGGCGGGACGGCTGGCGCCGCGCCACgc                                                                                    >  1:340829/1‑89 (MQ=255)
                                              gATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCCCGACAGGGAATTCGCTTTTTAAACGACGAGCACAACCGGGCGtt                                           >  2:85898/1‑90 (MQ=255)
                                                  gTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGGCCCGACAGGGAATTTGCCTTTTTAACCACGAGCACAACCGGGCGGTGCtt                                       >  2:491148/1‑90 (MQ=255)
                                                    gCCCTCTTACGGCGGCGGCGGTACGGCTGGGGCGGGGCCCCGGCAGGGAAATCGCTTTTTTAACGGCGGGCACAACCGGGCGGTTCTTcg                                     >  2:389557/1‑90 (MQ=255)
                                                    gCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGGGCCACGACAGGGAATTCGCTTTTTAAACGACGAGCACAACCGGGCGGTGCTTcg                                     >  2:130850/1‑90 (MQ=255)
                                                     cccTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgc                                    <  1:274934/90‑1 (MQ=255)
                                                      ccTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGCAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg                                   <  1:337453/90‑1 (MQ=255)
                                                        tctTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGGTGCTTCgcggc                                 >  1:489533/1‑90 (MQ=255)
                                                                 ggcggcggGACGGCTGGGGCGGGGCCACGACAGGGAATTCGCTTTTTTAACGACGAGCACAACCGGGGGGTGCTTCGCGGGGGTATTGCt                        >  2:204998/1‑90 (MQ=255)
                                                                   cggcggTACGGCTGGGGCGGGGGCACGACAGGGAATTCGCTTTTTAAAAGACGGGCACAACCGGGGGTTTCTTCGCGGCGGTATttcttc                      >  2:84175/1‑90 (MQ=255)
                                                                                  cgcgGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACtt       <  2:417915/90‑1 (MQ=255)
                                                                                       tGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACTTTGGTg  <  1:385540/90‑1 (MQ=255)
                                                                                       |                                                                                         
CGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTGCTTCTGGTAATCGCAACTTTGGTG  >  NZ_CP009273/2794409‑2794585

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