Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273765,6750AC53.8% 8.5 / 12.2 13T748P (ACC→CCC) mngBmannosylglycerate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base C (0/7);  total (4/9)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAA  >  NZ_CP009273/765589‑765750
                                                                                      |                                                                           
cTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCg                                                                          <  2:4795/90‑1 (MQ=255)
       gcCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGcc                                                                   >  2:67240/1‑90 (MQ=255)
                              tCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTcagc                                            <  1:67240/90‑1 (MQ=255)
                                ccGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTcagcag                                          <  1:145038/90‑1 (MQ=255)
                                               cGTGGTCTGTTTTTTTTCCGCCTAATTTTATTGATTTATCCCGGTACCCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTg                           <  1:306332/90‑1 (MQ=255)
                                               cGTGGTCTGTTTTCTTGTCGCCTAATTTTATTGATTTATCCCGGTCCGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTg                           <  1:66837/90‑1 (MQ=255)
                                               cGTGGTCTGTTTTCTTGTCGCCTAATTTTATTGAGTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTg                           <  1:21489/90‑1 (MQ=255)
                                                 tgttCTGTTTTCTTGTCGCCTAAGTTTATTGATTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGAc                         <  1:334825/87‑1 (MQ=255)
                                                   gTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTc                       >  1:476802/1‑90 (MQ=255)
                                                         tttCTTGTCGCTTAAGTTTATTGATTTATCCCGGCCCCCCACCCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTAc                 <  2:454343/90‑1 (MQ=255)
                                                             ttGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAgtg             >  1:547456/1‑90 (MQ=255)
                                                                 cGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTac         >  1:500009/1‑90 (MQ=255)
                                                                        gTTTATTGAGTTATCCCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACaaaa  <  1:143885/90‑1 (MQ=255)
                                                                                      |                                                                           
CTGGGCGGCCTTCGGGAATTAAAATGCCAGTCCCGGACTCACAACTACGTGGTCTGCTTTCTTGTCGCCTAAGTTTATTGAGTTATACCGGTACGCCAACCGCCGCTGGTGTAGCTCAGCAGGCGCGAGCATGGCTGACTCCAGTACAGTGTTACAACAAAA  >  NZ_CP009273/765589‑765750

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