Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,647,443 T→G 15.0% intergenic (+123/‑150) asd → / → glgB aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,647,4430TG15.0% 220.1 / 29.8 113intergenic (+123/‑150)asd/glgBaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (51/45);  new base G (7/10);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 113 total aligned reads displayed.

AGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGC  >  minE/2647380‑2647507
                                                               |                                                                
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:980173/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:749196/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:357030/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:3250644/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:2399360/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:2367653/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:2313565/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATa                                                             >  1:1687917/1‑69 (MQ=255)
aGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAACGGTGGCTGAATa                                                             >  1:1392643/1‑69 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:2430463/69‑1 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:578465/69‑1 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:548393/69‑1 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:2773068/69‑1 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:1997096/69‑1 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:1658248/69‑1 (MQ=255)
    cAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:137789/69‑1 (MQ=255)
     aGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATg                                                         <  1:2287680/68‑1 (MQ=255)
        tGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTa                                                     <  1:1919283/69‑1 (MQ=255)
        tGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTa                                                     <  1:2388023/69‑1 (MQ=255)
            tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                   >  1:741135/1‑67 (MQ=255)
            tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                   >  1:1410454/1‑67 (MQ=255)
            tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                   >  1:546999/1‑67 (MQ=255)
            tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                   >  1:505059/1‑67 (MQ=255)
            tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAtt                                                   >  1:494705/1‑67 (MQ=255)
            tCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTAt                                                    >  1:2980491/1‑66 (MQ=255)
                gATTTACCGTGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                             <  1:63403/69‑1 (MQ=255)
                gATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                             <  1:1693502/69‑1 (MQ=255)
                gATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                             <  1:2783638/69‑1 (MQ=255)
                gATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                             <  1:189035/69‑1 (MQ=255)
                  tttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                           >  1:1511177/1‑69 (MQ=255)
                  tttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                           >  1:2607754/1‑69 (MQ=255)
                  tttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCt                                           >  1:1226141/1‑69 (MQ=255)
                   ttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          <  1:3124220/69‑1 (MQ=255)
                   ttACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          <  1:1068668/69‑1 (MQ=255)
                    tACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          <  1:477158/68‑1 (MQ=255)
                        gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                     <  1:2502947/69‑1 (MQ=255)
                        gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                     <  1:118461/69‑1 (MQ=255)
                        gggAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTg                                     <  1:866998/69‑1 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGa                                    >  1:2151902/1‑60 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                           >  1:1214709/1‑69 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                           >  1:2035908/1‑69 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                           >  1:1707256/1‑69 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                           >  1:2224077/1‑69 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                           >  1:29352/1‑69 (MQ=255)
                                  tAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGgac                           >  1:1599283/1‑69 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:3321864/1‑67 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:272010/1‑67 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:2243395/1‑67 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:3283615/1‑67 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:542727/1‑67 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:1525475/1‑67 (MQ=255)
                                          tGGCTATTCTTTAAGGGTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCa                     >  1:1158496/1‑67 (MQ=255)
                                           ggCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagag                  >  1:2121939/1‑69 (MQ=255)
                                             cTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                             <  1:1386743/40‑1 (MQ=255)
                                                ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa              >  1:3233670/1‑68 (MQ=255)
                                                ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa              >  1:2720662/1‑68 (MQ=255)
                                                ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa              >  1:1518100/1‑68 (MQ=255)
                                                ttCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGa              >  1:2485125/1‑68 (MQ=255)
                                                 tCTTTAAGGGTGGCTGAATACATGAGTATTCACAGc                                             <  1:401763/36‑1 (MQ=255)
                                                  ctgtAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:1465274/4‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          >  1:1476779/1‑38 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          >  1:650529/1‑38 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          >  1:305415/1‑38 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          >  1:223759/1‑38 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCtt                                          >  1:2184567/1‑38 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:2466199/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:76871/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:1551166/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:685394/64‑1 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:1185947/64‑1 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:1681289/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:524567/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:1872992/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:461030/64‑1 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:1906661/64‑1 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:1391759/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:1218048/64‑1 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:2259811/1‑64 (MQ=255)
                                                  cTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                >  1:2978180/1‑64 (MQ=255)
                                                       aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca        <  1:408580/67‑1 (MQ=255)
                                                       aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca        <  1:2870167/67‑1 (MQ=255)
                                                       aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca        <  1:1420721/67‑1 (MQ=255)
                                                       aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca        <  1:161673/67‑1 (MQ=255)
                                                       aGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGcaca        <  1:624999/67‑1 (MQ=255)
                                                         ggTGGGGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:2978625/68‑1 (MQ=255)
                                                         ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:2294250/68‑1 (MQ=255)
                                                         ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:2750321/68‑1 (MQ=255)
                                                         ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:8437/68‑1 (MQ=255)
                                                         ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:817173/68‑1 (MQ=255)
                                                         ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:1189603/68‑1 (MQ=255)
                                                         ggTGGCGGAATACACGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACagag     <  1:2744820/68‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:3257593/56‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:3091221/56‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:1400197/56‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:3097602/56‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:1726701/56‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:321531/56‑1 (MQ=255)
                                                          gTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCagagag                <  1:2718412/56‑1 (MQ=255)
                                                           tGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgc  >  1:1320555/1‑69 (MQ=255)
                                                           tGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGtgc  >  1:1747538/1‑69 (MQ=255)
                                                               |                                                                
AGCGCAGATGTTTCATGATTTACCGGGAGTTAAATAGAGCATTGGCTATTCTTTAAGGGTGGCTGAATACATGAGTATTCACAGCCTTACCTGAAGTGAGGACGACGCAGAGAGGATGCACAGAGTGC  >  minE/2647380‑2647507

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: