Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,433,227 T→C 38.1% L360L (CTA→CTG mqo ← malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,433,2270TC38.1% 37.1 / 48.8 76L360L (CTA→CTGmqomalate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (28/19);  new base C (16/13);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAC  >  minE/1433170‑1433292
                                                         |                                                                 
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:1551264/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:2095028/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:2284516/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:19870/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:2418889/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:1788731/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:1621809/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:812980/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:155792/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:3388460/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:1488436/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:378131/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:1158938/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:117150/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:1372145/1‑68 (MQ=255)
aTTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCcaca                                                         >  1:93039/1‑68 (MQ=255)
                  catcatCGGCATCACGTTAGAGGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGaa                                      <  1:1475709/69‑1 (MQ=255)
                  catcatCGGCATCACGTTAGAGGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGaa                                      <  1:3166572/69‑1 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:824609/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:2184293/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:2338106/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:894268/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:195407/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:567214/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:415430/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:490371/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:1409279/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:1514857/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:1441139/68‑3 (MQ=255)
                        cGGCATCACGTTAGAGGTGGTGGTGGAACTCATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:429851/68‑3 (MQ=255)
                                      agGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCa                   <  1:2503205/68‑1 (MQ=255)
                                      agGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCa                   <  1:2763331/68‑1 (MQ=255)
                                      agGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCa                   <  1:1074684/68‑1 (MQ=255)
                                      agGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCa                   <  1:660577/68‑1 (MQ=255)
                                      agGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCa                   <  1:1473788/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:2080293/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:2006934/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:777308/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:1914602/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:1258431/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:1446664/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:1707543/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:2987222/68‑1 (MQ=255)
                                       ggtggtggtGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCaa                  <  1:3007740/68‑1 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCTCAATGAACCGTTTTTGAGGAATTTGGTTGAGAATGTGGCAAAt                >  1:2848542/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGTAAAt                >  1:1513510/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:853172/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:772644/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:577587/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:1459096/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:2403129/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:2489773/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:1744135/1‑69 (MQ=255)
                                        gtggtggtgGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAAt                >  1:2832202/1‑69 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:407241/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:1849097/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:3032728/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:2898328/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:3355456/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:1567940/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:3233967/1‑68 (MQ=255)
                                           gtggtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATgg              >  1:822255/1‑68 (MQ=255)
                                              gtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCa          >  1:2755994/1‑69 (MQ=255)
                                              gtggAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCa          >  1:1608314/1‑69 (MQ=255)
                                                  aaCTAATCAGATCCCACAATGAACCGTTTTTGAGGAAtttag                                 <  1:2952085/42‑3 (MQ=38)
                                                   aCTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTTGTTGAGAAGGTGGCAAATGGCCCAAACAg     <  1:897819/69‑1 (MQ=255)
                                                   aCTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAg     >  1:1101436/1‑69 (MQ=255)
                                                   aCTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAg     <  1:3349037/69‑1 (MQ=255)
                                                   aCTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAg     <  1:3179364/69‑1 (MQ=255)
                                                      cATTAGATCCCTCAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAc  >  1:1965585/1‑69 (MQ=255)
                                                      cATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAACCAGCAc  >  1:97800/1‑69 (MQ=255)
                                                      cATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAc  >  1:2727975/1‑69 (MQ=255)
                                                      cATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAc  >  1:2031334/1‑69 (MQ=255)
                                                      cATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAc  >  1:1374488/1‑69 (MQ=255)
                                                      cATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAc  >  1:1283397/1‑69 (MQ=255)
                                                      cATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAc  >  1:959421/1‑69 (MQ=255)
                                                         |                                                                 
ATTATCCAGCCCGACGTGCATCATCGGCATCACGTTAGAGGTGGTGGTGGAACTCATTAGATCCCACAATGAACCGTTTTTGAGGAATTTGGTTGAGAAGGTGGCAAATGGCCCAAACAGCAC  >  minE/1433170‑1433292

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: