Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,008,482 C→T 51.8% L180L (CTG→CTA hybO ← hydrogenase 2, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,008,4820CT51.8% ‑2.2 / 75.1 83L180L (CTG→CTAhybOhydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base C (37/3);  new base T (24/19);  total (61/22)
Fisher's exact test for biased strand distribution p-value = 1.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GAGGAAGTTGTGCGGGTTCGGCGGGCAGCCCGGAATATTGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACC  >  minE/2008425‑2008538
                                                         |                                                        
gAGGAAGTTGTGCGGGTTCGGCGGGCAGCCCGGAATATTGATAACGGTTTTGCCTGGCAGAACTTCTTg                                               <  1:1897480/69‑1 (MQ=255)
gAGGAAGTTGTGCGGGTTCGGCGGGCAGCCCGGAATATTGATAACGGTTTTGCCTGGCAGAACTTCTTg                                               <  1:1901309/69‑1 (MQ=255)
gAGGAAGTTGTGCGGGTTCGGCGGGCAGCCCGGAATATTGATAACGGTTTTGCCTGGCAGAACTTCTTg                                               <  1:1370793/69‑1 (MQ=255)
                              aggAATATTGATATCGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2716013/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGACCTTCTTGCAGGCTGa                                       >  1:2932339/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2265294/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:3087551/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2847357/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2843458/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:278818/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:3308924/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2779446/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2725373/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2719676/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2356452/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2553821/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:2645189/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:1835236/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:3349371/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:823171/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:770057/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:612400/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:381929/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:1637732/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:426411/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:588223/2‑47 (MQ=255)
                              aggAATATTGATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGa                                       >  1:1315513/2‑47 (MQ=255)
                                  atatTGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGc              >  1:54026/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:963914/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:709885/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:478763/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:308409/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:35639/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:3243558/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:355814/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1061380/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:129223/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:2690779/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1950297/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1886122/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:17092/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1677261/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1522513/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:2023390/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1226696/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:1137748/1‑68 (MQ=255)
                                     ttGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGc           >  1:2534686/1‑68 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:54421/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:3209604/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:428708/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:1069579/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:905631/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:917890/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:3246664/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:1020980/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2330669/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2326164/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:3092540/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2299334/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2186187/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2826791/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2179964/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:29812/69‑1 (MQ=255)
                                       gATAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:2886562/69‑1 (MQ=255)
                                        aTAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:3092570/68‑1 (MQ=255)
                                        aTAACGGTTTTGCCCGGTAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGc        <  1:1139934/68‑1 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:463710/1‑69 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:998048/1‑69 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:2031726/1‑69 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:1415688/1‑69 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:3101713/1‑69 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:1726993/1‑69 (MQ=255)
                                          aaCGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAac     >  1:3105792/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCa                >  1:1714977/1‑55 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:299287/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:2993621/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:2013215/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:1084463/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:1965241/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:843527/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:2775354/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:3183930/1‑69 (MQ=255)
                                             ggTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACAcc  >  1:3168596/1‑69 (MQ=255)
                                                         |                                                        
GAGGAAGTTGTGCGGGTTCGGCGGGCAGCCCGGAATATTGATAACGGTTTTGCCTGGCAGAACTTCTTGCAGGCTGACTGCGCCAGTTGGGTTAACTCCAGCTGCGGCAACACC  >  minE/2008425‑2008538

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