Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273536,0330GT53.8% 7.0 / 20.7 13A5S (GCA→TCA) ybbYuracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (7/0);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCC  >  NZ_CP009273/535955‑536106
                                                                              |                                                                         
tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTAAGGGGGGTTTATTTTCCATTTTTCAGGCAGCcgc                                                                >  1:197731/1‑90 (MQ=255)
       tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGcc                                                         <  1:250733/90‑1 (MQ=255)
        aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCCGGGTAAGGGGGGGTTTTGTTCAATTTTTCCGGCAGCCGCGGAAGCCt                                                        >  2:250733/1‑90 (MQ=255)
             ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGGGGGTTTTTGTTTAATTTTTCAGTCAGCCGCGAAAGCCTGtttt                                                   >  2:248771/1‑88 (MQ=255)
                  cTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCGGGTAAGGGGGGTTTTTTTTTTATTTTTTAGTGAGCCCCGGAAACCTGTTTTCAggg                                              >  1:356255/1‑89 (MQ=255)
                  cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTTAGGGGGGTTTTTGTTTAATTTTTTAGTCAGCCGCGGAAACCCGTTTTTAggg                                              >  1:593217/1‑89 (MQ=255)
                         gCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGt                                       <  1:297975/90‑1 (MQ=255)
                            cGTCCTTGCGGCGGGCATTATTCGGTTAAGGGGGGTTTTTGTTTAATTTTTCAGTCGGCCGGGGAAGCCTGTTTTCCGGGTTTTTGTGGt                                    >  1:565428/1‑90 (MQ=255)
                                   gcggcgGGCATTCTCCGGTTAAGGGGGGTTTATGTTTAATTTTTCAGGCAGCCGCGAAAGCCTGTTATCAGGATTTTAGTGGGTtttttt                             >  2:297975/1‑90 (MQ=255)
                                                  cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac              <  1:88114/90‑1 (MQ=255)
                                                        aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTtgt        >  2:239426/1‑90 (MQ=255)
                                                            gtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt    <  2:565428/90‑1 (MQ=255)
                                                              gtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTtcc  >  1:9215/1‑90 (MQ=255)
                                                                              |                                                                         
TCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCC  >  NZ_CP009273/535955‑536106

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