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

TCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT  >  NZ_CP009273/535937‑536103
                                                                             |                                                                                         
tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTc                                                                               <  2:526310/90‑1 (MQ=255)
                   cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGGGTTTTTGTTTAATTTTTCAGTCAGCcgcg                                                            >  2:79673/1‑90 (MQ=255)
                      cagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGaaa                                                         <  1:79673/90‑1 (MQ=255)
                          aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGTGTTTATGTTTAATTTTTCAGGGAGTCGCGAAAGccg                                                     >  2:354395/1‑89 (MQ=255)
                          aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTTTGTTCAATTTTTCAGGCAGCCGCGAAAGCCt                                                     >  2:101145/1‑90 (MQ=255)
                               ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGGGGTTATGTTCAATTTTGCAGGCAGCCCCGAAAGCCTGTTAt                                                >  1:705162/1‑90 (MQ=255)
                                gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATc                                               <  2:466323/90‑1 (MQ=255)
                                    cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa                                           <  1:409421/90‑1 (MQ=255)
                                    cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAggg                                           >  1:87167/1‑89 (MQ=255)
                                    cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTTTGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTTTCAGGa                                           >  1:454665/1‑90 (MQ=255)
                                         aTGCCCGTCCTTGCGGCGGGGATTTTCCGGGTAAGGGGGGGTTTTGTTTAATTTTTTGGGCCGGCGCGGAAGGCTGTTTTCAGGGtttta                                      >  2:410733/1‑88 (MQ=255)
                                         aTGCCCGTCCTTGCGGCGGGCATTTTCCCGGTAAGGGGGGGTTTTGTTCAATTTTTCAGGCAGCCGCGAAAGCCTGTTATCAGGAtttta                                      >  2:572391/1‑88 (MQ=255)
                                           gCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTTTTTTTAATTTTTCAGTCAGCCGCGGAAACCTGTTTTTAGGGTTTTAGt                                    >  2:520931/1‑90 (MQ=255)
                                              cGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTTAGTggg                                 >  1:183067/1‑89 (MQ=255)
                                               gTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTTTTTAATTTTTTCGTCAGCCGCGAAAACCCGTTATCAGGGTTTTAGTgggt                                >  2:597687/1‑88 (MQ=255)
                                                          ggggATTCTCCGGGTTAGGGGTGTTTTTGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATTAGGGTTTTAGTGGGtttttttttttt                     >  2:441097/1‑90 (MQ=255)
                                                                    cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac           <  2:141871/90‑1 (MQ=255)
                                                                    cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac           >  2:349075/1‑90 (MQ=255)
                                                                    cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac           <  1:379326/90‑1 (MQ=255)
                                                                             tgtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt  >  2:48393/1‑90 (MQ=255)
                                                                             |                                                                                         
TCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT  >  NZ_CP009273/535937‑536103

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