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,3770AG93.1% 79.0 / ‑0.4 29F222S (TTC→TCC) nifJpyruvate:ferredoxin (flavodoxin) oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (16/11);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

TGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCA  >  NZ_CP009273/1434292‑1434454
                                                                                     |                                                                             
tGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATa                                                                           <  2:685350/90‑1 (MQ=255)
tGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATa                                                                           <  1:696589/90‑1 (MQ=255)
               aTCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGc                                                            <  1:281062/90‑1 (MQ=255)
                    tcatCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACg                                                       <  2:33877/90‑1 (MQ=255)
                          ccTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACcgc                                                 <  1:741189/90‑1 (MQ=255)
                          ccTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACcgc                                                 <  2:279225/90‑1 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                              >  1:124244/1‑90 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                              <  2:531401/90‑1 (MQ=255)
                                cAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCa                                           <  1:645317/90‑1 (MQ=255)
                                  aCATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCAcc                                         >  2:365256/1‑90 (MQ=255)
                                     tGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                      <  1:641088/90‑1 (MQ=255)
                                                cgcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGtt                           >  1:741630/1‑90 (MQ=255)
                                                cgcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGtt                           >  2:101345/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                      <  1:731037/78‑1 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                      >  2:731037/1‑78 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:441303/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:381197/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  2:531591/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  1:460662/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                          >  1:41059/1‑90 (MQ=255)
                                                      gTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGc                     >  1:280484/1‑90 (MQ=255)
                                                           aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCggg                >  1:774212/1‑90 (MQ=255)
                                                           aTGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGAATCACCGGATGTTCCGGGTTGAGTGCCCggg                >  1:175243/1‑90 (MQ=255)
                                                                ttGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGa           >  1:752244/1‑90 (MQ=255)
                                                                ttGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGa           >  1:713502/1‑90 (MQ=255)
                                                                ttGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGa           >  1:396009/1‑90 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc   <  2:713502/90‑1 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc   <  2:819766/90‑1 (MQ=255)
                                                                         tCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCa  >  1:317418/1‑90 (MQ=255)
                                                                                     |                                                                             
TGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCA  >  NZ_CP009273/1434292‑1434454

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