Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 504,034 C→T 70.1% Y342Y (TAC→TAT modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE504,0340CT70.1% 73.4 / 36.8 80Y342Y (TAC→TATmodCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (9/15);  new base T (39/17);  total (48/32)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGC  >  minE/503980‑504102
                                                      |                                                                    
tATCAGCCCGTGGGCCAGGGATTAACTGGCTATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:998861/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTaaaa                                                         >  1:2431730/1‑68 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:2523485/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:700417/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:498776/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:2993509/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:1688043/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:1548977/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:143100/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:1391343/1‑69 (MQ=255)
tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAg                                                        >  1:1118009/1‑69 (MQ=255)
    aGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                    >  1:3000635/1‑69 (MQ=255)
      cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                  <  1:50659/69‑1 (MQ=255)
      cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                  <  1:2859924/69‑1 (MQ=255)
      cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                  <  1:451842/69‑1 (MQ=255)
      cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                  <  1:1187402/69‑1 (MQ=255)
      cccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                  <  1:2977815/69‑1 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGTCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:476987/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2474779/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:754001/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:1838531/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2684932/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2679441/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2545498/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:565089/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2450000/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2393483/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2341866/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:2879517/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:3127224/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:3308276/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:399573/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                       >  1:1358077/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTATGCGCAAATTAAAAGTCTGTCGATAACCGCCt                                       >  1:1773389/1‑69 (MQ=255)
                 gggATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                  <  1:3143036/58‑1 (MQ=255)
                   gATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                      >  1:11239/1‑68 (MQ=255)
                   gATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                      >  1:714100/1‑68 (MQ=255)
                   gATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                      >  1:2157743/1‑68 (MQ=255)
                      gAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                  <  1:1961791/69‑1 (MQ=255)
                      gAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                  <  1:168271/69‑1 (MQ=255)
                      gAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                  <  1:1970632/69‑1 (MQ=255)
                      gAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                  <  1:2013026/69‑1 (MQ=255)
                       aaCTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                 >  1:2266783/1‑69 (MQ=255)
                                 cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                 >  1:2931790/1‑59 (MQ=255)
                                 cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                 >  1:1599617/1‑59 (MQ=255)
                                 cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                 >  1:655455/1‑59 (MQ=255)
                                 cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                 >  1:2399178/1‑59 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTc                                                  >  1:3084927/1‑35 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATACCCGCCTGATTAAATCAGGTAGc                       >  1:1832403/1‑62 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:769909/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:3007367/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:2625911/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:1097661/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:1578078/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:2742835/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:3407090/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaaa                 >  1:2900669/1‑68 (MQ=255)
                                        ggCCTGTGGCTGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATaa                  >  1:2426430/1‑67 (MQ=255)
                                           cTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg               <  1:730289/67‑1 (MQ=255)
                                           cTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg               <  1:1750990/67‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAAGTGGGCa      <  1:2193938/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTggga       <  1:79775/68‑2 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:1521803/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:1425753/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:677404/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:680359/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:3260968/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:869912/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:3226460/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:3150022/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:924128/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:2718810/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:1965996/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:783028/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:2512398/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTGTGACGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:2074594/69‑1 (MQ=255)
                                                  tGTATGCGCAAATTAAAAGTATGTCGATAACCGCCTGATTAAATCAGGTAGCTATAAATGAACTGGGCa      <  1:308231/69‑1 (MQ=255)
                                                  tGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc       <  1:1656917/68‑1 (MQ=255)
                                                  tGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc       <  1:1493329/68‑1 (MQ=255)
                                                  tGTACGCGCAAATTAAAAGGGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc       <  1:1139105/68‑1 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1088767/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:938763/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:777737/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:791869/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1365890/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:980224/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:845355/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:56132/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:2816304/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1076296/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:2245751/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:2239018/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:2045327/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1942610/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1693942/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1144374/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1657840/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1613608/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1442991/1‑69 (MQ=255)
                                                      cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:1406161/1‑69 (MQ=255)
                                                      |                                                                    
TATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGC  >  minE/503980‑504102

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: