Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,852,4770GA60.0% ‑3.9 / 10.0 15E676K (GAA→AAA) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (0/9);  total (4/11)
Fisher's exact test for biased strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT  >  NZ_CP009273/2852394‑2852555
                                                                                   |                                                                              
gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCg                                                                          >  2:155750/1‑90 (MQ=255)
    tATCGATCGCATTTTTCCCCGCGTAGGCGCGGCAGATGCCCTGGCGCCCGGGCGCTCAACCTTTATGGGGGAGATGACTAAAACCGCCaa                                                                      <  2:67360/90‑1 (MQ=255)
                        gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCg                                                  >  2:349208/1‑90 (MQ=255)
                             ggCGCGGCAGTTGCCCGGGCGCCCGGCCCCTAACCCTTTAGGGTGGAGAGGACGAAAACCGCCAATATTTTACATAACGCCACCGAATAc                                             <  2:435924/90‑1 (MQ=255)
                                 cgGCAGATGCCCTGGCGCCGGGGCCCCCACCCTTTATGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc                                         <  1:155750/90‑1 (MQ=255)
                                  ggCAGACGCCCGGGCATCGGGCCCCCCAACCTTTCGGGGGGGGTGGCCTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt                                        <  1:254537/90‑1 (MQ=255)
                                    cAGAGGCCCTGGCGCCGGGCCGCCCACCTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTgg                                      <  2:237375/90‑1 (MQ=255)
                                        ggcccTGGCGTCGGGGCGCTCAACTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGtt                                  <  2:204747/87‑1 (MQ=255)
                                                     gggCCCCCACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc                     <  1:359075/90‑1 (MQ=255)
                                                          cTCACCTTTTAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCg                <  1:45717/90‑1 (MQ=255)
                                                               cttttATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa           <  2:90401/88‑1 (MQ=255)
                                                                   tgtgggggAGTTGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc       <  2:182455/88‑1 (MQ=255)
                                                                        tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt  >  2:377617/1‑90 (MQ=255)
                                                                                 cTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc          <  1:38791/73‑1 (MQ=255)
                                                                                 cTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc          >  2:38791/1‑73 (MQ=255)
                                                                                   |                                                                              
GACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT  >  NZ_CP009273/2852394‑2852555

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: