Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,379,0450TG64.7% 4.0 / 13.6 17G224G (GGT→GGGycjXYcjX family protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (11/0);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 2.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATTATCTCCACCACTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGCATACCAATGCC  >  NZ_CP009273/1378965‑1379132
                                                                                |                                                                                       
gATTATCTCCACCACTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCcgcg                                                                                >  2:292393/1‑90 (MQ=255)
          accacTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAAttct                                                                      >  2:357473/1‑90 (MQ=255)
           ccacTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGGCCCGCGCTGCAAttttt                                                                     >  2:354850/1‑90 (MQ=255)
           ccacTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGGCCCCCGCCGCCAttctt                                                                     >  1:451192/1‑90 (MQ=255)
              cTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGGCCCGCGCTGCCATTCTTccc                                                                  >  1:553969/1‑90 (MQ=255)
                gTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGt                                                                <  1:357473/90‑1 (MQ=255)
                gTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGGAATTCTTCCCgg                                                                >  2:52166/1‑89 (MQ=255)
                   aaGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCTGCAATTCTTTCCGTGGc                                                             >  1:97107/1‑90 (MQ=255)
                       gCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCCCCCGCGCTGCCATGTTTTCCGCTGCCGGa                                                         >  1:262904/1‑90 (MQ=255)
                        cAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCCATTCTTCCCCTGGCCgggt                                                        >  1:69472/1‑88 (MQ=255)
                               tGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGTCGCCCGCGCTGGAATTCTTCCCCTGGCCGGGAGTCGATa                                                 >  2:360219/1‑90 (MQ=255)
                                           aGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCCGCAATTTTTCCCCTGGCCGGGTGGTGGTACCTGGGGGGAgg                                     >  1:427703/1‑89 (MQ=255)
                                            gCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCTGCAATTCTTCCCCTTGCCGGGCGTCGATACCTCGGGCGAGTc                                    >  2:486896/1‑90 (MQ=255)
                                                         gTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAg                       <  2:451192/90‑1 (MQ=255)
                                                             tGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCg                   <  2:206824/90‑1 (MQ=255)
                                                                            cAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGGCCAAACTGGCGCAGGCCGATAAGCATACCAATg    >  1:240447/1‑90 (MQ=255)
                                                                              ggTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGCATACAAATGcc  >  2:113747/1‑90 (MQ=255)
                                                                                |                                                                                       
GATTATCTCCACCACTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGCATACCAATGCC  >  NZ_CP009273/1378965‑1379132

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