Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 37,432 G→A 29.8% T211T (ACC→ACT caiA ← crotonobetaine reductase subunit II, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE37,4320GA29.8% 108.8 / 47.8 107T211T (ACC→ACTcaiAcrotonobetaine reductase subunit II, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (31/44);  new base A (8/24);  total (39/68)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 110 total aligned reads displayed.

TGACGCGGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGG  >  minE/37365‑37497
                                                                   |                                                                 
tGACGCGGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAg                                                                   >  1:2944513/1‑68 (MQ=255)
tGACGCGGTTAAAGCCGTCACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAg                                                                   >  1:101666/1‑68 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:2327954/68‑1 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:2189252/68‑1 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:2165581/68‑1 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:1601182/68‑1 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:356851/68‑1 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:1187588/68‑1 (MQ=255)
   cgcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:565332/68‑1 (MQ=255)
    gcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:1078281/67‑1 (MQ=255)
    gcgGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTg                                                                <  1:1047134/67‑1 (MQ=255)
         tAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTcac                                                         >  1:426293/1‑69 (MQ=255)
         tAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTcac                                                         >  1:2144217/1‑69 (MQ=255)
         tAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTcac                                                         >  1:2782151/1‑69 (MQ=255)
         tAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTcac                                                         >  1:1576008/1‑69 (MQ=255)
         tAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTcac                                                         >  1:3154593/1‑69 (MQ=255)
         tAAAGCCCTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTcac                                                         >  1:1569869/1‑69 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:444650/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:2761992/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:253246/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:2411992/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:1632656/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:3264355/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:38832/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:1301974/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:1271663/69‑1 (MQ=255)
                gTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:1267153/69‑1 (MQ=255)
                 ttACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:426277/68‑1 (MQ=255)
                 ttACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCt                                                  <  1:286080/68‑1 (MQ=255)
                   aCCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATc                                               <  1:1631394/69‑1 (MQ=255)
                   aCCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATc                                               <  1:2755165/69‑1 (MQ=255)
                   aCCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATc                                               >  1:2131286/1‑69 (MQ=255)
                      ttCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCa                                             <  1:1043098/68‑1 (MQ=255)
                      ttCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCa                                             <  1:265709/68‑1 (MQ=255)
                      ttCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCa                                             <  1:9419/68‑1 (MQ=255)
                      ttCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCa                                             <  1:3090367/68‑1 (MQ=255)
                      ttCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCa                                             <  1:3249974/68‑1 (MQ=255)
                      ttCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCa                                             <  1:50306/68‑1 (MQ=255)
                                    tGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCt                              >  1:1655793/1‑69 (MQ=255)
                                    tGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCt                              >  1:1768941/1‑69 (MQ=255)
                                    tGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCt                              >  1:3017195/1‑69 (MQ=255)
                                    tGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCt                              >  1:2940352/1‑69 (MQ=255)
                                    tGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCt                              >  1:281711/1‑69 (MQ=255)
                                      tCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATc                                               >  1:961263/1‑50 (MQ=255)
                                      tCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCtt                             >  1:1179893/1‑68 (MQ=255)
                                      tCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCtt                             >  1:2739384/1‑68 (MQ=255)
                                      tCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCtt                             >  1:2708981/1‑68 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCttt                            >  1:2226269/1‑67 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCttt                            >  1:750660/1‑67 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCttt                            >  1:2911759/1‑67 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCttt                            >  1:2420628/1‑67 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCttt                            >  1:1603682/1‑67 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCttt                            >  1:2770404/1‑67 (MQ=255)
                                        tttCTCGTCCAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCt                              >  1:840955/1‑65 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTGACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1160932/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:3368508/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:3311848/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:328709/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1385630/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1127574/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1351664/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:3225633/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:2290164/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:507732/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:896166/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:47650/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1845813/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1922517/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:1963319/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:2068281/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:2964624/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:2072325/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:786625/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:720287/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:755419/68‑1 (MQ=255)
                                                ccAGTTCCACGTCGTCAAAAGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTg                   <  1:61714/68‑1 (MQ=255)
                                                     tCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:839973/1‑69 (MQ=255)
                                                     tCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:1203361/1‑69 (MQ=255)
                                                     tCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:2150578/1‑69 (MQ=255)
                                                     tCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:5804/1‑69 (MQ=255)
                                                      ccACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGa                                     >  1:3300899/1‑44 (MQ=255)
                                                      ccACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:3424906/1‑68 (MQ=255)
                                                      ccACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:324876/1‑68 (MQ=255)
                                                      ccACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACt             >  1:288898/1‑68 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:1666485/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:753535/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:1590749/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2770957/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2673906/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2665358/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2559783/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2445440/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2420517/68‑1 (MQ=255)
                                                        aCGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2389828/68‑1 (MQ=255)
                                                         cgtcgtCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2332457/67‑1 (MQ=255)
                                                              tcAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACttt           <  1:2444106/62‑1 (MQ=255)
                                                                aaaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCgg  >  1:1230968/1‑69 (MQ=255)
                                                                aaaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCgg  >  1:2598965/1‑69 (MQ=255)
                                                                aaaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCgg  >  1:2246102/1‑69 (MQ=255)
                                                                aaaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCgg  >  1:2396291/1‑69 (MQ=255)
                                                                   |                                                                 
TGACGCGGTTAAAGCCGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGG  >  minE/37365‑37497

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: