Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,006,379 G→A 39.4% T182T (ACC→ACT 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,3790GA39.4% 40.9 / 52.2 81T182T (ACC→ACThybBpredicted hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base G (20/29);  new base A (18/14);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACA  >  minE/2006325‑2006443
                                                      |                                                                
aTGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAgcgc                                                    <  1:2753730/69‑1 (MQ=255)
aTGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAgcgc                                                    <  1:1707761/69‑1 (MQ=255)
            gCCCGCCTAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCcat                                        >  1:1863398/1‑69 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                               <  1:1045145/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                               <  1:1644915/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                               <  1:1616519/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                               <  1:1345899/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                               <  1:2961699/69‑1 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGgcag                                                          >  1:1777331/1‑41 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATaaagaa                               >  1:2298012/1‑68 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1985450/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:2952348/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:2966846/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:2911097/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:80504/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:2208772/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:2267352/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:584893/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1360952/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1095666/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1138461/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:899179/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1441565/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1734151/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                              >  1:1445248/1‑69 (MQ=255)
                      ataatcaGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATaaagaa                               >  1:2473503/4‑68 (MQ=255)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:2901254/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:157013/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:477388/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:287622/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:2300039/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:3349989/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:3215600/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:3192870/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:8186/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:790864/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:795274/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:2018400/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:80399/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                            <  1:3074979/37‑1 (MQ=37)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:328097/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        <  1:402178/69‑1 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:258815/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:1193713/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:1226388/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:1510496/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:1836348/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:1918847/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:197059/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:2160748/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:596133/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:2323968/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:591350/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:2506644/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                        >  1:2595053/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAc                         >  1:3317296/1‑68 (MQ=255)
                             gCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATtaagaa                               <  1:1412031/61‑1 (MQ=255)
                                        gAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGc            >  1:656510/1‑69 (MQ=255)
                                        gAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGc            >  1:25827/1‑69 (MQ=255)
                                        gAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGc            >  1:1155155/1‑69 (MQ=255)
                                        gAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGc            >  1:3046225/1‑69 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:3314243/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:3353415/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:301749/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:405202/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:474401/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:2069024/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:2655817/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:2179907/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc      <  1:635449/67‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca  <  1:2420345/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca  <  1:1988584/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca  <  1:2140330/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca  <  1:3263811/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca  <  1:563337/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac   <  1:2255329/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac   <  1:2806253/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac   <  1:1857192/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac   <  1:1314024/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac   <  1:3196616/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac   <  1:986471/68‑1 (MQ=255)
                                                      |                                                                
ATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACA  >  minE/2006325‑2006443

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: