Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,502,0090CG56.2% 9.9 / 19.7 16A388A (GCG→GCCyhfWphosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base G (9/0);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTG  >  NZ_CP009273/3501926‑3502092
                                                                                   |                                                                                    
gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCgg                                                                               >  2:139358/1‑90 (MQ=255)
       gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGgaaaaat                                                                        >  2:207715/1‑88 (MQ=255)
        tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAAtt                                                                        >  2:2224/1‑90 (MQ=255)
            aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacac                                                                   <  2:95928/90‑1 (MQ=255)
                     cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAACACACCGGAGcc                                                          >  2:158007/1‑90 (MQ=255)
                     cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATAACACACAGGAGcc                                                          >  1:106393/1‑90 (MQ=255)
                          ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATTACACAACGGAGccccccc                                                     >  1:50246/1‑88 (MQ=255)
                           gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAAACACACACCGGAGCCCCccca                                                    >  1:167782/1‑90 (MQ=255)
                             gtCTCCTGCAAACCGGAGGGAGGAAAAAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATTACACACGGGAGCCCCcccacc                                                  >  1:311257/1‑88 (MQ=255)
                                   tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGCGGGAAAAATAAAACACCGTAACCCCCGCACCAAAAAg                                            >  1:166258/1‑90 (MQ=255)
                                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                          <  1:224419/90‑1 (MQ=255)
                                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                        <  1:341165/90‑1 (MQ=255)
                                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                        <  2:335935/90‑1 (MQ=255)
                                                 aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCacac                              <  2:167782/90‑1 (MQ=255)
                                                      agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTTACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                         <  1:225745/90‑1 (MQ=255)
                                                             agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                  <  2:50246/90‑1 (MQ=255)
                                                                             cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg  <  1:150950/90‑1 (MQ=255)
                                                                                   |                                                                                    
GTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTG  >  NZ_CP009273/3501926‑3502092

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