Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,502,0210TA53.8% 3.5 / 16.9 13E384D (GAA→GATyhfWphosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base A (7/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
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.

CATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA  >  NZ_CP009273/3501947‑3502103
                                                                          |                                                                                  
cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATTACACACCGGAGcc                                                                     >  2:22734/1‑90 (MQ=255)
       tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                              <  2:57621/90‑1 (MQ=255)
        gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACACCGGAGCCCCCCCATc                                                             >  1:131194/1‑90 (MQ=255)
        gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATTACACACCGGAGCCCCCCCATc                                                             >  1:71759/1‑90 (MQ=255)
              tGCAAACCGGAGGAAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGGAGCCCCCCCATCAGAAAg                                                       >  2:168439/1‑90 (MQ=255)
                cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAAAACACGGGAGCCCCCCCCACAGaaaagg                                                     >  1:119308/1‑87 (MQ=255)
                                        agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                             <  1:141360/90‑1 (MQ=255)
                                         gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                            <  2:288011/90‑1 (MQ=255)
                                           gACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAAAGGGGGGGGCACACCCAGCTGCGGAGCGa                          >  1:69901/1‑90 (MQ=255)
                                           gACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCCCCTCAGAAAACGGGGGGGGCACACCCAGCCGCGGAGCGa                          >  1:23221/1‑90 (MQ=255)
                                            aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                         <  2:131194/90‑1 (MQ=255)
                                                    gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                 <  2:71759/90‑1 (MQ=255)
                                                                   ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa  <  2:64325/90‑1 (MQ=255)
                                                                          |                                                                                  
CATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA  >  NZ_CP009273/3501947‑3502103

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