Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 468,885 A→T 26.2% A303A (GCA→GCT sucB → dihydrolipoyltranssuccinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE468,8850AT26.2% 183.8 / 68.2 122A303A (GCA→GCTsucBdihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base A (49/41);  new base T (11/21);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 6.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

TGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACC  >  minE/468834‑468948
                                                   |                                                               
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:691850/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:1000384/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:2680755/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:214759/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:2801388/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:2852877/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:3002209/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:1520390/1‑68 (MQ=255)
tGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:3135939/1‑68 (MQ=255)
tGTCAATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCTGTTAAAGGCCGTGACg                                                 >  1:1808233/1‑68 (MQ=255)
  tCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGc                                               >  1:907060/1‑68 (MQ=255)
       tcccTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                         >  1:354706/2‑69 (MQ=255)
       aCCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                         >  1:959400/1‑69 (MQ=255)
       aCCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                         >  1:2385584/1‑69 (MQ=255)
       aCCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                         >  1:2381631/1‑69 (MQ=255)
       aCCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                         >  1:2433913/1‑69 (MQ=255)
       aCCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                         >  1:3235860/1‑69 (MQ=255)
          cTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGt                                                     >  1:2812673/1‑54 (MQ=255)
            cGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGt                                    <  1:263998/69‑1 (MQ=255)
            cGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGt                                    <  1:3409189/69‑1 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGAcc                                      >  1:2719872/1‑63 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGAcc                                      >  1:7170/1‑63 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGAcc                                      >  1:635981/1‑63 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGAcc                                      >  1:410009/1‑63 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGAcc                                      >  1:1960190/1‑63 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGAcc                                      >  1:1873606/1‑63 (MQ=255)
                aTGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGCCCGTGACGGCAAGCTGAc                                       >  1:1781561/1‑62 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTTACGGCa                                              >  1:2854849/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:3380813/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2772115/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              >  1:2724672/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              >  1:2704679/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2696252/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              >  1:619069/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              >  1:2882963/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2447632/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:3250020/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              >  1:2432898/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:828202/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2599002/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              >  1:2373117/1‑52 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2334865/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2029454/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2197469/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:1322510/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2201701/52‑1 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGAt                             >  1:943191/1‑69 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGAt                             >  1:1061835/1‑69 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGAt                             >  1:731965/1‑69 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGAt                             >  1:601214/1‑69 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGAt                             >  1:3335682/1‑69 (MQ=255)
                   gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTCACGGCa                                              <  1:1272039/52‑1 (MQ=255)
                    cAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCa                                              <  1:2643778/51‑1 (MQ=255)
                     aGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCa                                                             <  1:1095056/35‑1 (MQ=255)
                     aGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCa                                                             <  1:3099029/35‑1 (MQ=255)
                     aGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCa                                                             <  1:1824293/35‑1 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACtt            >  1:3050523/1‑67 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:1471688/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:3013890/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:3405461/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:1510229/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:721142/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:1727959/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:2870345/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:2441474/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:2251091/1‑68 (MQ=255)
                                      tCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtc           >  1:2243755/1‑68 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:2282806/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:1034059/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:1245343/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:2555404/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:1363721/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:194396/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:2407891/69‑1 (MQ=255)
                                        aaaGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTCACCGTTGAAGATCTGACCGGTGGTAACTtcacc        <  1:1229416/69‑1 (MQ=255)
                                          agagCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGt                                    <  1:1739757/39‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAATCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:1856544/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:1859582/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2769601/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2336190/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:932571/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:103853/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2525069/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2539807/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:414539/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:3404601/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2653212/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2733122/69‑1 (MQ=255)
                                              cTGGCTGTTAAAGGCCGTGACGGCAAGATGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:2901684/69‑1 (MQ=255)
                                               tGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:593063/68‑1 (MQ=255)
                                               tGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:195486/68‑1 (MQ=255)
                                               tGGCTGTTAAAGGCCGTGACGGCAAGCTGACGGTTGAAGATCTGACCGGTGGTAACTtcaccatcacc  <  1:1965936/68‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:1924962/62‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:1245295/62‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:2475712/62‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:2440764/62‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:833679/62‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:3360326/62‑1 (MQ=255)
                                                 gCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:2072261/62‑1 (MQ=255)
                                                 gCAGGCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTtcaccat      <  1:537265/62‑1 (MQ=255)
                                                   |                                                               
TGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACC  >  minE/468834‑468948

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: