Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,216,058 A→G V145A (GTC→GCC)  nfeF ← NADPH‑dependent ferric chelate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,216,0580AG84.3% 131.4 / 14.5 51V145A (GTC→GCC) nfeFNADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base G (21/22);  total (25/26)
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.00e+00

CCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATC  >  NC_000913/3215923‑3216192
                                                                                                                                       |                                                                                                                                      
ccATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATa                                                                                                                                    <  1:391575/140‑1 (MQ=255)
  aTCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACa                                                                                                                                  >  1:975057/1‑140 (MQ=255)
   tCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAg                                                                                                                                 >  1:884859/1‑140 (MQ=255)
       ggTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGa                                                                                                                             <  2:1018656/140‑1 (MQ=255)
          gcgcGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAc                                                                                                                          <  1:804696/140‑1 (MQ=255)
           cgcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATAcg                                                                                                                         >  2:755188/1‑140 (MQ=255)
           cgcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATAcg                                                                                                                         >  2:461665/1‑140 (MQ=255)
             cgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGa                                                                                                                             <  1:144123/134‑1 (MQ=255)
             cgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGa                                                                                                                             <  1:303364/134‑1 (MQ=255)
             cgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGa                                                                                                                             >  2:144123/1‑134 (MQ=255)
             cgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGa                                                                                                                             >  2:303364/1‑134 (MQ=255)
               agaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTa                                                                                                                     >  2:804477/1‑140 (MQ=255)
                 aTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAAt                                                                                                                   <  2:826200/140‑1 (MQ=255)
                   aaTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCt                                                                                                                 >  1:444299/1‑140 (MQ=255)
                    aTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCACAGGCATACAGCTGATACGCGTAATCtt                                                                                                                <  1:99548/140‑1 (MQ=255)
                          aCAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGca                                                                                                          >  2:586845/1‑140 (MQ=255)
                           cAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCgggcc                                                                                                         <  2:444299/140‑4 (MQ=255)
                           cAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGcac                                                                                                         <  2:701333/140‑1 (MQ=255)
                            aGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGcacc                                                                                                        >  1:221661/1‑140 (MQ=255)
                               cGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGcac                                                                                                         <  1:706678/136‑1 (MQ=255)
                               cGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGcac                                                                                                         >  2:706678/1‑136 (MQ=255)
                               cGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCaccacc                                                                                                     >  2:770831/1‑140 (MQ=255)
                                         cACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgc                                                                                           >  2:956135/1‑140 (MQ=255)
                                         cACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgc                                                                                           >  1:481145/1‑140 (MQ=255)
                                         cACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgc                                                                                           >  1:335695/1‑140 (MQ=255)
                                                aCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAcc                                                                                    >  1:1015275/1‑140 (MQ=255)
                                                  cAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTg                                                                                  <  1:11922/140‑1 (MQ=255)
                                                  cAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTg                                                                                  <  1:695247/140‑1 (MQ=255)
                                                  cAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTg                                                                                  <  2:486717/140‑1 (MQ=255)
                                                   aGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCt                                                                                   <  1:934195/138‑1 (MQ=255)
                                                   aGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCt                                                                                   >  2:934195/1‑138 (MQ=255)
                                                   aGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGc                                                                                 <  1:802573/140‑1 (MQ=255)
                                                    gcgcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGcc                                                                                <  2:858201/140‑1 (MQ=255)
                                                        aCTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCg                                                                            >  2:103374/1‑140 (MQ=255)
                                                        aCTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCg                                                                            >  2:624868/1‑140 (MQ=255)
                                                                      tAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGcc                                                                                <  1:803325/122‑1 (MQ=255)
                                                                      tAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGcc                                                                                >  2:803325/1‑122 (MQ=255)
                                                                         cGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGt                                                           <  1:62292/140‑1 (MQ=255)
                                                                                 ttGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTg                                                   <  2:884859/140‑1 (MQ=255)
                                                                                    cTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCAt                                                <  1:861726/140‑1 (MQ=255)
                                                                                         cGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGccc                                           <  2:335695/140‑1 (MQ=255)
                                                                                                aGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCt                                    <  1:586845/140‑1 (MQ=255)
                                                                                                     gcggcgCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGa                               <  1:103374/140‑1 (MQ=255)
                                                                                                     gcggcgCAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGa                               <  1:624868/140‑1 (MQ=255)
                                                                                                       ggcgCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGAcc                             >  2:570708/1‑140 (MQ=255)
                                                                                                          gcAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCa                          <  1:755188/140‑1 (MQ=255)
                                                                                                          gcAATGCAGGCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCa                          <  1:630510/140‑1 (MQ=255)
                                                                                                                  ggCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTg                  >  2:789184/1‑140 (MQ=255)
                                                                                                                  ggCATTCCGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTg                  <  1:770831/140‑1 (MQ=255)
                                                                                                                        ccGGATTCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATaaa            >  2:101687/1‑140 (MQ=255)
                                                                                                                          ggATTCAGCGCAGGCATACAGCTGATACGCGTAATCTTCGGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGa          <  2:341016/140‑1 (MQ=255)
                                                                                                                             ttCATCGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAaat       >  2:334396/1‑140 (MQ=255)
                                                                                                                                  cGCAGGCATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAaatcaatc  >  2:716152/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGTGAATAAAGAAATCAATC  >  NC_000913/3215923‑3216192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑