Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,434,377 A→G F222S (TTC→TCC)  nifJ ← pyruvate:ferredoxin (flavodoxin) oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,434,3770AG96.2% 77.7 / ‑3.7 26F222S (TTC→TCC) nifJpyruvate:ferredoxin (flavodoxin) oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (17/8);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CCTGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCA  >  NZ_CP009273/1434290‑1434454
                                                                                       |                                                                             
ccTGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGgaa                                                                             >  1:554424/1‑90 (MQ=255)
   gTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATaa                                                                          <  2:508342/90‑1 (MQ=255)
   gTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATaa                                                                          >  1:55501/1‑90 (MQ=255)
   gTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATaa                                                                          >  2:36107/1‑90 (MQ=255)
    tCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAg                                                                         <  1:228122/90‑1 (MQ=255)
             agaAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGAt                                                                >  1:437812/1‑90 (MQ=255)
                 aTCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGc                                                            <  2:437812/90‑1 (MQ=255)
                               gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTc                                                                     <  1:423366/67‑1 (MQ=255)
                               gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTc                                                                     >  2:423366/1‑67 (MQ=255)
                               gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                              >  1:368694/1‑90 (MQ=255)
                                       tGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                      <  2:313923/90‑1 (MQ=255)
                                                 ccGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGt                            >  2:400837/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                      <  1:426757/78‑1 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                      >  2:426757/1‑78 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:508623/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:454921/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:545242/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:206833/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  1:668420/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:638451/1‑90 (MQ=255)
                                                   gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:671042/1‑90 (MQ=255)
                                                             aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCgg                              <  2:285304/76‑1 (MQ=255)
                                                             aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCgg                              >  1:285304/1‑76 (MQ=255)
                                                             aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCggg                >  2:526701/1‑90 (MQ=255)
                                                                  ttGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGa           >  2:589225/1‑90 (MQ=255)
                                                                     gtggCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGa        <  2:632256/90‑1 (MQ=255)
                                                                           tCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGGTCCGGGTTGAGTGCCCGGGGGCGATGAGCATCa  >  1:327203/1‑90 (MQ=255)
                                                                                       |                                                                             
CCTGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCA  >  NZ_CP009273/1434290‑1434454

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