Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,329,9450TG58.8% 10.7 / 17.1 17T132P (ACC→CCC) adiAarginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (10/0);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 3.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGC  >  NZ_CP009273/4329861‑4330031
                                                                                    |                                                                                      
gcgcTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGGAGCGGGGCATcg                                                                                   >  1:253097/1‑90 (MQ=255)
               atTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGGCATCGCGGCAACGGcgcggc                                                                    >  1:65619/1‑88 (MQ=255)
               atTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGGCATCGCGGCAACGGCGCGTc                                                                    >  2:126930/1‑90 (MQ=255)
               atTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGGAGCGGGGCATCCCGGCCACGgcgccgc                                                                    >  1:359466/1‑87 (MQ=255)
                  cATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCgg                                                                 <  2:55282/90‑1 (MQ=255)
                    tGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGGCATCGCGGCAACGGCGCGGCCGGCg                                                               >  1:35573/1‑90 (MQ=255)
                                   ttCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTcggcggga                                                >  1:107729/1‑88 (MQ=255)
                                    tcatcaGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGGGGGGCATCGCGGCAAAGGGGCGGCCGGCGATAAAAGcggcggggt                                               >  1:3577/1‑87 (MQ=255)
                                      atcaGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCt                                             <  1:200443/90‑1 (MQ=255)
                                         aGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGGCATCGCGGCAACGGGGCGTCCCGCGCTAAAATCGGCGGGAACTTcc                                          >  2:254376/1‑90 (MQ=255)
                                         aGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGGAGCGGGGCATCGCGGCAACGGGGCGTCCCGCGATAAAGTCGGCGGGATCTTcc                                          >  1:352257/1‑90 (MQ=255)
                                                      gTGGCGGCAACAGCTGCTGGCGGGAGCGGGGCATCCCGGCAACGGGGCGGCCGGGGAAAAAGGCGGGGGGATCTTCCCGAAACCCggggg                             >  1:345513/1‑87 (MQ=255)
                                                           ggcAACAGCTGCTGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCa                        >  2:99509/1‑90 (MQ=255)
                                                           ggcAACAGCTGCTGGCGGTAGCGGGGCATCGCGGCAACGGCGCGTCCGGCGATAAAGGCGGCGGGATCTTCCCGAATCCAGGCGAATTCa                        >  1:132014/1‑90 (MQ=255)
                                                                  gctgctGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACaa                 <  2:107729/90‑1 (MQ=255)
                                                                  gctgctGGCGGTAGCGAGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACaa                 <  1:296980/90‑1 (MQ=255)
                                                                                 gggTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTcgc  <  1:232816/90‑1 (MQ=255)
                                                                                    |                                                                                      
GCGCTGCCCAGGAATATTCATGGATGTCACTATATTTCATCAGCGCGCTGAACAGTGGCGGCAACAGCTGCTGGCGGTAGCGGGTCATCGCGGCAACGGCGCGTCCGGCGATAAAGTCGGCGGTATCTTCCAGAATCCAGGCGAATTCATCGACAAGCTCCAGCAGGTCGC  >  NZ_CP009273/4329861‑4330031

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