Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,126,6610TG56.2% 3.8 / 18.9 16H219P (CAC→CCC) yghSdTMP kinase
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (9/0);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGT  >  NZ_CP009273/3126583‑3126748
                                                                              |                                                                                       
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                              <  1:496240/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                              <  2:275924/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                              <  2:391840/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGGCCCTcc                                                                              >  1:283345/1‑90 (MQ=255)
      aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCCCCCCGTTcc                                                                        >  1:64536/1‑89 (MQ=255)
          aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                    <  1:696548/90‑1 (MQ=255)
           ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGcc                                                                                   <  2:74111/74‑1 (MQ=255)
           ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGcc                                                                                   >  1:74111/1‑74 (MQ=255)
                gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAGGGGCtttt                                                              >  2:312926/1‑89 (MQ=255)
                         aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGTTCAAGGATCTTTTCGCCCttta                                                     >  2:730800/1‑88 (MQ=255)
                            gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCACCCAGGTCAAAGAGCTTTGGGCCATTGTATg                                                  >  2:646897/1‑90 (MQ=255)
                              gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCCCCCCTTCAAGGGGCTTTTGGCCCTTTGATGTc                                                >  1:418216/1‑90 (MQ=255)
                                 gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCTCCCGTTTAAAGATTTTTTCGCCCTTTAATGGCAGt                                             >  2:696111/1‑90 (MQ=255)
                                                                       ggcggGGGGGCGCCCCTCCCGTTCAAAGATCTTTGCGCCCTTGAATGTCAGGTGCGGCGTAACGGCGGGTTTTTCCTGTAACGCCGCCAg       >  2:636941/1‑90 (MQ=255)
                                                                           gggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGACAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtgg   <  2:346250/90‑1 (MQ=255)
                                                                            ggTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtggt  <  2:418216/90‑1 (MQ=255)
                                                                              |                                                                                       
TCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGT  >  NZ_CP009273/3126583‑3126748

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