Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,360,8100AG55.6% 14.1 / 23.5 18G87G (GGA→GGGarnCundecaprenyl‑phosphate 4‑deoxy‑4‑formamido‑L‑arabinose transferase
Reads supporting (aligned to +/- strand):  ref base A (0/8);  new base G (10/0);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 2.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGC  >  NZ_CP009273/2360723‑2360899
                                                                                       |                                                                                         
ggTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGtt                                                                                         >  2:1065917/1‑90 (MQ=255)
  tCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGttta                                                                                       >  1:779585/1‑88 (MQ=255)
  tCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCa                                                                                       <  2:1163805/90‑1 (MQ=255)
    gAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGtttatt                                                                                     >  2:1461189/1‑86 (MQ=255)
     aaGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTc                                                                                    <  1:730038/90‑1 (MQ=255)
                        cAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAat                                                                 <  2:504881/90‑1 (MQ=255)
                          gCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTTGTCCCGTTATTGGGGACTTAattt                                                               >  1:1464979/1‑89 (MQ=255)
                                  gtgtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGGTTTATCACGTTTCTTGCGGCTTTATTTTTTTCCtt                                                       >  2:516196/1‑90 (MQ=255)
                                    gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCACCGTTCTGGCGGCTTAATTTTTTccccttt                                                     >  1:1039881/1‑86 (MQ=255)
                                    gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCACCGTACTGGCGGGTTAATTTTTATCCtttt                                                     >  2:775837/1‑88 (MQ=255)
                                    gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGGTTAGGTACGTTACTTGCGGCTTAATTATTTACCTTGa                                                     >  1:969333/1‑90 (MQ=255)
                                          tttGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGa                                               <  1:1285118/90‑1 (MQ=255)
                                          tttGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGa                                               <  1:451358/90‑1 (MQ=255)
                                                     gCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCCCGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATc                                    <  1:1120012/90‑1 (MQ=255)
                                                             gggCAACATTCAGCGATTATGGCGGGGTTTAGTCCCCTTACTGGCGGCCTAATTATTACCCCTTATGCCGCTCTCCCGAATCCCCCcgaa                            >  1:1582811/1‑90 (MQ=255)
                                                             gggCAACATTCAGCGATTATGGCGGGGATCAGTCACGCTACTTGCGGCCTAATTTTTATCCTTTATGCCGGTCTCCACAAACCCCCacaa                            >  2:1527006/1‑90 (MQ=255)
                                                                              tATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCtggt           <  2:271870/90‑1 (MQ=255)
                                                                                  gCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCa       <  2:78282/90‑1 (MQ=255)
                                                                                       aTTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGc  <  2:739609/90‑1 (MQ=255)
                                                                                       |                                                                                         
GGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGC  >  NZ_CP009273/2360723‑2360899

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