Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,497,667 G→A 47.3% G124G (GGC→GGT nuoE ← NADH:ubiquinone oxidoreductase, chain E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,497,6670GA47.3% 10.8 / 83.7 100G124G (GGC→GGTnuoENADH:ubiquinone oxidoreductase, chain E
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (22/25);  total (47/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATCATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACC  >  minE/1497602‑1497729
                                                                 |                                                              
tCGATCATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGc                                                               <  1:3118108/67‑1 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:967532/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:189207/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:1970548/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:1590477/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:261903/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:878309/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaa                                                         >  1:684912/1‑68 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:540062/68‑1 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:320884/68‑1 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:2432143/68‑1 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:1685023/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:1147377/1‑69 (MQ=255)
                                ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:600004/1‑69 (MQ=255)
                                ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:2562124/1‑69 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2431471/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2839287/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2804792/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2148497/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2180161/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2641393/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2889864/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3229037/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3236985/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3349765/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:406970/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1488758/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:448603/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:490298/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1456516/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:677137/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:727351/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:767449/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:850382/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1094676/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:290515/1‑68 (MQ=255)
                                        agcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTtc                      >  1:2187861/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:1424356/1‑54 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:975358/1‑54 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:1996598/1‑54 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTtttt                        >  1:858978/1‑53 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1475128/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1296784/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:866241/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:147423/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1462393/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:671304/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:14375/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:131063/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1140323/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1543217/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:921597/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:185666/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1090874/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:3101697/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2679694/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2173791/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1887750/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2887291/1‑55 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCTCCGCCTg        <  1:1814104/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:2215189/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:780585/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:52228/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:3402343/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:3144283/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:1884037/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:912095/69‑1 (MQ=255)
                                                           aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:1645941/69‑1 (MQ=255)
                                                           aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:963642/69‑1 (MQ=255)
                                                           aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:3052137/69‑1 (MQ=255)
                                                           aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATAcc  <  1:2824926/69‑1 (MQ=255)
                                                                 |                                                              
TCGATCATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACC  >  minE/1497602‑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: