Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,433,3420AG65.2% 22.0 / 21.9 23F567S (TTC→TCC) nifJpyruvate:ferredoxin (flavodoxin) oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (8/7);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAACTTCTTGCGGCAAGCGCGAC  >  NZ_CP009273/1433221‑1433476
                                                                                                                         |                                                                                                                                      
cagcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATg                                                                                                                       <  2:451316/139‑1 (MQ=255)
        gagCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTAt                                                                                                                >  1:447027/1‑138 (MQ=255)
                 aGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGgc                                                                                                      <  2:224757/139‑1 (MQ=255)
                                      ggCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGa                                                                                                                      >  2:339400/1‑102 (MQ=255)
                                      ggCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGa                                                                                                                      <  1:339400/102‑1 (MQ=255)
                                       gCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGc                                                                                 >  1:97503/1‑138 (MQ=255)
                                            tGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGAtttt                                                                           >  1:407822/1‑139 (MQ=255)
                                            tGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGAtttt                                                                           >  2:17071/1‑139 (MQ=255)
                                            tGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGAtttt                                                                           >  1:252538/1‑139 (MQ=255)
                                              cTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTcg                                                                         <  2:66239/139‑1 (MQ=255)
                                                  tGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTcgccgc                                                                     >  2:310/1‑139 (MQ=255)
                                                      aCTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGtt                                                                  <  1:310/138‑1 (MQ=255)
                                                      aCTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGtt                                                                  <  2:294887/138‑1 (MQ=255)
                                                       cTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTaa                                                                <  1:17071/139‑1 (MQ=255)
                                                           tGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGAc                                                                                                                     >  2:245246/1‑82 (MQ=255)
                                                           tGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGAc                                                                                                                     <  1:245246/82‑1 (MQ=255)
                                                                                   cGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGtt                                    >  2:447814/1‑139 (MQ=255)
                                                                                     aGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGa                                                                                                                      >  2:221612/1‑55 (MQ=255)
                                                                                     aGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGa                                                                                                                      <  1:221612/55‑1 (MQ=255)
                                                                                       ggCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAac                                <  1:341915/139‑1 (MQ=255)
                                                                                                      cAGAATTTGCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAActtct                 >  1:101246/1‑139 (MQ=255)
                                                                                                              gCGTCAGATGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAACTTCTTGCGGCa          <  2:447027/138‑1 (MQ=255)
                                                                                                                     aTGGGAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAACTTCTTGCGGCAAGCGCGAc  <  2:101246/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                      
CAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAACTTCTTGCGGCAAGCGCGAC  >  NZ_CP009273/1433221‑1433476

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