Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,849,435 T→G I139L (ATC→CTC)  hycB ← hydrogenase 3, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,849,4350TG98.2% 194.1 / ‑2.8 55I139L (ATC→CTC) hycBhydrogenase 3, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (31/23);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AAGTCCGTGTTAAAGGTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGA  >  NC_000913/2849287‑2849581
                                                                                                                                                    |                                                                                                                                                  
aaGTCCGTGTTAAAGGTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGcg                                                                                                                                                  <  2:680017/151‑1 (MQ=255)
       tgtTAAAGGTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGaa                                                                                                                                            >  2:287960/1‑150 (MQ=255)
             aGGTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTgg                                                                                                                                     >  2:582130/1‑151 (MQ=255)
              ggTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGc                                                                                                                                    >  1:560877/1‑151 (MQ=255)
               gTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACTGCGAGCGCGCGAATACCTGGc                                                                                                                                    >  2:335348/1‑150 (MQ=255)
                  aGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCAccc                                                                                                                                >  1:342448/1‑151 (MQ=255)
                  aGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCAccc                                                                                                                                >  2:281719/1‑151 (MQ=255)
                   gCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCAccc                                                                                                                                >  2:418404/1‑150 (MQ=255)
                      cACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTc                                                                                                                            <  2:719831/151‑1 (MQ=255)
                      cACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTc                                                                                                                            >  2:152210/1‑151 (MQ=255)
                      cACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTc                                                                                                                            >  2:138508/1‑151 (MQ=255)
                         gCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCa                                                                                                                           >  1:112747/1‑149 (MQ=255)
                                     tGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGac                                                                                                             >  1:237710/1‑151 (MQ=255)
                                        cgcgGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTTGCACCCAGTCAAGCAATGTGCTGTCAcg                                                                                                          <  1:736863/151‑1 (MQ=255)
                                        cgcgGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGacac                                                                                                           <  2:536333/150‑1 (MQ=255)
                                               gATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCgg                                                                                                    >  2:546428/1‑150 (MQ=255)
                                               gATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGa                                                                                                   <  2:46241/151‑1 (MQ=255)
                                               gATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGa                                                                                                   <  2:154305/151‑1 (MQ=255)
                                                atatCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCgg                                                                                                    >  1:69898/1‑149 (MQ=255)
                                                        gtTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGt                                                                                          <  2:549763/151‑1 (MQ=255)
                                                             ccaccaGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCgg                                                                                     >  2:248675/1‑151 (MQ=255)
                                                                   gATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGgcgc                                                                               <  1:620061/151‑1 (MQ=255)
                                                                    aTGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCt                                                                              <  1:693346/151‑1 (MQ=255)
                                                                        aGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCttt                                                                            >  1:78739/1‑149 (MQ=255)
                                                                         gggCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCgg                                                                         >  2:432734/1‑151 (MQ=255)
                                                                                   cGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGGATTGGCg                                                               >  2:29316/1‑151 (MQ=255)
                                                                                    gggCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGt                                                              <  1:503527/151‑1 (MQ=255)
                                                                                     ggCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGtt                                                             <  1:691681/151‑1 (MQ=255)
                                                                                        acaTCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTg                                                           <  1:609434/150‑1 (MQ=255)
                                                                                         caTCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGcc                                                         >  2:802886/1‑151 (MQ=255)
                                                                                          aTCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGGATTGGCGTTTGCCg                                                        >  2:371481/1‑151 (MQ=255)
                                                                                            ccGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCgg                                                       <  1:571182/150‑1 (MQ=255)
                                                                                                  gcaGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATcc                                                >  2:810398/1‑151 (MQ=255)
                                                                                                     ggCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGc                                             <  2:446018/151‑1 (MQ=255)
                                                                                                     ggCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGc                                             <  2:461635/151‑1 (MQ=255)
                                                                                                          gACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACg                                        >  2:546685/1‑151 (MQ=255)
                                                                                                            ccTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGc                                      <  2:801027/151‑1 (MQ=255)
                                                                                                                 ttCATCAAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGGATTGGCGTTTGCCGGAATATCCAGGGGGCGGCTGCCg                                 >  1:616970/1‑151 (MQ=255)
                                                                                                                  tcatcaAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCg                                 <  1:589970/150‑1 (MQ=255)
                                                                                                                    atcaAAGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGaa                              <  1:40833/151‑1 (MQ=255)
                                                                                                                         aGCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTc                         <  2:874742/151‑1 (MQ=255)
                                                                                                                          gCTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGaaa                             >  1:256540/1‑146 (MQ=255)
                                                                                                                           cTACAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCa                        >  1:892515/1‑150 (MQ=255)
                                                                                                                             aCAAAGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCaa                       <  2:560895/149‑1 (MQ=255)
                                                                                                                               aaaGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGAAGGCTGCCGGAAAATTCAAtt                     >  2:283302/1‑149 (MQ=255)
                                                                                                                               aaaGGTCACATTTGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTAATGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTg                    <  2:440766/150‑1 (MQ=255)
                                                                                                                                         tttGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACggg         >  2:294819/1‑151 (MQ=255)
                                                                                                                                          ttGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACggg         >  2:8685/1‑150 (MQ=255)
                                                                                                                                           tGACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCGTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGGAAATTCAATTGCGCCAAACCGGc        >  1:221106/1‑150 (MQ=255)
                                                                                                                                            gACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCa       <  1:587359/150‑1 (MQ=255)
                                                                                                                                            gACGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAg      <  1:248669/151‑1 (MQ=255)
                                                                                                                                              cGGCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGc    <  1:193664/151‑1 (MQ=255)
                                                                                                                                               ggCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCg   >  2:420569/1‑151 (MQ=255)
                                                                                                                                               ggCGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGGGGCGCTTTCGGGGTATTGGCGTTTGCCGGGATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAgg     >  2:410282/1‑149 (MQ=255)
                                                                                                                                                 cGAGCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGa  >  1:445210/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAGTCCGTGTTAAAGGTCAGCTCACGCTTACGTTTGCTGACGCGGGCGATATCGGTGTTATCCACCAGATGCAGGGCTTTAGTCGGGCACATCCGCACGCAGGCCGGACCTTGTTCATCAAAGCTACAAAGGTCACATTTGACGGCGATCGCGCGAATACCTGGCACCCAGTCAAGCAATGTGCTGACACGCGCCGGAGCAGGCGGTGCCGGTGGCGCTTTCGGGGTATTGGCGTTTGCCGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGA  >  NC_000913/2849287‑2849581

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