Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,322,4130TG60.0% 7.1 / 13.8 15T68P (ACC→CCC) cobOcob(I)yrinic acid a,c‑diamide adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (3/6);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCATGTTGCCAGACTTCGCGGCAGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGTGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCC  >  NZ_CP009273/1322279‑1322546
                                                                                                                                      |                                                                                                                                     
tGCATGTTGCCAGACTTCGCGGCAGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGTGCCt                                                                                                                                   >  2:259356/1‑139 (MQ=255)
tGCATGTTGCCAGACTTCGCGGCAGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGGGc                                                                                                                                     <  1:307672/137‑1 (MQ=255)
                 cgcgGCAGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGTGCCTTTAATAAACTGCacgac                                                                                                                  >  1:375030/1‑139 (MQ=255)
                      cAGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGGGCCTTTAATAAACTGCACGACGCCTa                                                                                                             <  1:505625/139‑1 (MQ=255)
                       aGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGGGCCTTTAATAAACTGCACGACGCCTAc                                                                                                            >  2:264181/1‑139 (MQ=255)
                                         cgcgGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGGGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAAc                                                                                          >  2:83497/1‑139 (MQ=255)
                                                  gtgtATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGGGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCgtg                                                                                  <  2:437492/138‑1 (MQ=255)
                                                                                          aaCGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGGGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTg                                         >  2:409844/1‑139 (MQ=255)
                                                                                           aCGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGTGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGc                                        >  1:297562/1‑139 (MQ=255)
                                                                                                   tgGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGTGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAg                                <  1:259356/139‑1 (MQ=255)
                                                                                                                    gTTCGCCATTAGGCCAGGGGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACcgcg               <  1:241976/139‑1 (MQ=255)
                                                                                                                           aTTAGGCCAGGTGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATcc        <  2:297562/139‑1 (MQ=255)
                                                                                                                                 ccAGGTGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATcctgggcc  <  1:46448/139‑1 (MQ=255)
                                                                                                                                 ccAGGGGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATcctgggcc  <  1:264181/139‑1 (MQ=255)
                                                                                                                                 ccAGGGGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATcctgggcc  <  1:409844/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TGCATGTTGCCAGACTTCGCGGCAGGCGGCGGTATCAGACTCGCGGTTTTGTGTATCCCAGGTAAAGCCCGTTGCCATCACCTGAAACTCAACGCCATGTGGCTCCAGCAGATTGCGTTCGCCATTAGGCCAGGTGCCTTTAATAAACTGCACGACGCCTACTTTTTTTCCGTGACCAACTGCGCGTGTTGCCGTACCAAATGCCGCGGTGGTTTTGCCTTTTCCATTGCCGGTAAAGACGATGATAATACCGCGTTCATCCTGGGCC  >  NZ_CP009273/1322279‑1322546

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