Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,008,479 T→C 51.8% P181P (CCA→CCG hybO ← hydrogenase 2, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,008,4790TC51.8% 3.2 / 72.5 83P181P (CCA→CCGhybOhydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base T (37/3);  new base C (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 = 5.90e-01

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