Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,004,302 G→A 51.3% S483F (TCT→TTT)  hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,004,3020GA51.3% ‑6.1 / 121.7 78S483F (TCT→TTT) hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (15/25);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

AGAGTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTC  >  minE/2004232‑2004367
                                                                      |                                                                 
agagTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgag                                                                   <  1:657532/71‑1 (MQ=255)
agagTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgag                                                                   <  1:2592813/71‑1 (MQ=255)
agagTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgag                                                                   <  1:3098401/71‑1 (MQ=255)
  agTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgagag                                                                 <  1:1333030/71‑1 (MQ=255)
   gTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTgaaaga                                                                <  1:1795757/71‑1 (MQ=255)
    ttCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGaaagag                                                               <  1:3398011/71‑1 (MQ=255)
    ttCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGaaagag                                                               <  1:3558389/71‑1 (MQ=255)
      ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCa                                                             >  1:896900/1‑71 (MQ=255)
      ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCa                                                             >  1:112286/1‑71 (MQ=255)
      ccAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCa                                                             >  1:336003/1‑71 (MQ=255)
       cAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCAt                                                            <  1:3002693/71‑1 (MQ=255)
       cAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCAt                                                            <  1:1292300/71‑1 (MQ=255)
       cAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCAt                                                            <  1:3475969/71‑1 (MQ=255)
          gTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCa                                                             <  1:1664546/67‑1 (MQ=255)
            tGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACcg                                                        >  1:2894769/1‑70 (MQ=255)
               tGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACcgcgcg                                                    >  1:3975785/1‑71 (MQ=255)
               tGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACcgcgc                                                     >  1:3367490/1‑70 (MQ=255)
                 gAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCgcggc                                                  >  1:1966547/1‑71 (MQ=255)
                   acaacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGgcgc                                                >  1:666882/1‑71 (MQ=255)
                   acaacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGgcgc                                                >  1:2757327/1‑71 (MQ=255)
                   acaacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGgcgc                                                >  1:233844/1‑71 (MQ=255)
                   acaacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGgcgc                                                >  1:149713/1‑71 (MQ=255)
                    caacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCt                                               >  1:3366590/1‑71 (MQ=255)
                    caacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCt                                               >  1:1355308/1‑71 (MQ=255)
                    caacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCt                                               <  1:3552034/71‑1 (MQ=255)
                    caacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCt                                               <  1:2415114/71‑1 (MQ=255)
                    caacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCt                                               <  1:3146843/71‑1 (MQ=255)
                     aacCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATAC‑‑CGCGGCGCTTCg                                            <  1:2802897/71‑1 (MQ=255)
                      acCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCACGCGGCGCTTc                                             <  1:1358542/71‑1 (MQ=255)
                      acCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTc                                             <  1:2064683/71‑1 (MQ=255)
                      acCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTc                                             <  1:1902575/71‑1 (MQ=255)
                       ccGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCg                                            <  1:2853343/71‑1 (MQ=255)
                       ccGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCg                                            <  1:1892518/71‑1 (MQ=255)
                       ccGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTc                                             <  1:3500845/70‑1 (MQ=255)
                       ccGCCTGGTAGTTGCTGATGATACCGTCTTTAATAAACATCCAGTGAGAGAGCATACCGCGCGGCGCTTCg                                            <  1:59807/71‑1 (MQ=255)
                        cGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTc                                             >  1:1807128/1‑69 (MQ=255)
                               tAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:3964480/71‑1 (MQ=255)
                               tAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:3756188/71‑1 (MQ=255)
                               tAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:1673158/71‑1 (MQ=255)
                               tAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:4252717/71‑1 (MQ=255)
                                  ttGCTGATGATACCGTCTTTAATAACCATCCAGTGaaagag                                                               <  1:669759/41‑1 (MQ=255)
                                  ttGCTGATGATACCGTCTTTAATAACCATCCAGTGaaagag                                                               <  1:3058984/41‑1 (MQ=255)
                                    gCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAaca                                >  1:961349/1‑70 (MQ=255)
                                    gCTGATGATACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAaca                                <  1:3830207/70‑1 (MQ=255)
                                        atgatACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCttt                           >  1:3156050/1‑71 (MQ=255)
                                        atgatACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCttt                           >  1:442972/1‑71 (MQ=255)
                                            tACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAAt                       >  1:3510646/1‑71 (MQ=255)
                                            tACCGTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAAt                       >  1:2507850/1‑71 (MQ=255)
                                              ccGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAAtt                      >  1:3396703/1‑70 (MQ=255)
                                                gTCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCAc                   <  1:3926249/71‑1 (MQ=255)
                                                 tCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:936907/53‑1 (MQ=255)
                                                 tCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:3955265/53‑1 (MQ=255)
                                                 tCTTTAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:3115646/53‑1 (MQ=255)
                                                  cTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCg                                            <  1:3479314/44‑1 (MQ=255)
                                                   tttAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCaa                                  >  1:415231/1‑53 (MQ=255)
                                                   tttAATAACCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCaa                                  >  1:2410470/1‑53 (MQ=255)
                                                      aataaCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGtt              >  1:572628/1‑70 (MQ=255)
                                                       ataaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCa                    >  1:2131068/1‑63 (MQ=255)
                                                       ataaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCa                    >  1:3600432/1‑63 (MQ=255)
                                                        taaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTc                                             <  1:2971467/37‑1 (MQ=255)
                                                        taaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTc                                             <  1:3860320/37‑1 (MQ=255)
                                                        taaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTc                                             <  1:252302/37‑1 (MQ=255)
                                                        taaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTc                                             <  1:2650152/37‑1 (MQ=255)
                                                        taaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGc            >  1:3840970/1‑70 (MQ=255)
                                                        taaCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCAACACCttt                           >  1:2407506/1‑54 (MQ=255)
                                                        taaCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTc                                             <  1:3061568/37‑1 (MQ=255)
                                                        taaCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTc                                             <  1:920400/37‑1 (MQ=255)
                                                         aaCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGaa                                       >  1:1252218/1‑42 (MQ=255)
                                                         aaCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGaa                                       >  1:1757765/1‑42 (MQ=255)
                                                         aaCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGaa                                       >  1:840531/1‑42 (MQ=255)
                                                          aCCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGa                                           >  1:153332/1‑37 (MQ=255)
                                                          aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCgg         <  1:3944594/71‑1 (MQ=255)
                                                          aCCATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCg          >  1:175698/1‑70 (MQ=255)
                                                           ccATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGa        <  1:2606955/71‑1 (MQ=255)
                                                           ccATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGa        <  1:95639/71‑1 (MQ=255)
                                                           ccATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGa        <  1:1418984/71‑1 (MQ=255)
                                                           ccATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGa        <  1:1752829/71‑1 (MQ=255)
                                                           ccATCCAGTGAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGa        <  1:3821176/71‑1 (MQ=255)
                                                           ccATCCAGTAAAAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGa        <  1:1716682/71‑1 (MQ=255)
                                                             aTCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAAt      <  1:2345342/71‑1 (MQ=255)
                                                              tCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGcc                                    <  1:2593495/40‑1 (MQ=255)
                                                               ccAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGt    >  1:4124662/1‑71 (MQ=255)
                                                               ccAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGt    >  1:1918869/1‑71 (MQ=255)
                                                                 aGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTc  <  1:523973/71‑1 (MQ=255)
                                                                 aGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTc  <  1:1588384/71‑1 (MQ=255)
                                                                 aGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTc  <  1:2182353/71‑1 (MQ=255)
                                                                 aGTGAGAGAGCATACAGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTc  <  1:171824/71‑1 (MQ=255)
                                                                      |                                                                 
AGAGTTCCAGGTTGATGGAACAACCGCCTGGTAGTTGCTGATGATACCGTCTTTAATAACCATCCAGTGAGAGAGCATACCGCGCGGCGCTTCGAGGAAGCCAACACCTTTGAATTCACCCGTTGCCGGAATGTTC  >  minE/2004232‑2004367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: