Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 113,344 C→T 18.6% N231N (AAC→AAT lpd → lipoamide dehydrogenase, E3 component is part of three enzyme complexes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE113,3440CT18.6% 107.4 / 13.5 69N231N (AAC→AATlpdlipoamide dehydrogenase, E3 component is part of three enzyme complexes
Reads supporting (aligned to +/- strand):  ref base C (48/8);  new base T (5/8);  total (53/16)
Fisher's exact test for biased strand distribution p-value = 1.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACAAAGACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTATGTGAC  >  minE/113292‑113403
                                                    |                                                           
cTGACAAAGACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAAcc                                              <  1:1251862/68‑1 (MQ=255)
       aGACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTAc                                      >  1:2980249/1‑69 (MQ=255)
       aGACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTAc                                      >  1:2670490/1‑69 (MQ=255)
        gACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTAcc                                     <  1:725020/69‑1 (MQ=255)
        gACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTAcc                                     <  1:1975201/69‑1 (MQ=255)
              tttAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:2206612/2‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCt                                                           <  1:1436795/41‑1 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTAccg                                    >  1:1834962/1‑64 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:486261/1‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:2310481/1‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:848910/1‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:1909035/1‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:376101/1‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:2900557/1‑68 (MQ=255)
              gTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACAAAAGTTACCGACGt                                >  1:886006/1‑68 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:2487515/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:476139/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:2585660/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:1925869/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:196043/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:861812/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:1285394/1‑66 (MQ=255)
                      cTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:3009382/1‑66 (MQ=255)
                      cTTCAACAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGc                          >  1:875025/1‑66 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGaa                       >  1:2463671/1‑68 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGaa                       >  1:2702375/1‑68 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGaa                       >  1:2986121/1‑68 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGaa                       >  1:1723736/1‑68 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGaa                       >  1:2132664/1‑68 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGa                        >  1:2063176/1‑67 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGa                        >  1:1500379/1‑67 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTCGAAGCGaa                       >  1:2468111/1‑68 (MQ=255)
                       ttCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTAGAAGCGaa                       >  1:875666/1‑68 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:3085880/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:105359/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:1274571/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:3050142/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:1706712/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:2406753/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:239705/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaagaag                  >  1:1889795/1‑69 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaag                     <  1:1436987/66‑1 (MQ=255)
                           ccAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCCTTGAAGCGAaaga                    >  1:1324204/1‑67 (MQ=255)
                            cAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAaag                     <  1:1308768/65‑1 (MQ=255)
                                  tATCAGCAAGAAATTCAACCTGATGCTGGAAACAAAAGTTACCGCCGTTGAAGCGAaag                     <  1:1221560/59‑1 (MQ=255)
                                    tCAGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGt                                >  1:2331188/1‑46 (MQ=255)
                                    tCAGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:378339/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:380132/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:205067/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAtt          >  1:2955649/1‑68 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:2650171/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:1357588/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:880406/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:2014545/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:355935/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:2984255/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:540673/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAttt         >  1:756616/1‑69 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAtt          >  1:2011696/1‑68 (MQ=255)
                                    tCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCAtt          >  1:470651/1‑68 (MQ=255)
                                      aGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:232314/68‑1 (MQ=255)
                                      aGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:2154506/68‑1 (MQ=255)
                                      aGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:2202652/68‑1 (MQ=255)
                                      aGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:3192202/68‑1 (MQ=255)
                                      aGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:3071529/68‑1 (MQ=255)
                                      aGCAAGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:3058338/68‑1 (MQ=255)
                                          aGAAATTCAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:2568111/64‑1 (MQ=255)
                                           gAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTATGTGAc  <  1:1626657/69‑1 (MQ=255)
                                                 cAATCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTa        <  1:947668/57‑1 (MQ=255)
                                                    |                                                           
CTGACAAAGACATCGTTAAAGTCTTCACCAAGCGTATCAGCAAGAAATTCAACCTGATGCTGGAAACCAAAGTTACCGCCGTTGAAGCGAAAGAAGACGGCATTTATGTGAC  >  minE/113292‑113403

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: