Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,224,1710AC54.5% 14.3 / 29.4 22P162P (CCA→CCCpgiglucose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/9);  new base C (12/0);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 2.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGC  >  NZ_CP009273/4224083‑4224254
                                                                                        |                                                                                   
cGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCAt                                                                                    <  2:435798/90‑1 (MQ=255)
 gATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATa                                                                                   <  1:344887/90‑1 (MQ=255)
          cGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGGCCTCCGCCCCCACATgggggg                                                                          >  2:279130/1‑86 (MQ=255)
               aGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTTTTACCCCCGCCCCCACCTGGTGGCCGAAg                                                                     >  1:29353/1‑90 (MQ=255)
                         ttATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTcc                                                           <  1:122228/90‑1 (MQ=255)
                          tataCCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTTACCTCCGCCCCCACATGGTGACCGAAGGTCTTCGTCCg                                                          >  2:311970/1‑90 (MQ=255)
                                 gCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTGACCCCCGCCCCCACAAGGTGACCCGAACATTGCGCCCGTGCaaaa                                                   >  1:389717/1‑90 (MQ=255)
                                 gCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCTACATTGGGAACGGAGCTCTGCGGCCGTACaaaa                                                   >  2:12906/1‑90 (MQ=255)
                                 gCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTCTGACCCCGGCCCCCACATTGGGGCCGAAGGGCTGCGTCCGCGCaaaa                                                   >  1:314955/1‑90 (MQ=255)
                                                aCGTAGTGAACATCGGGATCGGCGGTTCTTACCCCCGCCCCTACATGGGGGCCGAAGCTCTGCGTCCCTACAAAAAACACCCGAAAAtac                                    >  2:299427/1‑88 (MQ=255)
                                                aCGTAGTGAACATCGGGATCGGCGGTTCTGACCCCCGCCCCTACATGGGGGCCGAAGCTCTGCGTCCCTACAAAAAACACCCGAAAATGc                                    >  2:216664/1‑90 (MQ=255)
                                                    aGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACtt                                <  1:38524/90‑1 (MQ=255)
                                                     gTGAACATCGGGATCGGCGGTTCTGACCCCGGCCCCTACATGGGGACCGAAGCTCTGCGTCCGTACAAAAAACACCCGAACATGCACttt                               >  2:38524/1‑90 (MQ=255)
                                                         aCATCGGGATCGGCGGTTCTGACCTCGGCCCCTACATGGTGACCGAAGCTCTGCGTCCGTACAAAAAACACCCGAACATGCACTtttttt                           >  2:33960/1‑90 (MQ=255)
                                                         aCATCGGGATCGGCGGTTCTGACCCCGGCCCCTACATGGTGAACGAAGCTCTGCGTCCGGACAAAAAACACCCCAACAAGCACTtttttt                           >  1:163797/1‑90 (MQ=255)
                                                                    ggcggGTCTGACCCCCGCCCCAACATGGGGACCCAAGCTCTTCGTCCGTACAAAAAACACCTCAACATGCCATTTTTTTTTTACGTCCAt                >  2:571615/1‑90 (MQ=255)
                                                                     gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg               <  1:571615/90‑1 (MQ=255)
                                                                     gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg               <  2:34467/90‑1 (MQ=255)
                                                                     gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg               <  1:311970/90‑1 (MQ=255)
                                                                     gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg               <  2:575135/90‑1 (MQ=255)
                                                                                cTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATc    <  2:265416/90‑1 (MQ=255)
                                                                                  cGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATcgc  >  1:355360/1‑90 (MQ=255)
                                                                                        |                                                                                   
CGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGC  >  NZ_CP009273/4224083‑4224254

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