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

TCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTA  >  NZ_CP009273/3126583‑3126733
                                                                              |                                                                         
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                <  1:341688/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                <  1:770279/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                <  1:883784/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                <  2:37086/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                <  2:736504/90‑1 (MQ=255)
tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCTcc                                                                >  2:11791/1‑90 (MQ=255)
  cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAg                                                             <  2:403986/90‑1 (MQ=255)
      aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                          <  1:898409/74‑1 (MQ=255)
      aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                          >  2:898409/1‑74 (MQ=255)
      aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGGTCa                                                         >  1:905157/1‑90 (MQ=255)
                              gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                          <  2:169395/50‑1 (MQ=255)
                              gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg                                                                          >  1:169395/1‑50 (MQ=255)
                                 gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAATTCAAGGAACTTTTGGCCCTTTAATGGCaat                              >  2:402965/1‑88 (MQ=255)
                                     gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTTAATGTTAGTTggg                          >  2:266378/1‑88 (MQ=255)
                                          gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCCCCTTTAAAGGGCTTTTTGCCCTTTAATGGCAGTTgggggggg                     >  2:969391/1‑86 (MQ=255)
                                           cgcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCCCCCACGTTCAAGAACTTTTTTCCCTTGTATGTGTGTTGTGGCGTaa                    >  2:939875/1‑90 (MQ=255)
                                           cgcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGTTTAAGGGGCTTTTTGGCCTTTTATGTTAGTTGGGGCGtga                    >  2:978881/1‑88 (MQ=255)
                                                       cTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGGTCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGGGAgttttt        >  2:17093/1‑90 (MQ=255)
                                                             gAATTTCATCGGCGGGGGGGTGCCCATCCCGGTTAAGGGTCTTTTCGCCCTTTGATGTCAGTTGGGGGGGAAAGGCGGTTTTTTTCtttt  >  2:832951/1‑87 (MQ=255)
                                                                              |                                                                         
TCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTA  >  NZ_CP009273/3126583‑3126733

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