Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,074,9240AG52.2% 17.9 / 30.4 23D55G (GAT→GGT) putPsodium/proline symporter PutP
Reads supporting (aligned to +/- strand):  ref base A (0/11);  new base G (12/0);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 7.40e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTG  >  NZ_CP009273/1074837‑1075006
                                                                                       |                                                                                      
ggCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGata                                                                                      <  2:1614200/90‑1 (MQ=255)
                   tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAAt                                                                 <  1:1293231/90‑1 (MQ=255)
                                  gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGGGGGTATTAGCGGGTGGGTGGTTATTGGGGTGCCGGgggc                                                  >  2:1138575/1‑90 (MQ=255)
                                  gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTGTGGGGGGGGGGGTGGTGATGGGGGTGGCgggggg                                                  >  1:999745/1‑87 (MQ=255)
                                  gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGTGGGTGGTTATTGGGGTGGCgggggc                                                   >  2:1250857/1‑88 (MQ=255)
                                  gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGGGGGGGGTTATGGGGGTGGCgggggg                                                  >  2:773293/1‑87 (MQ=255)
                                         acaggTCTTGGGCCATTCGTGGCGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGttttt                                         <  2:1834899/86‑1 (MQ=255)
                                            aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGGGCGGGTGGGTGTTTATGGGGGTTGCGGGCGGCGtttttttt                                      >  2:646750/1‑90 (MQ=255)
                                                 tGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGGGGCTGGCTGTTTATGGGGGTGGCGGGGGCGGTtttttttttccg                                 >  2:186146/1‑86 (MQ=255)
                                                 tGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTAGCGGGGGGGTGGTTATTGGGGTGGCGGGGGGGGTtttttttttttg                                 >  2:909329/1‑86 (MQ=255)
                                                 tGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGGGGGGGGTTATTGGGGTGGCCGGGGGGGGGTtttttttcgt                                  >  2:1110626/1‑88 (MQ=255)
                                                 tGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATTAGGGGGTGGGTGTTTATGGGGGTGGCGGGGGCGGTTTTTTTTTTCgg                                 >  2:1424549/1‑90 (MQ=255)
                                                 tGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTGTGCTGGGCTGGCTGGTTATGGGGGTGTTCGGCGCTGTGtttttttttgg                                 >  2:58505/1‑86 (MQ=255)
                                                       aTTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATGAGGGGCGGGGTGTTTATGGGGGTGGCGGGGGCTGTTTTTTTTTTCCGGGtttt                            >  2:101079/1‑90 (MQ=255)
                                                       aTTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGGGGGCGGCTGTTTATGGGGGTGGCGGGGGCTGGTTTTTTTTTCTGGGGttt                            >  2:2054387/1‑90 (MQ=255)
                                                         tCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCg                         <  2:252115/90‑1 (MQ=255)
                                                          cGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGa                        <  2:571961/90‑1 (MQ=255)
                                                                         ggcggGTGCGTCGGATATGAGCGGCTGGCTGTTTATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAt         <  1:147557/90‑1 (MQ=255)
                                                                          gcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAtt        <  2:1798399/90‑1 (MQ=255)
                                                                              gTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGcc    <  1:1250857/90‑1 (MQ=255)
                                                                              gTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGcc    <  2:1863675/90‑1 (MQ=255)
                                                                              gTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGcc    <  2:803062/90‑1 (MQ=255)
                                                                                gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTg  <  2:2084094/90‑1 (MQ=255)
                                                                                gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTg  <  2:178524/90‑1 (MQ=255)
                                                                                       |                                                                                      
GGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTG  >  NZ_CP009273/1074837‑1075006

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: