Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,852,4630TG53.8% 2.7 / 10.4 13M671R (ATG→AGG) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (0/7);  total (3/10)
Fisher's exact test for biased strand distribution p-value = 6.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAA  >  NZ_CP009273/2852377‑2852546
                                                                                      |                                                                                   
acaAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggt                                                                                  >  2:54535/1‑90 (MQ=255)
        tCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGa                                                                          >  2:163788/1‑90 (MQ=255)
                          atcgCATCTTTCCCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAATCCGCCAATAttt                                                        <  2:72977/90‑1 (MQ=255)
                                    tccccGCGCGGGGGCGGCAGAGGCCCTGCCGGCGGGGCCCTAAACTTTTTTGGTGGAGTTTACTAAAACCGCCAATATTTTACATAACGc                                              <  2:187487/88‑1 (MQ=255)
                                    tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGc                                              >  1:204679/1‑90 (MQ=255)
                                                  cgGCAGATGCCCTGGCGTCCGGGCGCTCAACTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc                                <  1:335058/90‑1 (MQ=255)
                                                  cgGCAGATGCCCGGGCGTCGGGGCGCTCACCTTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc                                <  2:6748/90‑1 (MQ=255)
                                                   gccAGATGACCGGGCGCCCGGCCCCTCACCTTTTAGGGGGGAGATGACGAAACCCGCCAATATTTTACATAACGCCACCGAATACAGTCt                               <  2:281333/88‑1 (MQ=255)
                                                           aCCGGGGGTCGGGGCGCTCAACTTTTAGGGGGGAGATGTCAGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa                       <  2:204679/90‑1 (MQ=255)
                                                             cgggCGCCCGGGCCCCCAACCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATg                     <  1:396464/88‑1 (MQ=255)
                                                                    ccGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATgaga              >  1:376815/1‑90 (MQ=255)
                                                                      gggCGCCAACCCTTTAGGGGGGAGATGAGGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc            <  2:357412/90‑1 (MQ=255)
                                                                       gccgcTCAACCTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCg           <  2:376815/88‑1 (MQ=255)
                                                                                cttttAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa  <  2:25926/88‑1 (MQ=255)
                                                                                      |                                                                                   
ACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAA  >  NZ_CP009273/2852377‑2852546

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