Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,514,1100GT53.8% 7.5 / 15.9 13T90K (ACA→AAA) hofPDNA utilization protein HofP
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (7/0);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGC  >  NZ_CP009273/3514023‑3514187
                                                                                        |                                                                             
gTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCCGTTTTTc                                                                             >  1:165558/1‑90 (MQ=255)
    cAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCgg                                                                         <  1:317009/90‑1 (MQ=255)
        aTCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGGTTTTTCCCggacc                                                                      >  2:140740/1‑87 (MQ=255)
         tCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTAccc                                                                    <  1:140740/90‑1 (MQ=255)
         tCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTAccc                                                                    <  2:59537/90‑1 (MQ=255)
                  gtGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTTTCAGTTTGTGCCTGTAACCCGCGATact                                                           >  2:7637/1‑88 (MQ=255)
                       gCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCCGTTTTTCCCCGTACCCCGCGTTAGTACGTc                                                      >  2:274798/1‑90 (MQ=255)
                               ttGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTTTCCCGGTACCCCGCGTTAGTACGTc                                                      >  1:278043/1‑82 (MQ=255)
                               ttGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTc                                                      <  2:278043/82‑1 (MQ=255)
                                    tCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTTTCCCCGTACCCCACGGTGGTACGCCTGGCGTCCGCTGt                                         >  2:210219/1‑90 (MQ=255)
                                                  ttGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTTTCCCCGTACCCCGCGTTTGTACGTCTGGGGGCCGCTGTAAAAATTTCCCGcc                           >  2:308424/1‑90 (MQ=255)
                                                                aCCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACAt             <  2:182100/90‑1 (MQ=255)
                                                                           ggTTCGCAGTTTTTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTtgctgc  >  1:259177/1‑90 (MQ=255)
                                                                                        |                                                                             
GTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAGTTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGTAAAATTGTCCAGCCGTTTTCCAGCACATCGTTTTGCTGC  >  NZ_CP009273/3514023‑3514187

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