Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,927,6610AT52.4% 7.3 / 29.6 21S401T (TCA→ACA) eddphosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base T (0/11);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTG  >  NZ_CP009273/1927581‑1927741
                                                                                 |                                                                                
cAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGc                                                                         <  2:154327/90‑1 (MQ=255)
       ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                  >  2:476149/1‑90 (MQ=255)
       ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                  >  1:244728/1‑90 (MQ=255)
         gCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGAGCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCg                                                                >  1:389565/1‑90 (MQ=255)
           ttAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTAGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGcc                                                              >  1:531382/1‑90 (MQ=255)
            tAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGcca                                                             >  2:316954/1‑90 (MQ=255)
                   tttGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGccagtccctg                                                      >  1:398410/1‑87 (MQ=255)
                                   gagagaAGGGTTGTTCGAGGGAGGCGATCACTTTGCTGTCGAGTGTTTTTCCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGcc                                      <  2:259349/90‑1 (MQ=255)
                                     gaaaaaGTTTTTTCGAAGGAAGCGATCACTTTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                    <  2:79106/87‑1 (MQ=255)
                                     gaaaaaGTTTGTTCGAAGAACGCAACCACTTTGCTGTGGATTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                    <  2:601914/87‑1 (MQ=255)
                                     gaaaaaGTTTGTCCGAAGGAAGCGATCACTTTGTTGTGGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                    <  1:310668/87‑1 (MQ=255)
                                      gaaaaGGTTTTTTCGAAGGAAGCGATCACTTTGCTGTCGATTTTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                    <  2:244728/88‑1 (MQ=255)
                                           ggTTGTTCAAGGGAAGCGATCCCTTTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa                              <  2:115965/90‑1 (MQ=255)
                                            tttgttCAAAGAAAGCGATCCTATTGCTGTCGAGTGTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa                             <  1:558223/89‑1 (MQ=255)
                                               gttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAggg                          >  1:269873/1‑90 (MQ=255)
                                                 ttCGAAGGAAGCAACCACTTTGCTGTCGGGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGt                         <  2:269873/90‑1 (MQ=255)
                                                           gCGATCCCTTTGTTGTGGAGTTTTTTTTCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac               <  1:476149/90‑1 (MQ=255)
                                                                  ctttgttgTCGCGTTTTTTTCCCGCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCaa        <  1:520131/88‑1 (MQ=255)
                                                                        tGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTCCAAGGGTATAACGAGACAGACCAAAACCTg  >  1:34491/1‑90 (MQ=255)
                                                                            gAGTGTTTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac               <  1:81954/73‑1 (MQ=255)
                                                                            gAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac               >  2:81954/1‑73 (MQ=255)
                                                                                 |                                                                                
CAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTG  >  NZ_CP009273/1927581‑1927741

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