Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,497,661 A→G 49.5% F126F (TTT→TTC nuoE ← NADH:ubiquinone oxidoreductase, chain E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,497,6610AG49.5% ‑5.9 / 88.8 95F126F (TTT→TTCnuoENADH:ubiquinone oxidoreductase, chain E
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (22/25);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 101 total aligned reads displayed.

TCGATCATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACC  >  minE/1497602‑1497729
                                                           |                                                                    
tCGATCATCATGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGc                                                               <  1:3118108/67‑1 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTaa                                                                     >  1:3093417/1‑56 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:1590477/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:967532/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:261903/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:189207/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:878309/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaaa                                                        >  1:1970548/1‑69 (MQ=255)
     catcatGTTTGGCCCTTTATCACAGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCaa                                                         >  1:684912/1‑68 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:1685023/68‑1 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:540062/68‑1 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:320884/68‑1 (MQ=255)
                         cacaGTTCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTg                                     <  1:2432143/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:2562124/1‑69 (MQ=255)
                                ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:600004/1‑69 (MQ=255)
                                ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAg                             >  1:1147377/1‑69 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3229037/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1094676/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2431471/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3349765/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2641393/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:850382/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1456516/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2180161/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:406970/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2148497/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2804792/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:448603/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3236985/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:767449/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:727351/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:490298/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2839287/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2889864/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:1488758/68‑3 (MQ=255)
                                   cAGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:677137/68‑3 (MQ=255)
                                    aGGCAGCAAGTTGGCAGCAGCGTGAAGCGACCGGCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:3341713/67‑3 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:2904740/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:297614/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:1072642/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:3379001/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:3397819/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:935599/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:1361994/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:2048088/68‑1 (MQ=255)
                                       cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCttttt                       <  1:1486282/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:3275549/1‑68 (MQ=255)
                                                    gcagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCgccgcc          >  1:971651/1‑68 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAgtt                           <  1:2583933/50‑3 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:975358/1‑54 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:1424356/1‑54 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGtttttt                       >  1:1996598/1‑54 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTtttt                        >  1:858978/1‑53 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:671304/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1090874/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1140323/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1296784/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:131063/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:921597/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:866241/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:14375/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:3101697/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:185666/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2887291/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2679694/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:1887750/1‑55 (MQ=255)
                                                     cagcGTGAAGCGACCGTCAAATGTCGTTTGCCCTGGTTTGATGTTCAGTTTTTtc                      >  1:2173791/1‑55 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCTCCGCCTg        <  1:1814104/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:3144283/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:2215189/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:3402343/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:780585/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:1884037/69‑1 (MQ=255)
                                                     cagcGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTg        <  1:52228/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:2824926/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)
                                                           |                                                                    
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: