Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,488,482 G→A 63.7% G182G (GGC→GGT nuoM ← NADH:ubiquinone oxidoreductase, membrane subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,488,4820GA63.7% 89.8 / 92.7 146G182G (GGC→GGTnuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base G (31/22);  new base A (48/45);  total (79/67)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 162 total aligned reads displayed.

GTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTT  >  minE/1488437‑1488549
                                             |                                                                   
gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                                               >  1:1021568/1‑68 (MQ=255)
gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                                               >  1:1831518/1‑68 (MQ=255)
gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                                               >  1:414919/1‑68 (MQ=255)
  aGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATgaag                                             <  1:2225037/68‑1 (MQ=255)
   gTGAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACT‑‑CCTGGGTGTAAATGAAGAAc                                          <  1:405860/68‑1 (MQ=255)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:13136/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:1894943/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:592725/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:147032/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:1410425/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:1409243/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:2564492/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:2858643/68‑1 (MQ=39)
     gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                                          <  1:2981725/68‑1 (MQ=39)
      aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACtt                                        >  1:738888/1‑69 (MQ=255)
      aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACtt                                        >  1:2726750/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTgg                                                         >  1:463226/1‑47 (MQ=37)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:2273504/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:3013666/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:2304703/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:2379688/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:381033/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:1851817/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:1790918/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:6562/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:1084670/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:1746668/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCCCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:1562615/1‑69 (MQ=255)
           aaaCCAGCGCCAGGATGGCGATCAACATCACAAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                                   >  1:730549/1‑69 (MQ=255)
              ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg                                >  1:1287423/1‑69 (MQ=255)
              ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg                                >  1:1000563/1‑69 (MQ=255)
              ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg                                >  1:2835546/1‑69 (MQ=255)
              ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgcc                                 >  1:936533/1‑68 (MQ=255)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:3389345/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:2801593/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:1126520/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:3378263/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:1186067/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:2479723/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:3031068/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:2936488/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:764556/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:2255036/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:1631176/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:1506893/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:1601037/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:1914240/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:682596/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:561947/55‑1 (MQ=38)
                        gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                    <  1:215467/55‑1 (MQ=38)
                        gATGGCGAGCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGAtt                                    <  1:2309748/55‑1 (MQ=37)
                          tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcgcg                    >  1:2749390/1‑69 (MQ=255)
                          tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcgcg                    >  1:1972794/1‑69 (MQ=255)
                          tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGCGATAcgcg                    >  1:1612076/1‑69 (MQ=255)
                            gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGa                          <  1:129519/61‑1 (MQ=255)
                                tCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAc              <  1:541461/69‑1 (MQ=255)
                                tCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAc              <  1:1940732/69‑1 (MQ=255)
                                tCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAc              <  1:2542963/69‑1 (MQ=255)
                                tCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAc              <  1:806771/69‑1 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGGCGTGATACGCGt                   >  1:867784/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:356774/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:496623/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:3249244/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:266813/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:120453/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:2102240/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:1960801/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:2200817/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:1230529/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:2342133/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:2601122/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:2879555/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:2898114/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:1247674/1‑64 (MQ=255)
                                tCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                   >  1:1704808/1‑64 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:2835625/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:3400974/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:912336/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             <  1:3263674/69‑1 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:845129/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:820578/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:2195086/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTAcc             >  1:937421/1‑69 (MQ=255)
                                 cAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGATTTAcc             >  1:1263246/1‑69 (MQ=255)
                                          cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg     <  1:2920091/68‑1 (MQ=255)
                                          cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg     <  1:1107659/68‑1 (MQ=255)
                                          cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg     <  1:655450/68‑1 (MQ=255)
                                          cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg     <  1:1562238/68‑1 (MQ=255)
                                          cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTAGTTGCCGCCGTGATACGCGTTTTACCGTCAGAgg     <  1:966231/68‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:2127005/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:1443244/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:6074/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:439724/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:439464/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:225216/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:2697603/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:3428985/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:3171229/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:3096052/69‑1 (MQ=255)
                                            ggCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt  <  1:3069322/69‑1 (MQ=255)
                                             |                                                                   
GTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTT  >  minE/1488437‑1488549

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: