Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,004,952 G→A 26.3% A266A (GCC→GCT hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,004,9520GA26.3% 77.2 / 24.4 68A266A (GCC→GCThybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base G (14/36);  new base A (6/12);  total (20/48)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

TTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGA  >  minE/2004889‑2005017
                                                               |                                                                 
ttCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAAt                                                              <  1:681748/69‑1 (MQ=255)
ttCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAAt                                                              <  1:2713024/69‑1 (MQ=255)
ttCCGGCACGCTCAGGCAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAAt                                                              <  1:2954267/69‑1 (MQ=255)
 tCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATa                                                             <  1:2136456/69‑1 (MQ=255)
 tCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATa                                                             <  1:3144859/69‑1 (MQ=255)
  ccGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATa                                                             <  1:2699665/68‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:1128247/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:1099036/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:1824446/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:1809579/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:1716324/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:1620580/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:14598/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:576341/69‑1 (MQ=255)
               gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCAGGGTAGAACGCGGCAATAACTGCGGTATCAAc                                               <  1:2995742/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:2981779/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:2845334/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:2805585/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:3297834/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:3348679/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:402036/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:2493082/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:2647577/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:2013427/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:1711127/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:1769668/69‑1 (MQ=255)
                 aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                             <  1:1446477/69‑1 (MQ=255)
                  gTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCt                                             >  1:1642217/1‑68 (MQ=255)
                  gTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCt                                             >  1:18617/1‑68 (MQ=255)
                              tttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                >  1:3134356/1‑69 (MQ=255)
                              tttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                >  1:522012/1‑69 (MQ=255)
                               ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTca                               <  1:1943969/69‑1 (MQ=255)
                               ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTca                               <  1:2082039/69‑1 (MQ=255)
                               ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                >  1:3082179/1‑68 (MQ=255)
                               ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                >  1:3040105/1‑68 (MQ=255)
                               ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                >  1:1537764/1‑68 (MQ=255)
                               ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                >  1:2923365/1‑68 (MQ=255)
                                    gcgTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTAAACaaa                          <  1:4214/69‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:1196647/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:606306/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:619481/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:1465276/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                >  1:1144006/1‑42 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:684557/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:1940135/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:2875406/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:2092575/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                >  1:1784392/1‑42 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                >  1:2099435/1‑42 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:2364995/42‑1 (MQ=255)
                                         tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                <  1:2904076/42‑1 (MQ=255)
                                           aGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTc                   >  1:1505399/1‑69 (MQ=255)
                                             ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                 >  1:939371/1‑69 (MQ=255)
                                             ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                 >  1:2857288/1‑69 (MQ=255)
                                             ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                 >  1:279401/1‑69 (MQ=255)
                                             ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                 >  1:412944/1‑69 (MQ=255)
                                             ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                 >  1:558746/1‑69 (MQ=255)
                                             ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTAAg                 >  1:660698/1‑69 (MQ=255)
                                               aTTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               >  1:2629886/1‑69 (MQ=255)
                                               aTTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               >  1:1531582/1‑69 (MQ=255)
                                                       tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTTCTCAACAAAGTCGCTCAGtt               <  1:967186/61‑1 (MQ=255)
                                                       tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:1868856/61‑1 (MQ=255)
                                                       tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:564773/61‑1 (MQ=255)
                                                       tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:1898696/61‑1 (MQ=255)
                                                       tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:1449322/61‑1 (MQ=255)
                                                       tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:2730785/61‑1 (MQ=255)
                                                        aGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:3397203/60‑1 (MQ=255)
                                                        aGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt               <  1:2060273/60‑1 (MQ=255)
                                                            cgcgGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGa  <  1:183095/69‑1 (MQ=255)
                                                               |                                                                 
TTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGA  >  minE/2004889‑2005017

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: