Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273688,3940AC56.2% 8.2 / 10.2 16G147G (GGT→GGGybeZPhoH family protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (0/9);  total (3/13)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCCGGACGAGTCAGCAGAATACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAAT  >  NZ_CP009273/688312‑688483
                                                                                  |                                                                                         
cTGCCGGAGGAGTCAGAAAAATACGGCGAATTCCCTGCCGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtacccgtac                                                                                    <  2:408057/90‑1 (MQ=255)
cTGCCGGACGAGTCAGCAGAATACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtacccgtac                                                                                    >  1:425154/1‑90 (MQ=255)
  gCCGGACGAGTCAGCAGAATACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtacccgtaccc                                                                                  >  1:292547/1‑90 (MQ=255)
  gCCGGACGAGTCAGCAGAATACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtacccgtaccc                                                                                  >  2:164005/1‑90 (MQ=255)
           gTCAGCAAAATACGGCGATTTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCa                                                                         <  2:565106/90‑1 (MQ=255)
            tCAGCAGAATACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCaa                                                                        <  2:32045/90‑1 (MQ=255)
                   aacaCGGCGATTTTCCTGGCCCCCCGGGGCAAAAACTCCCGCACCCCTTCCCAGGTAGGTTTTCCCCGCACCCGCCGGGCCAACGCCGaa                                                                 <  1:111729/87‑1 (MQ=255)
                   aaTCCGGGGAATTTCCTGGCGCCCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGaa                                                                 <  2:578330/90‑1 (MQ=255)
                   aaTACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGaa                                                                 <  2:59622/90‑1 (MQ=255)
                    aTACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAg                                                                >  2:389425/1‑90 (MQ=255)
                       cGGCGATTTTCCTGCCGCCCCAGGGCATAACCTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTa                                                             <  2:425154/90‑1 (MQ=255)
                       cGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTa                                                             >  1:277163/1‑90 (MQ=255)
                       cGGCAAATTCCCTGACCCCCCGGGGCATCACCTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTa                                                             <  2:580464/90‑1 (MQ=255)
                              tttCCTGACGCTCCAGGGATTCACCTGCAGCGTCCACTGCCGGGTAGGTTTTCCCCGTCCCCGCCGGGCCAACGCCGAAGGTAATGTCAt                                                      <  1:120815/90‑1 (MQ=255)
                                         tCCAGGGCATCACCTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAata                                           <  1:261182/90‑1 (MQ=255)
                                                     aCTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTAc                               <  1:154130/90‑1 (MQ=255)
                                                                          taggtTTTCCCCGTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTgcgc          <  2:379531/90‑1 (MQ=255)
                                                                          taggtTTTCCCCGTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTgcgc          <  1:139985/90‑1 (MQ=255)
                                                                             gtTTTCCCCGTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGc       <  2:83315/90‑1 (MQ=255)
                                                                              ttttACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCt      >  1:426175/1‑90 (MQ=255)
                                                                                  acccgtacccgCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAAt  >  1:616174/1‑90 (MQ=255)
                                                                                  |                                                                                         
CTGCCGGACGAGTCAGCAGAATACGGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAAT  >  NZ_CP009273/688312‑688483

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