Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273108,7890TC57.1% 8.0 / 23.5 21S100G (AGC→GGC) zapDcell division protein ZapD
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base C (0/12);  total (7/14)
Fisher's exact test for biased strand distribution p-value = 3.10e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCA  >  NZ_CP009273/108706‑108858
                                                                                   |                                                                     
gATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgc                                                                 <  2:713788/90‑1 (MQ=255)
   gCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc                                                              >  1:705582/1‑90 (MQ=255)
        gccGCTGACGCACCAGAGCAATCAAACGCTCTCCCCGCAAAAATTCCCCGATACGCGGCGGGAAAATTAATCCCCCCCCCGCCGCTTTTa                                                         <  1:468241/88‑1 (MQ=255)
                      agagCAATCAAACGTTCTCCCCGCAGAAATTGCCCGATACGGGGGGGGAAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAAt                                           <  1:502853/90‑1 (MQ=255)
                      agagCAATCAAACGATCTCCCCGAAAAAATTCCCCGATCCGCGGCGCGGAAATTAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAAt                                           <  1:572031/90‑1 (MQ=255)
                           aaTCAAACGATCTTCACGCAAAAATTCCCCGATCCGCGGCGCGGAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTc                                      <  2:533560/90‑1 (MQ=255)
                              cAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAAt                                   >  2:668104/1‑90 (MQ=255)
                                aaCGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc                                 <  2:625687/90‑1 (MQ=255)
                                    aTCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGct                             >  1:765309/1‑90 (MQ=255)
                                     tCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGctc                            <  2:123168/90‑1 (MQ=255)
                                      cTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGctct                           >  1:264527/1‑90 (MQ=255)
                                        tCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTgg                         >  1:738538/1‑90 (MQ=255)
                                          aCGCAGAAATTGCCGGATCCGCGGCGCGAAAATTAATACCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTc                       <  2:476178/90‑1 (MQ=255)
                                               aaaaTTCCCCGATACGCGGCGCGAAAATAAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGc                  <  1:708471/89‑1 (MQ=255)
                                                    tccccGATCCGCGGCGCGAAAATTAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGc             <  2:725622/88‑1 (MQ=255)
                                                    tccccGATCCGCGGCGCGAAAATTAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGc             <  1:241594/88‑1 (MQ=255)
                                                         gATACGGGGCGAGAAAATTAACCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGcc        <  1:447805/90‑1 (MQ=255)
                                                              ggggcgggaAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc   <  1:331528/88‑1 (MQ=255)
                                                               cggcgcggAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCa  >  2:68696/1‑90 (MQ=255)
                                                                ggcgcggAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc   <  2:143394/88‑1 (MQ=255)
                                                                ggcgcggAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc   >  1:143394/1‑88 (MQ=255)
                                                                                   |                                                                     
GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCA  >  NZ_CP009273/108706‑108858

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