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

TGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCAATTTCGAC  >  NZ_CP009273/1434292‑1434462
                                                                                     |                                                                                     
tGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATa                                                                                   >  2:795757/1‑90 (MQ=255)
tGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATa                                                                                   >  1:141859/1‑90 (MQ=255)
      ggCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTc                                                                             <  2:574392/90‑1 (MQ=255)
      ggCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTc                                                                             <  1:1074563/90‑1 (MQ=255)
                    tcatCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACg                                                               >  2:1299814/1‑90 (MQ=255)
                      atCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTa                                                             >  2:1480008/1‑90 (MQ=255)
                         gCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACcg                                                          >  1:214871/1‑90 (MQ=255)
                            tGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCgg                                                       <  2:974827/90‑1 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAg                                                                                 >  1:1289904/1‑63 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAg                                                                                 <  2:1289904/63‑1 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                                      >  2:1058364/1‑90 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                                      >  2:335326/1‑90 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                                      >  2:159446/1‑90 (MQ=255)
                             gTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGa                                                      >  1:249109/1‑90 (MQ=255)
                                     tGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGa                                              <  2:583461/90‑1 (MQ=255)
                                       gTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAAt                                                                                    >  2:1161487/1‑50 (MQ=255)
                                       gTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAAt                                                                                    <  1:1161487/50‑1 (MQ=255)
                                       gTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATg                                            >  2:281940/1‑90 (MQ=255)
                                       gTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATg                                            >  2:1051431/1‑90 (MQ=255)
                                               ccGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGt                                    <  1:349700/90‑1 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTAGCGGGGATCACCGGGTGTTCCGGGTTg                                  >  1:1473726/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCCCCGGATGTTCCGGGGTg                                  >  1:931252/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGGTGGTCCGGGTTg                                  >  1:1412309/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                                  >  1:170714/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                                  >  1:1411241/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                                  >  2:372240/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                                  >  1:767317/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                                  >  1:407944/1‑90 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTg                                  <  1:297320/90‑1 (MQ=255)
                                                 gcgTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGGTCCGGGTTg                                  >  1:1316310/1‑90 (MQ=255)
                                                        aCCATGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGccc                           <  2:1411241/90‑1 (MQ=255)
                                                           aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCggg                        >  2:872787/1‑90 (MQ=255)
                                                           aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCggg                        >  1:154421/1‑90 (MQ=255)
                                                           aTGGGTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCggg                        >  2:1486087/1‑90 (MQ=255)
                                                             gggTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGgcg                      <  1:1397349/90‑1 (MQ=255)
                                                              ggTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGgcgc                     >  2:146434/1‑90 (MQ=255)
                                                              ggTTGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGgcgc                     <  1:335326/90‑1 (MQ=255)
                                                                ttGGTGGCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGa                   >  2:317309/1‑90 (MQ=255)
                                                                   gtggCTTCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGa                <  1:1120623/90‑1 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc           <  2:214871/90‑1 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc           >  2:1118673/1‑90 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc           <  1:1233058/90‑1 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc           >  2:78216/1‑90 (MQ=255)
                                                                        ttCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATc           <  1:566714/90‑1 (MQ=255)
                                                                         tCGCGAGACTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCa          >  2:465539/1‑90 (MQ=255)
                                                                         tCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCgg                                                       <  1:1397562/45‑1 (MQ=255)
                                                                         tCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCgg                                                       >  2:1397562/1‑45 (MQ=255)
                                                                               gaCTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCAATTTCg    <  2:636947/90‑1 (MQ=255)
                                                                                 cTGGGAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCAATTTCGAc  >  2:1108555/1‑90 (MQ=255)
                                                                                     |                                                                                     
TGTCGCGGCAGAGAAATCATTCATCGCCTGTTCAACATGGTCATAGACCGCGTTGTACCATGGGTTGGTGGCTTCGCGAGACTGGAAATAAGTGTCAGGATTGGCGGACGTACCGCGGATCACCGGATGTTCCGGGTTGAGTGCCCGGGCGCGATGAGCATCAATTTCGAC  >  NZ_CP009273/1434292‑1434462

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