Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,342,8900GC53.8% 6.3 / 13.3 13G93G (GGC→GGGmtgAmonofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (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 = 1.13e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACC  >  NZ_CP009273/3342801‑3342971
                                                                                         |                                                                                 
gtctgtTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAg                                                                                   >  2:193411/1‑90 (MQ=255)
              cGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAgg                                                                     >  1:151065/1‑90 (MQ=255)
                                             cgcTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGcc                                      <  1:193411/90‑1 (MQ=255)
                                             cgcTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGcc                                      >  2:281182/1‑90 (MQ=255)
                                               ctctTTGTGCGCCGGGGTTTTCTCAAGGGAAGCGACATCAAACCCCCAGTGCTCAGGAAATTTCTGATCTTCGGCGGCAATCACCGCCAg                                    <  2:95548/90‑1 (MQ=255)
                                               cTCGTTGTGCCCCGGGGTTTTCTAAAGGGAAGCGAAATCAAACCCCCAGTGCTCAGGAAATTTCTGATTTTCGGCGGCAATCACCGCCAg                                    <  2:63573/90‑1 (MQ=255)
                                                           cAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACgg                        >  2:224427/1‑90 (MQ=255)
                                                             ggtcTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCg                      <  1:313670/87‑1 (MQ=255)
                                                             gggCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCg                      >  1:66633/1‑90 (MQ=255)
                                                                  ttCTCAAGGGAAGCGACATCAAACCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGAtt                 <  1:346889/90‑1 (MQ=255)
                                                                  ttCTCAAGGGAAGCGACATCAAACCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGAtt                 <  2:70976/90‑1 (MQ=255)
                                                                       aaTGGAAGCGACATCAAACCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATc            <  1:204208/90‑1 (MQ=255)
                                                                                gACATCAAACCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAACCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGAc   <  1:310336/90‑1 (MQ=255)
                                                                                 aCATAAAACCCCCAGTGCTCAGGAATTTTCTGTTCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGAcc  <  1:367306/90‑1 (MQ=255)
                                                                                         |                                                                                 
GTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCTTTCTCAATGGAAGCGACATCAAAGCCCCAGTGCTCAGGAAATTTCTGATCTTCTGCGGCAATCACCGCCAGTCCCATCCACGGCGAGATTTGATCCATACTGACC  >  NZ_CP009273/3342801‑3342971

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