Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,500,204 A→G 39.9% S132S (TCT→TCC nuoB ← NADH:ubiquinone oxidoreductase, chain B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,500,2040AG39.9% 37.4 / 60.1 89S132S (TCT→TCCnuoBNADH:ubiquinone oxidoreductase, chain B
Reads supporting (aligned to +/- strand):  ref base A (22/31);  new base G (21/14);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TATACACATCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAC  >  minE/1500137‑1500265
                                                                   |                                                             
tataCACATCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCag                                                              <  1:775015/69‑1 (MQ=255)
tataCACATCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCag                                                              <  1:2695696/69‑1 (MQ=255)
    cacaTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:811938/1‑69 (MQ=255)
    cacaTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:11679/1‑69 (MQ=255)
    cacaTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:543947/1‑69 (MQ=255)
    cacaTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:825548/1‑69 (MQ=255)
    cacaTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:1690900/1‑69 (MQ=255)
    cacaTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:1561559/1‑69 (MQ=255)
    cacaACAACCGGGATGAATTTATCGACGCCATGCACAACGGAATAAATATCGTACATACCACCAGAGtt                                                          >  1:2031219/1‑69 (MQ=255)
      caTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        <  1:2485236/69‑1 (MQ=255)
      caTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        <  1:3419532/69‑1 (MQ=255)
      caTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        <  1:2097608/69‑1 (MQ=255)
       aTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGc                                                       <  1:213683/69‑1 (MQ=255)
       aTCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGc                                                       <  1:903588/69‑1 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:1498545/1‑67 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:2673453/1‑67 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:1697984/1‑67 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:174866/1‑67 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:2599284/1‑67 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:927608/1‑67 (MQ=255)
        tCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTg                                                         >  1:2389363/1‑66 (MQ=255)
          aaCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGcaca                                                    >  1:1198217/1‑69 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAGATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:3371597/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCTGAGTTGGcaca                                                    >  1:1671645/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1146993/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1490230/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:579264/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:2385277/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:877840/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:3269581/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:2650202/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:2665548/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1287482/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:2180922/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:2155061/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:895375/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1360806/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1387289/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:702060/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:719795/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1515003/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:996428/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1682109/1‑68 (MQ=255)
           aCCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGcaca                                                    >  1:1435288/1‑68 (MQ=255)
             cGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTgg                                                        >  1:362875/1‑62 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGTCAcc                                              <  1:3296583/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2906390/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:291113/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:1103971/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:292435/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:3018034/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:1471457/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:410024/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:471304/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:661863/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2747413/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2716665/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2664125/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2654355/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2600450/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2327899/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2568355/69‑1 (MQ=255)
                gATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:890208/69‑1 (MQ=255)
                 aTGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:1892969/68‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:994677/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:2802371/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:2023531/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:2102693/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:590756/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:233340/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:2387159/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:3432218/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:1235817/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:127144/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:2620625/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:2635908/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:31758/69‑1 (MQ=255)
                                cccTGCACAACGGAATAAATATCGTACATACCGCCGGAGTTGGCACAGGCACCCATTGAGATAACCCAt                              <  1:3006329/69‑1 (MQ=255)
                                    gCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTg                          >  1:999636/1‑69 (MQ=255)
                                    gCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTg                          >  1:599791/1‑69 (MQ=255)
                                    gCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTg                          >  1:3412331/1‑69 (MQ=255)
                                    gCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTg                          >  1:2103186/1‑69 (MQ=255)
                                            gAATAAATATCGTACATACCACCAGAGTTGGCACAGGCAcc                                              <  1:2570191/41‑1 (MQ=255)
                                                   tatCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGt           <  1:2264010/69‑1 (MQ=255)
                                                     tCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCa         <  1:1393710/69‑1 (MQ=255)
                                                            tACCACCAGAGTTGTCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAc  <  1:788980/69‑1 (MQ=255)
                                                            tACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGa   <  1:2512069/68‑1 (MQ=255)
                                                            tACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAc  >  1:8499/1‑69 (MQ=255)
                                                            tACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAc  <  1:2042890/69‑1 (MQ=255)
                                                            tACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAc  <  1:1887014/69‑1 (MQ=255)
                                                            tACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAc  >  1:111806/1‑69 (MQ=255)
                                                                   |                                                             
TATACACATCAACCGGGATGAATTTATCGACGCCCTGCACAACGGAATAAATATCGTACATACCACCAGAGTTGGCACAGGCACCCATTGAGATAACCCATTTTGGTTCCAGCATCTGGTCATACAGAC  >  minE/1500137‑1500265

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: