Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,006,358 C→T 57.7% G189G (GGG→GGA hybB ← predicted hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,006,3580CT57.7% 24.9 / 63.7 85G189G (GGG→GGAhybBpredicted hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base T (35/14);  total (51/34)
Fisher's exact test for biased strand distribution p-value = 1.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CAGCGGCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACC  >  minE/2006292‑2006421
                                                                  |                                                               
cAGCGGCAACGTTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGAccc                                                               <  1:1931734/69‑1 (MQ=255)
cAGCGGCAACATTTCATAGCTCTGTCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGAccc                                                               <  1:1664553/69‑1 (MQ=255)
 aGCGGCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGAccc                                                               <  1:1130094/68‑1 (MQ=255)
 aGCGGCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGAccc                                                               <  1:3010105/68‑1 (MQ=255)
       aaCATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGaa                                                        <  1:431581/69‑1 (MQ=255)
       aaCATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGaa                                                        <  1:1358944/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGTTCATTAGCGACCCCATTGAAGACTgg                                                  <  1:942344/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:2921011/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:229705/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:2833842/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:1251173/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:941404/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:982879/69‑1 (MQ=255)
             tCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:2571759/69‑1 (MQ=255)
              cATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTgg                                                  <  1:1427328/68‑1 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCTAGATCATCAGCGATCCCa                                                             >  1:2794896/1‑38 (MQ=38)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGAtccc                                                              >  1:769297/1‑37 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:2185736/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1787661/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1835167/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:416548/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:404664/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:401917/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:3220663/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:2722577/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:2585026/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1109485/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:679001/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1487700/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:689878/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1288531/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1294363/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGATCCCa                                                             >  1:1393433/1‑38 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAgcgc                              <  1:1707761/69‑1 (MQ=255)
                                 aTGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAgcgc                              <  1:2753730/69‑1 (MQ=255)
                                             gCCCGCCTAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCcat                  >  1:1863398/1‑69 (MQ=255)
                                                      gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa         <  1:2961699/69‑1 (MQ=255)
                                                      gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa         <  1:1045145/69‑1 (MQ=255)
                                                      gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa         <  1:1345899/69‑1 (MQ=255)
                                                      gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa         <  1:1616519/69‑1 (MQ=255)
                                                      gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa         <  1:1644915/69‑1 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGgcag                                    >  1:1777331/1‑41 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATaaagaa         >  1:2298012/1‑68 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:2952348/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:899179/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:2911097/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:2966846/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:80504/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:584893/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1095666/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1138461/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1734151/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1445248/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1441565/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:2267352/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:2208772/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1360952/1‑69 (MQ=255)
                                                       atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc        >  1:1985450/1‑69 (MQ=255)
                                                       ataatcaGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATaaagaa         >  1:2473503/4‑68 (MQ=255)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:3192870/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:157013/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:477388/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:790864/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:795274/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:3349989/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:80399/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:8186/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:3215600/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:3074979/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:2018400/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:2300039/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:2901254/37‑1 (MQ=37)
                                                         catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                      <  1:287622/37‑1 (MQ=37)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  <  1:402178/69‑1 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:1226388/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:1193713/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:1510496/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:596133/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:591350/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:1836348/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:1918847/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:2595053/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:197059/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:328097/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:2160748/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:2323968/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:2506644/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc  >  1:258815/1‑69 (MQ=255)
                                                             aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAc   >  1:3317296/1‑68 (MQ=255)
                                                              gCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATtaagaa         <  1:1412031/61‑1 (MQ=255)
                                                                  |                                                               
CAGCGGCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACC  >  minE/2006292‑2006421

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: