Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,006,382 C→T 35.5% P181P (CCG→CCA 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,3820CT35.5% 64.0 / 51.1 90P181P (CCG→CCAhybBpredicted hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base T (18/14);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

ATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTC  >  minE/2006325‑2006448
                                                         |                                                                  
aTGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAgcgc                                                         <  1:1707761/69‑1 (MQ=255)
aTGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAgcgc                                                         <  1:2753730/69‑1 (MQ=255)
            gCCCGCCTAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCcat                                             >  1:1863398/1‑69 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                                    <  1:1345899/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                                    <  1:1616519/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                                    <  1:1644915/69‑1 (MQ=255)
                     gaTCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATgaagaa                                    <  1:1045145/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:2267352/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1095666/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:2208772/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:584893/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1985450/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:2911097/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1734151/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:2952348/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:2966846/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1360952/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1138461/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:899179/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:80504/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1445248/1‑69 (MQ=255)
                      atcatcAGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATAAAGAAc                                   >  1:1441565/1‑69 (MQ=255)
                      ataatcaGCGATCCCATTGAAGACTGGTGCATAGTTGGCAGCAGCGCACCAAGCGCGATGATaaagaa                                    >  1:2473503/4‑68 (MQ=255)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:3215600/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:3349989/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:80399/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:2300039/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:8186/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:3074979/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:3192870/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:287622/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:2901254/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:2018400/37‑1 (MQ=37)
                        catcaGCGATCCCATTGAAGACTGGTGCATAGTTGgc                                                                 <  1:795274/37‑1 (MQ=37)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:2595053/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:1510496/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:1193713/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:1836348/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:596133/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:1918847/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:1226388/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:328097/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             <  1:402178/69‑1 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:258815/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:197059/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:2506644/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:591350/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:2323968/1‑69 (MQ=255)
                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCAcc                             >  1:2160748/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:301749/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:635449/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:2069024/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:2655817/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:3314243/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:3353415/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:2179907/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:405202/67‑1 (MQ=255)
                                                gTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGc           <  1:474401/67‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca       <  1:2140330/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca       <  1:1988584/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca       <  1:563337/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca       <  1:2420345/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGaca       <  1:3263811/69‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac        <  1:2255329/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac        <  1:986471/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac        <  1:1314024/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac        <  1:3196616/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac        <  1:1857192/68‑1 (MQ=255)
                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac        <  1:2806253/68‑1 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:1837560/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:1459064/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:1457206/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:523858/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:123812/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:1234585/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:2824776/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:952732/1‑69 (MQ=255)
                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTc  >  1:453550/1‑69 (MQ=255)
                                                         |                                                                  
ATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTC  >  minE/2006325‑2006448

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: