Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,004,964 G→A 26.7% T262T (ACC→ACT hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,004,9640GA26.7% 75.5 / 23.8 67T262T (ACC→ACThybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base G (14/35);  new base A (6/12);  total (20/47)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

GTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTA  >  minE/2004904‑2005026
                                                            |                                                              
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:1446477/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:2805585/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:2647577/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:2845334/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:2493082/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:2981779/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:2013427/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:3297834/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:1769668/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:402036/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:3348679/69‑1 (MQ=255)
  aGTTCACCGCACCTTTACCGCGCGTCAGCCATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCt                                                      <  1:1711127/69‑1 (MQ=255)
   gTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCt                                                      >  1:1642217/1‑68 (MQ=255)
   gTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCt                                                      >  1:18617/1‑68 (MQ=255)
               tttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                         >  1:522012/1‑69 (MQ=255)
               tttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                         >  1:3134356/1‑69 (MQ=255)
                ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTca                                        <  1:1943969/69‑1 (MQ=255)
                ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTca                                        <  1:2082039/69‑1 (MQ=255)
                ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                         >  1:3040105/1‑68 (MQ=255)
                ttACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                         >  1:3082179/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:2904076/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:2875406/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:2364995/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:684557/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:619481/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:606306/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:1465276/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         >  1:1784392/1‑42 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         >  1:2099435/1‑42 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:2092575/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:1196647/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         <  1:1940135/42‑1 (MQ=255)
                          tCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCaa                                                         >  1:1144006/1‑42 (MQ=255)
                            aGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTc                            >  1:1505399/1‑69 (MQ=255)
                              ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                          >  1:558746/1‑69 (MQ=255)
                              ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAg                          >  1:412944/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)
                              ccATTCCGGGTAGAACGCAGCGATAACTGCAGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTAAg                          >  1:660698/1‑69 (MQ=255)
                                aTTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        >  1:1531582/1‑69 (MQ=255)
                                aTTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        >  1:2629886/1‑69 (MQ=255)
                                        tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTTCTCAACAAAGTCGCTCAGtt                        <  1:967186/61‑1 (MQ=255)
                                        tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        <  1:564773/61‑1 (MQ=255)
                                        tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        <  1:2730785/61‑1 (MQ=255)
                                        tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        <  1:1449322/61‑1 (MQ=255)
                                        tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        <  1:1868856/61‑1 (MQ=255)
                                        tAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGtt                        <  1:1898696/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)
                                                    aTAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATg    <  1:2137606/69‑1 (MQ=255)
                                                      aaCTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTa  <  1:1321820/69‑1 (MQ=255)
                                                      aaCTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTa  <  1:1163595/69‑1 (MQ=255)
                                                      aaCTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTa  <  1:623824/69‑1 (MQ=255)
                                                            |                                                              
GTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTA  >  minE/2004904‑2005026

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: