Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273536,0140TG56.2% 3.4 / 12.7 16intergenic (+50/‑7)allB/ybbYallantoinase AllB/uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base G (9/0);  total (13/3)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAAC  >  NZ_CP009273/535931‑536092
                                                                                   |                                                                              
gAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGTGttt                                                                          >  2:89038/1‑90 (MQ=255)
                 ttAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGGTTAGGGGGGTTTTTGTTTAATTTTTCAGt                                                         >  2:148256/1‑90 (MQ=255)
                        tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGTTTATGTTCAATTTTTCAGTCAGCcgc                                                  >  1:364956/1‑90 (MQ=255)
                               tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGTGTTTTTGTTTAATTTTTCAGGCAGCCGCGAAAGcc                                           >  1:265571/1‑90 (MQ=255)
                                aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTTTGTTCAATTTTGCAGTCAGCCGCGAAAGcc                                           >  1:8954/1‑89 (MQ=255)
                                aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGcc                                           <  2:8954/89‑1 (MQ=255)
                                     ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCt                                          >  1:84812/1‑85 (MQ=255)
                                     ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCt                                          <  2:84812/85‑1 (MQ=255)
                                          cTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCCGTTTAGGGGTGGTTTTGTTTTATTTTTGAGTCAGCCCCGGAAGCCCGTTTTCTggg                                >  2:181961/1‑89 (MQ=255)
                                          cTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGGGGTTTTGTTTCATTTTTCAGTCAGCCGCGAAAACCTGTTATCAGGa                                >  1:343976/1‑90 (MQ=255)
                                          cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGTGGGTTTTTTTTTAATTTTTTAGGCGGCCGCGAAAGCCTGTTATCAggg                                >  2:261903/1‑89 (MQ=255)
                                          cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTTTGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa                                >  1:35611/1‑90 (MQ=255)
                                                     gTCCTTGCGGCGGGCATTCTCCCGTTTAGGGGTGTTTTTGTTCAATTTTGCAGTCCACCCCGGAAGACTGGTTTAAGGGTTTTAGAGGtt                     >  1:327862/1‑90 (MQ=255)
                                                      tCCTTGCGGCGGGCGTTCTTCGGGTTAGGGGGGTTTTTATTTTATTTTTCAGTCAGGCCCGGAAACCCGTTATCAAGATTTCAGCGGgtt                    >  2:330971/1‑90 (MQ=255)
                                                                gggCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTTTTTAttt          >  2:90485/1‑90 (MQ=255)
                                                                        tcCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAac  <  1:188483/90‑1 (MQ=255)
                                                                                   |                                                                              
GAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAAC  >  NZ_CP009273/535931‑536092

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