Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,182,099 G→A 56.7% S222S (AGC→AGT yeaE ← predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,182,0990GA56.7% 21.7 / 70.0 90S222S (AGC→AGTyeaEpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base A (47/4);  total (69/21)
Fisher's exact test for biased strand distribution p-value = 1.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTGGCCGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCAT  >  minE/1182031‑1182154
                                                                    |                                                       
gTGGCCGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCcgcg                                                         >  1:3379011/1‑69 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACTCCAACAATACTTGTGCCGCGCTGa                                                     <  1:1219626/68‑1 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGa                                                     <  1:1015311/68‑1 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGa                                                     <  1:97159/68‑1 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGa                                                     <  1:3305028/68‑1 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGa                                                     <  1:3301035/68‑1 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGa                                                     <  1:1681913/68‑1 (MQ=255)
     cGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAACACTTGTGCCGCGCTGa                                                     <  1:3018766/68‑1 (MQ=255)
                tCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCAt                                         >  1:740929/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:3098987/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2707130/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2777391/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2822273/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2182983/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:655507/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2825852/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2880842/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2011405/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2931262/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2469043/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:76004/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:2962113/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:1216550/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:1568722/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:822879/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:1498737/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:878694/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:3336257/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:1458614/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:1416263/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:1375250/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:3336969/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:104812/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:944733/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCaa                                   >  1:3410226/1‑69 (MQ=255)
                      tCACACCCTGATGGCTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCa                                    >  1:41834/1‑68 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:2519285/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:2214367/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:2618102/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:2937524/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:2965320/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:2460942/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:1492857/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:918920/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:565678/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:1684448/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:1827685/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:686866/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:662101/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:618090/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:583037/1‑69 (MQ=255)
                        acacCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAAtt                                 >  1:132916/1‑69 (MQ=255)
                                  ggcTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGAc                         <  1:2074142/65‑1 (MQ=255)
                                  ggcTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGAc                         <  1:778585/65‑1 (MQ=255)
                                  ggcTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGAc                         <  1:782368/65‑1 (MQ=255)
                                  ggcTGATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGAc                         <  1:1576147/65‑1 (MQ=255)
                                  gaCTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTa                       >  1:3360401/1‑69 (MQ=255)
                                  gaCTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTa                       >  1:1219693/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:422723/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:909505/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:511373/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1015166/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:548558/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:960354/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1598507/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1094001/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1613059/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:3264033/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1740064/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1358801/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:2473201/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1472788/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1750411/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1476648/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:3420392/1‑69 (MQ=255)
                                      gATTACCCACGCCAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:1534049/1‑69 (MQ=255)
                                      gATTACCCACGACAACAATACTTGTGCCGCACTGATATTGTGAGCATGTGCAATTTCGTTGACTACcgc                   >  1:382444/1‑69 (MQ=255)
                                      gATAACCCACGCCAACAATACTTGTGCCGCACTGatat                                                  >  1:64642/1‑38 (MQ=38)
                                         cACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGc                                     >  1:2379668/1‑48 (MQ=255)
                                         cACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGc                                     >  1:1904559/1‑48 (MQ=255)
                                         cACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGc                                     >  1:1242873/1‑48 (MQ=255)
                                                gCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAAc         <  1:1226764/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:2025962/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:13183/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:73927/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:2430111/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:2695719/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:535693/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:3276648/69‑1 (MQ=255)
                                                    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTc     <  1:360270/69‑1 (MQ=255)
                                                       aTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCAt  <  1:2156010/69‑1 (MQ=255)
                                                                    |                                                       
GTGGCCGCTTTTGGAATCGCCATCACACCCTGATGACTGATCACCCACGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCAT  >  minE/1182031‑1182154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: