Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,869,5390GC66.7% 34.5 / 25.1 27P99A (CCG→GCG) nimRDNA‑binding transcriptional regulator NimR
Reads supporting (aligned to +/- strand):  ref base G (3/6);  new base C (9/9);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACGCTGAGTCGTCGTCCTGTTTGCCAGGGTTAAAATTAATTCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTATGTTCCACCCCGCCCG  >  NZ_CP009273/1869454‑1869622
                                                                                     |                                                                                   
gCACGCTGAGTCGTCGTCCTGTTTGCCAGGGTTAAAATTAATTCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCAtt                                                                                 >  2:743629/1‑90 (MQ=255)
gCACGCTGAGTCGTCGTCCTGTTTGCCAGGGTTAAAATTAATTCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCAtt                                                                                 >  2:716161/1‑90 (MQ=255)
gCACGCTGAGCCGTCGTCCTGTTTGCCAGGGTTAAAATTAATTCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCAtt                                                                                 <  1:453197/90‑1 (MQ=255)
                                 aaaaTTAATTCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATaaacagaaaacag                                                <  1:7627/90‑1 (MQ=255)
                                    gttaattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGt                                             <  2:781540/89‑1 (MQ=255)
                                      taattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTc                                           <  1:743629/90‑1 (MQ=255)
                                      taattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTc                                           <  1:426029/90‑1 (MQ=255)
                                      taattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTc                                           >  2:35797/1‑90 (MQ=255)
                                       aattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCa                                          <  1:60737/90‑1 (MQ=255)
                                       aattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCa                                          <  2:742442/90‑1 (MQ=255)
                                       aattCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCa                                          >  2:389768/1‑90 (MQ=255)
                                           cccGGCTTAAGGGCGATGTCTTAAAAAAGCAAAAACTTGTGGGCAGTGTCACTGTTCTGGGTTCGATAAAAAGAAACCAGAGTCCAACAt                                      <  2:756415/90‑1 (MQ=255)
                                             cGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTa                                    >  1:290539/1‑90 (MQ=255)
                                             cGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTa                                    >  1:366866/1‑90 (MQ=255)
                                                cATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCt                                 <  2:817479/90‑1 (MQ=255)
                                                cATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCt                                 <  2:259810/90‑1 (MQ=255)
                                                         gATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGa                        >  1:44228/1‑90 (MQ=255)
                                                         gATATTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGa                        >  1:131806/1‑90 (MQ=255)
                                                         gATATTTTCAATGTGCAACATGTTGACGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTc                                           <  2:322771/71‑1 (MQ=255)
                                                         gATATTTTCAATGTGCAACATGTTGACGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTc                                           >  1:322771/1‑71 (MQ=255)
                                                          atatTTTCAATGTGCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGAt                       <  1:738763/90‑1 (MQ=255)
                                                             ttttCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGc                    <  1:229670/90‑1 (MQ=255)
                                                               ttCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTa                  <  2:559430/90‑1 (MQ=255)
                                                                    tgtgCAACATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTATGTTc             >  1:757203/1‑90 (MQ=255)
                                                                          aCATGTTGTCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTATGTTCCAcccc       <  1:770527/90‑1 (MQ=255)
                                                                               ttgtCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTATGTTCCACcccgcccg  >  2:405132/1‑90 (MQ=255)
                                                                               ttgtCGCCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTATGTTCCACcccgcccg  >  1:681761/1‑90 (MQ=255)
                                                                                     |                                                                                   
GCACGCTGAGTCGTCGTCCTGTTTGCCAGGGTTAAAATTAATTCCCGGCATAAGGGCGATATTTTCAATGTGCAACATGTTGTCGGCATTGTCACTGCGCTGGGTTCGATAAACAGAAAACAGAGTTCAGCATTAGCTGTTACTTGATTGCTATGTTCCACCCCGCCCG  >  NZ_CP009273/1869454‑1869622

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