Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,432,4730AG47.2% 26.4 / 14.6 36S857P (TCT→CCT) nifJpyruvate:ferredoxin (flavodoxin) oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (16/3);  new base G (0/17);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.35e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCAAGGCCAAATTCGGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAA  >  NZ_CP009273/1432388‑1432559
                                                                                     |                                                                                      
ccAAGGCCAAATTCGGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGc                                                                                    <  2:1250201/90‑1 (MQ=255)
        aatttCGGCATTATCTTAAAAAAGAGAGTTCGCCCAGGCCGCCCCCCGCCCGTTGGCATGGGTGGTATACGGTGTAGGGGGCAGGTTAcc                                                                            <  1:768603/87‑1 (MQ=255)
              ggCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATa                                                                      >  2:424777/1‑90 (MQ=255)
                    aTCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTga                                                                >  2:170932/1‑90 (MQ=255)
                         cAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGGGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAAc                                                           >  2:860819/1‑90 (MQ=255)
                          aaaTAGAGATTTCCCCCATGCCGGCCCACGACCGTTGGCATGGGTGGTATCCGGTGTAGGGGGCAGGTTACCGCCATAAATTGAAGAACa                                                          <  2:1260000/90‑1 (MQ=255)
                               gagagTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAg                                                     >  2:523662/1‑90 (MQ=255)
                                agagTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGt                                                    >  2:1111325/1‑90 (MQ=255)
                                   gTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGc                                                 >  2:972075/1‑90 (MQ=255)
                                   gTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGc                                                 >  1:1259999/1‑90 (MQ=255)
                                         cggggccgccCCCCGCCCTTGGGCTTGGGTGGTATAGGGTGTGGGGGGCGGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGc                                           <  1:113306/86‑1 (MQ=255)
                                         cggggccgccCCACGCCCGTGGGAAGGGGTGTTATCGGGTGAAGGGGGGAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGc                                           <  1:1470789/86‑1 (MQ=255)
                                            tGCCGCCCCCCGCCCTTTGGCATGGGGGGTATACGGTGAAGGGGGAGGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGAt                                        <  2:578497/90‑1 (MQ=255)
                                            tGCCGCCCCACGCCCGTGGGCATGGGTGGTATAGGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGAt                                        <  2:709904/90‑1 (MQ=255)
                                             gccgccCCACGCCGGTGGGCATGGGGGTTATGGGGTGGGGGGGGCGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATc                                       <  1:795252/90‑1 (MQ=255)
                                                   ccACGACCGTGGGCATGGGTGGTATAGGGTGTAGGGGGCGGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATc                                 <  2:1643652/90‑1 (MQ=255)
                                                      cGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATcctg                              >  2:1470790/1‑88 (MQ=255)
                                                        ccggTTGGAATGGGTGGTATGGGGTGTAGGGGGCAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTc                            <  1:645829/87‑1 (MQ=255)
                                                        cccGTGGGCATGGGGGTTATCGGGTGAGGGGGGCAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTc                            <  1:1111324/89‑1 (MQ=255)
                                                        aCCGTGGGCAGGGGTGATATGGGGTGAAGGGGGGAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTc                            <  1:424777/90‑1 (MQ=255)
                                                           gTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGcc                         >  2:1047998/1‑90 (MQ=255)
                                                             tGGCATCGGTGGTATCCGGTGTAGAGGGAAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCAt                       <  1:366678/90‑1 (MQ=255)
                                                             tGGCACCGGTGTTACACGTTTTAGGGGGCAGGTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCAt                       <  2:119252/90‑1 (MQ=255)
                                                                   gggtggtATCGGGTGTAGGGGGCGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCt                 <  2:1428177/89‑1 (MQ=255)
                                                                    ggtggtATCGGGTGAGGGGGGAAGTTCCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTg                <  1:357209/90‑1 (MQ=255)
                                                                     gtggtATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGa               >  1:1205700/1‑90 (MQ=255)
                                                                        gtATGCGGTGAAGGGGGAGGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTc            <  2:1272539/90‑1 (MQ=255)
                                                                           tggggTGTAGGGGGGGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGt         <  2:827064/87‑1 (MQ=255)
                                                                               gtgtAGGGGGCGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt     <  2:1205701/90‑1 (MQ=255)
                                                                               gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACGGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt     >  1:1846685/1‑90 (MQ=255)
                                                                               gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCTTAGAGCTGAGTCAGTAAtt     >  1:210951/1‑90 (MQ=255)
                                                                               gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTaact     >  1:1428176/1‑88 (MQ=255)
                                                                               gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt     >  2:1413703/1‑90 (MQ=255)
                                                                               gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt     >  2:1449913/1‑90 (MQ=255)
                                                                                ggtAGGGGGCAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAttt    <  1:523662/89‑1 (MQ=255)
                                                                                 gtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTa   >  1:619671/1‑90 (MQ=255)
                                                                                  gagagGGCAGGTTACCGCCATAAATTAACGAACAGCAAGTGGGGGTAGAGATGAACATCCGGCCGCCATAGAGATGAGCCAGTAATTTaa  >  1:170060/2‑90 (MQ=255)
                                                                                     |                                                                                      
CCAAGGCCAAATTCGGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAA  >  NZ_CP009273/1432388‑1432559

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