Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,780,5860AG58.1% 5.9 / 40.6 31S262P (TCG→CCG) ppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base G (10/8);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCAC  >  NZ_CP009273/1780453‑1780718
                                                                                                                                     |                                                                                                                                    
gCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCa                                                                                                                                 >  2:16069/1‑139 (MQ=255)
  cTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATg                                                                                                                               >  1:264081/1‑139 (MQ=255)
          ccAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTgcggc                                                                                                                        <  1:391056/138‑1 (MQ=255)
            aGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCa                                                                                                                     <  2:378763/139‑1 (MQ=255)
                         cttcGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCa                                                                                                                     >  2:267794/1‑126 (MQ=255)
                         cttcGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCa                                                                                                                     <  1:267794/126‑1 (MQ=255)
                         cttcGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGc                                                                                                         >  1:160947/1‑138 (MQ=255)
                                 gTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTcgccg                                                                                                >  2:278902/1‑139 (MQ=255)
                                 gTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTcgccg                                                                                                >  1:72005/1‑139 (MQ=255)
                                  tCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTcgccgc                                                                                               <  2:264081/139‑1 (MQ=255)
                                                  aCGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATg                                                                                <  2:308544/138‑1 (MQ=255)
                                                        ttCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTa                                                                          >  2:394059/1‑138 (MQ=255)
                                                           cTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTaaaa                                                                       <  2:60915/138‑1 (MQ=255)
                                                                      cTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTcgccgc                                                                                               <  1:438308/103‑1 (MQ=255)
                                                                      cTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTcgccgc                                                                                               >  2:438308/1‑103 (MQ=255)
                                                                        tCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGGGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTaa                                                         >  2:372664/1‑139 (MQ=255)
                                                                              ttAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCa                                                    <  2:251689/138‑1 (MQ=255)
                                                                                  ccTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTg                                               <  2:318239/139‑1 (MQ=255)
                                                                                      cTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCAc                                            >  1:240102/1‑138 (MQ=255)
                                                                                      cTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCAc                                            >  2:365406/1‑138 (MQ=255)
                                                                                         gCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGAc                                         <  2:485921/138‑1 (MQ=255)
                                                                                                  cTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCAcc                                <  1:365406/138‑1 (MQ=255)
                                                                                                     ggTCGGCGCGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAgg                            <  1:16069/139‑1 (MQ=255)
                                                                                                          gcgcGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCa                       >  1:350520/1‑139 (MQ=255)
                                                                                                          gcgcGTAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGAACATCTCAACAAGGccccc                       >  2:350897/1‑138 (MQ=255)
                                                                                                          gcgcGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCGCGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCa                       >  1:353791/1‑139 (MQ=255)
                                                                                                          gcgcGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCa                       >  2:4365/1‑139 (MQ=255)
                                                                                                               tAAACCATGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCgg                  <  1:473885/139‑1 (MQ=255)
                                                                                                                     aTGCGGATTTTTTTCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGa            <  2:72005/139‑1 (MQ=255)
                                                                                                                     aTGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGa            >  2:100690/1‑139 (MQ=255)
                                                                                                                                tttCGGCCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAAcaccac  <  1:278902/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GCCTGTTTTGCCAGTTCCTGCACTTCTTCGTTGGTCAGCGAGAAGATGTCACGCTGTTCCTGCGGTACGTCTTCGATTTTAACCTGCTTGCCGTGCTCCTGGGTCGGCGCGTAAACCATGCGGATTTTTTTCGACCCCATGGTGCGGCGCACGATAGCCGGGCGATTCGCCGCCAGTGTCGGTTTATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCAC  >  NZ_CP009273/1780453‑1780718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑