Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,872,7300CG54.2% 2.5 / 37.3 24G253G (GGC→GGGdgcPdiguanylate cyclase DgcP
Reads supporting (aligned to +/- strand):  ref base C (7/4);  new base G (8/5);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGG  >  NZ_CP009273/1872648‑1872816
                                                                                  |                                                                                      
ataGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGtgatg                                                                                 <  2:424502/90‑1 (MQ=255)
ataGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGtgatg                                                                                 <  2:516821/90‑1 (MQ=255)
   gTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGt                                                                              >  1:819765/1‑90 (MQ=255)
   gTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGt                                                                              >  1:479055/1‑90 (MQ=255)
     ggggATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGttt                                                                            >  1:159440/4‑90 (MQ=255)
        gATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTa                                                                         >  1:310614/1‑90 (MQ=255)
        gATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTa                                                                         >  1:449467/1‑90 (MQ=255)
        gATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGGGATGAGGTTTTa                                                                         >  1:266223/1‑90 (MQ=255)
             gTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTAgttgt                                                                    <  2:819765/90‑1 (MQ=255)
                     ttCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTAgttgt                                                                    <  1:400089/82‑1 (MQ=255)
                     ttCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTAgttgt                                                                    >  2:400089/1‑82 (MQ=255)
                     ttCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTaa                                                            >  1:447035/1‑90 (MQ=255)
                      tCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTAGTTGTTTCACTaaa                                                           <  1:104471/90‑1 (MQ=255)
                         ggTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGGGATGAGGTTTTTGTTGTTTCACCAaacaa                                                        >  1:687506/1‑90 (MQ=255)
                                  gcgcCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGGGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGc                                               <  1:753828/90‑1 (MQ=255)
                                               cTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTtcgtcg                                  >  1:202872/1‑90 (MQ=255)
                                               cTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGGGGTGAGCTTTTAGTTGTTTTGCCTAATAAACAGAATTCCGAGATTtcgtcg                                  >  1:370647/1‑90 (MQ=255)
                                                   aGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTTGTTGTTTTACTAAAAAACGAGAATGCGGGTATTTCGTCGCGGc                              >  2:76441/1‑90 (MQ=255)
                                                           ggCGAAGTTATTGGTCGTCTCGGCGGTGGTGAGTTTTTAGTTGTTTCAGTAAACAACGAGAAAGCGGATATTTCGTCGCTGCGAGGACGc                      >  2:44727/1‑90 (MQ=255)
                                                             cGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCAt                    <  1:93073/90‑1 (MQ=255)
                                                                  ttATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTcagc               <  2:202872/90‑1 (MQ=255)
                                                                        gtcgtcTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCaaa         <  2:605926/90‑1 (MQ=255)
                                                                          cgtcTCGGGGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATa       >  1:3917/1‑90 (MQ=255)
                                                                           gtcTCGGGGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATAc      >  2:229646/1‑90 (MQ=255)
                                                                               cggcggTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTgg  <  2:447035/90‑1 (MQ=255)
                                                                                  |                                                                                      
ATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGG  >  NZ_CP009273/1872648‑1872816

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