Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,852,4660TG60.0% ‑1.6 / 13.5 20V672G (GTG→GGG) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (7/1);  new base G (0/12);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTC  >  NZ_CP009273/2852378‑2852550
                                                                                        |                                                                                    
caAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggtg                                                                                     <  1:211205/90‑1 (MQ=255)
                    tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa                                                                 >  2:436917/1‑90 (MQ=255)
                    tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa                                                                 >  2:114214/1‑90 (MQ=255)
                    tATCGACCGCATCTTTCCCCGCGTGGGCGCGGCAGATGCCCTGGCGCCCGGGCGCTCACCCTTTATGGGGGAGATGACTGAAACCGCCaa                                                                 <  1:81024/90‑1 (MQ=255)
                     atcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAat                                                                >  2:483085/1‑90 (MQ=255)
                                 tttCCCCGCGTGGGCGCGGCAGATGCCCGGGCGCCCGGGCGCTCAACTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAAc                                                    <  2:328279/90‑1 (MQ=255)
                                   tccccGCGCGGGCGCGCCAGGGCACCGGGCGCCCGGGCCCCCCACTTTTAGGGGGGAATTGACAAACACCGCCAATATTTTACATAACGc                                                  <  2:332516/88‑1 (MQ=255)
                                   cccccGCGTAGGCGCGGCAGATGCCCTGGCGTCCGGGCGCTCACCCTTTATGGGGGAGTTGACTAAAACCGCCAATATTTTACATAACGc                                                  <  1:543936/88‑1 (MQ=255)
                                      ccGCGTGGGCGCGGCAGATGACCTGGCGCCCGGGCCCTCAACCTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc                                               <  1:293988/90‑1 (MQ=255)
                                      ccGCGGGGGCGCGGCAGATGCCTGGGCGCCCGGGCGCTCAATTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc                                               <  2:322143/90‑1 (MQ=255)
                                      ccGCGAGGGCGCGGCAGATGCCCGGGCGCCCGGGCCCACACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc                                               <  2:280468/90‑1 (MQ=255)
                                       cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc                                              >  2:368784/1‑90 (MQ=255)
                                        gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCg                                             >  2:122214/1‑90 (MQ=255)
                                        gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCg                                             >  2:543936/1‑90 (MQ=255)
                                                  ggCAGATGCCCGGGCGCCGGGGCGCTCAACCTTTATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt                                   <  1:35/90‑1 (MQ=255)
                                                  ggCAGATGACCTGGCGTCCGGGCGCTCATCCTTTAtggt                                                                                      >  2:289235/1‑39 (MQ=255)
                                                  ggCAGATGACCTGGCGTCCGGGCGCTCATCCTTTAtggt                                                                                      <  1:289235/39‑1 (MQ=255)
                                                  ggCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt                                   >  1:29822/1‑90 (MQ=255)
                                                   gaaGATGCCCGGGCGCCGGGGCGCTCAACTTTTAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTg                                  <  1:122214/88‑1 (MQ=255)
                                                                     gggCCCTCACCTTTTATGGGGGAGATGACGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc                <  2:325967/90‑1 (MQ=255)
                                                                               cttttAGGGGGGAGATGATGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa      <  2:352247/88‑1 (MQ=255)
                                                                                  ttATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt   <  1:181425/90‑1 (MQ=255)
                                                                                   tATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc  <  1:436917/90‑1 (MQ=255)
                                                                                        |                                                                                    
CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTC  >  NZ_CP009273/2852378‑2852550

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