Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,497,670 A→G 47.5% D123D (GAT→GAC nuoE ← NADH:ubiquinone oxidoreductase, chain E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,497,6700AG47.5% 3.3 / 89.0 99D123D (GAT→GACnuoENADH:ubiquinone oxidoreductase, chain E
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (22/25);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACC  >  minE/1497607‑1497729
                                                               |                                                           
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:967532/1‑69 (MQ=255)
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:878309/1‑69 (MQ=255)
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:261903/1‑69 (MQ=255)
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:1970548/1‑69 (MQ=255)
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:189207/1‑69 (MQ=255)
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:1590477/1‑69 (MQ=255)
catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaa                                                         >  1:684912/1‑68 (MQ=255)
                    cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:1685023/68‑1 (MQ=255)
                    cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:2432143/68‑1 (MQ=255)
                    cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:320884/68‑1 (MQ=255)
                    cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:540062/68‑1 (MQ=255)
                     acaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGatttc                                    >  1:2181581/1‑66 (MQ=255)
                     acaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGatttc                                    >  1:2848761/1‑66 (MQ=255)
                          tCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCa                              >  1:2732552/1‑69 (MQ=255)
                          tCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTATGCCCTGGTTTGATGTTc                               >  1:787459/1‑68 (MQ=255)
                           ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCa                              >  1:2204880/1‑68 (MQ=255)
                           ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:600004/1‑69 (MQ=255)
                           ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:2562124/1‑69 (MQ=255)
                           ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:1147377/1‑69 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2804792/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2431471/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1094676/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1456516/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:677137/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:850382/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:490298/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2180161/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2641393/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2148497/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:448603/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3349765/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2839287/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1488758/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:406970/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:727351/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2889864/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:767449/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3229037/68‑3 (MQ=255)
                              cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3236985/68‑3 (MQ=255)
                               aGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGGCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3341713/67‑3 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGTCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:987417/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:3397819/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:1072642/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:935599/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:1361994/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:1486282/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:3379001/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:297614/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:2904740/68‑1 (MQ=255)
                                  cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:2048088/68‑1 (MQ=255)
                                   agcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCGCTGGTTTGATGTTCAGCTTTTtc                      >  1:3110327/1‑68 (MQ=255)
                                   agcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       >  1:1805057/1‑67 (MQ=255)
                                   agcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTtc                      >  1:1989886/1‑68 (MQ=255)
                                   agcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTtc                      >  1:2187861/1‑68 (MQ=255)
                                   agcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTtc                      >  1:290515/1‑68 (MQ=255)
                                   agcaAGTTGGAAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTtc                      >  1:517591/1‑68 (MQ=255)
                                      aaGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:1849653/50‑1 (MQ=255)
                                          tGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:925892/56‑3 (MQ=255)
                                               gcagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCgccgcc          >  1:971651/1‑68 (MQ=255)
                                               gcagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCgccgcc          >  1:3275549/1‑68 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2583933/50‑3 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:1996598/1‑54 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:1424356/1‑54 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:975358/1‑54 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTtttt                        >  1:858978/1‑53 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:131063/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1296784/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:866241/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1140323/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:921597/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1090874/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:14375/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:671304/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1462393/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:147423/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1475128/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1543217/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1887750/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2887291/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:185666/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2679694/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2173791/1‑55 (MQ=255)
                                                cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:3101697/1‑55 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCTCCGCCTg        <  1:1814104/69‑1 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:2215189/69‑1 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:52228/69‑1 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:1884037/69‑1 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:3144283/69‑1 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:780585/69‑1 (MQ=255)
                                                cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:3402343/69‑1 (MQ=255)
                                                 agcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1141257/49‑3 (MQ=255)
                                                  gcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2606248/48‑3 (MQ=255)
                                                    gTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATa    >  1:220777/1‑69 (MQ=255)
                                                    gTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATa    >  1:2226820/1‑69 (MQ=255)
                                                      aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:3052137/69‑1 (MQ=255)
                                                      aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:912095/69‑1 (MQ=255)
                                                      aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:1645941/69‑1 (MQ=255)
                                                      aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:963642/69‑1 (MQ=255)
                                                      aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:2824926/69‑1 (MQ=255)
                                                               |                                                           
CATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACC  >  minE/1497607‑1497729

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: