Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273169,1720CG53.3% 11.4 / 19.7 15G408G (GGC→GGGfhuBFe(3+)‑hydroxamate ABC transporter permease FhuB
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base G (8/0);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 1.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGT  >  NZ_CP009273/169083‑169252
                                                                                         |                                                                                 
aGCGGGGCGTTGCTCGAGGATTTAATGCGCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGc                                                                                   <  2:15396/90‑1 (MQ=255)
aGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGc                                                                                   <  1:107674/90‑1 (MQ=255)
aGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGc                                                                                   <  2:245476/90‑1 (MQ=255)
   ggggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGt                                                                                <  2:7347/90‑1 (MQ=255)
    gggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGta                                                                               >  2:45537/1‑90 (MQ=255)
             tCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGc                                                                      <  2:398050/90‑1 (MQ=255)
                   aTTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGTATTTTTTAGCGACTGa                                                                >  2:246531/1‑90 (MQ=255)
                         tGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGaa                                                         <  2:139111/90‑1 (MQ=255)
                          gCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGGGTCATGCTGGCGGGGGGGGGGTGTTTTTTTTAGCGGCTGACCGGgaa                                                        >  1:175598/1‑90 (MQ=255)
                              tGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGCGGGGTGTATTATTCAGCGACTGACCGGAAAACCg                                                    >  2:199167/1‑90 (MQ=255)
                                gcgcTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGTATTATTCAGCGGCCGACCGGAAAccccgt                                                  >  2:249042/1‑87 (MQ=255)
                                    tGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCa                                              <  1:87325/90‑1 (MQ=255)
                                    aggCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCa                                              <  2:50873/89‑1 (MQ=255)
                                         gcgAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGCGGGGTGTATTATTCAGCGACTGACCGGAAACCCCATGGGAAAccc                                         >  2:319486/1‑90 (MQ=255)
                                               tATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGt                                   <  1:83397/90‑1 (MQ=255)
                                                aTGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGTATTTTTTAGCGACTGACCGGGAAACCCCTGGGCAGCCCCGGAGGt                                   >  1:50873/1‑90 (MQ=255)
                                                      gcgcTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGGATTTTTTAGCGGCTGACCGGAAACCCGATGGCCAAGCCGGGAGGGGTgggg                            >  1:92580/1‑90 (MQ=255)
                                                                                ggtggcggGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCgt  <  2:319748/90‑1 (MQ=255)
                                                                                         |                                                                                 
AGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGT  >  NZ_CP009273/169083‑169252

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