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

GCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTG  >  NZ_CP009273/1927573‑1927741
                                                                                        |                                                                                
gCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGAt                                                                                 >  2:328255/1‑90 (MQ=255)
gCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGAt                                                                                 >  2:19737/1‑90 (MQ=255)
         aGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGcc                                                                        <  1:328255/90‑1 (MQ=255)
               ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                  >  2:577574/1‑90 (MQ=255)
               ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc                                                                  >  2:124878/1‑90 (MQ=255)
                   taaaCATTTTTGCCCCCCCATGAGGAGAAAAAGGTTTTTCGAAGGAAGCGATCACTTTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGcc                                                              <  1:19737/88‑1 (MQ=255)
                             tgccccACCATGATGAGAAAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTGGAGTGTTTTTTCCCCCCCTCCCCGCCAGTCCAGTTc                                                    <  2:402967/87‑1 (MQ=255)
                                      atgatgAGAAAAAGTTTGTTCGAAGGAAGCGATCACTTTGTTGTCGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCattatt                                           <  1:242047/90‑1 (MQ=255)
                                        gatgaGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCa                                         >  2:582614/1‑90 (MQ=255)
                                             gagaAGGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                    <  1:189686/90‑1 (MQ=255)
                                             gagaAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTTTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt                                    <  1:43050/90‑1 (MQ=255)
                                                    tttgtcCAAAGGACGCGACCACCTTGCTCTGAGGTTTTTTTCCCCCCCCTGCCCGCCAGTCCAGTTCACCATTCTTCAGCCATGGTTCaa                             <  1:452086/89‑1 (MQ=255)
                                                    gttgttCGAAGGAAGCGATCACTTTGCTGTGGAGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa                             <  2:582761/90‑1 (MQ=255)
                                                       gttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAggg                          >  1:180937/1‑90 (MQ=255)
                                                        ttCGAAGGAAGCGACCACTTTGTTGTGGGTTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGt                         <  1:280446/90‑1 (MQ=255)
                                                        tccaaaggaagCGATCACATGGTGGTCGAGTTTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGt                         <  1:297380/86‑1 (MQ=255)
                                                                   cGATCACTTTGCTGTCGATTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagaca              <  1:124878/90‑1 (MQ=255)
                                                                               tGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTg  >  1:490907/1‑90 (MQ=255)
                                                                                        |                                                                                
GCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTG  >  NZ_CP009273/1927573‑1927741

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