Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,046 A→G 21.3% N341N (AAT→AAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0460AG21.3% 180.4 / 43.6 94N341N (AAT→AACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (44/30);  new base G (8/12);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAG  >  minE/1014978‑1015115
                                                                    |                                                                     
ttAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                                                      >  1:1242021/1‑70 (MQ=255)
 tAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGttg                                                                    >  1:515249/1‑70 (MQ=255)
 tAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGttg                                                                    >  1:4670912/1‑70 (MQ=255)
 tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                                                    >  1:4273965/1‑71 (MQ=255)
  aaGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGttgt                                                                   >  1:836348/1‑69 (MQ=255)
  aaGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTAt                                                                   >  1:3408586/1‑71 (MQ=255)
        aaCGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTc                                                               >  1:4136100/1‑69 (MQ=255)
           gcgcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtt                                                          >  1:1228562/1‑71 (MQ=255)
            ttcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                                                         >  1:3514730/3‑71 (MQ=255)
             tcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                                                         <  1:2587755/69‑1 (MQ=255)
             tcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                                                         <  1:3939654/69‑1 (MQ=255)
             gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                        <  1:1686862/71‑1 (MQ=255)
              cACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                        <  1:2858984/70‑1 (MQ=255)
                 cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                    >  1:788805/1‑71 (MQ=255)
                 cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                    >  1:2906151/1‑71 (MQ=255)
                 cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                    >  1:189977/1‑71 (MQ=255)
                 cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                    >  1:5163860/1‑71 (MQ=255)
                  gCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                                                   >  1:5046768/1‑71 (MQ=255)
                  gCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGAt                                                   <  1:3122600/71‑1 (MQ=255)
                   cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                  <  1:3000213/71‑1 (MQ=255)
                   cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                  <  1:4133638/71‑1 (MQ=255)
                   cTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATc                                                  <  1:4118388/71‑1 (MQ=255)
                      cGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCa                                                <  1:290012/70‑1 (MQ=255)
                      cGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCa                                                <  1:1328439/70‑1 (MQ=255)
                       gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:2193112/70‑1 (MQ=255)
                       gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:2727331/70‑1 (MQ=255)
                       gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:3180425/70‑1 (MQ=255)
                       gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:5063043/70‑1 (MQ=255)
                        ccTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGAt                                             <  1:2532187/71‑1 (MQ=255)
                          ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGa                                                                       >  1:5271891/1‑43 (MQ=255)
                          ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                                                      <  1:1672559/44‑1 (MQ=255)
                          ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                                                    >  1:3055530/1‑46 (MQ=255)
                          ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTAGATTATGTTCCAGTTTTTCGATc                                                  <  1:796523/64‑1 (MQ=255)
                           tctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATccc                                          >  1:376745/1‑71 (MQ=255)
                            ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcccc                                         >  1:2551854/1‑71 (MQ=255)
                               ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                        <  1:323718/69‑1 (MQ=255)
                                  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATg                                                                  <  1:3716480/40‑1 (MQ=255)
                                  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATg                                                                  <  1:5297305/40‑1 (MQ=255)
                                  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                   <  1:3582069/71‑1 (MQ=255)
                                  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                   <  1:139382/71‑1 (MQ=255)
                                  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                   >  1:1005824/1‑71 (MQ=255)
                                  cTTTGCGCAGCTCTTCAGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                   >  1:5533592/1‑71 (MQ=255)
                                   tttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGAt                                             >  1:2850554/1‑60 (MQ=255)
                                    tggcgcAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGtt                                 >  1:3669571/3‑71 (MQ=255)
                                     tGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGtt                                 >  1:1342141/1‑70 (MQ=255)
                                      gcgcAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCa                                                <  1:2989331/54‑1 (MQ=255)
                                       cgcAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGa                               >  1:2236303/1‑70 (MQ=255)
                                       cgcAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGa                               <  1:4711421/70‑1 (MQ=255)
                                        tctgCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGa                               <  1:1433269/66‑1 (MQ=255)
                                        gcAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGAtt                             <  1:458306/71‑2 (MQ=255)
                                         cAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCt                            <  1:2214529/71‑1 (MQ=255)
                                            ctctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCt                            >  1:748674/1‑68 (MQ=255)
                                             tctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGc                        >  1:2771865/1‑71 (MQ=255)
                                              ctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                        <  1:4661426/54‑1 (MQ=255)
                                              ctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCt                       >  1:3193210/1‑71 (MQ=255)
                                              ctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCt                       >  1:1180815/1‑71 (MQ=255)
                                                 ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                                                   >  1:3977174/1‑40 (MQ=255)
                                                 ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                                                   >  1:5459567/1‑40 (MQ=255)
                                                 ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATACCCTGACGGTTGATCTTCGCt                       <  1:5277373/68‑1 (MQ=255)
                                                  cGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                   <  1:190477/71‑1 (MQ=255)
                                                      atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               >  1:1319784/1‑39 (MQ=255)
                                                      atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                     >  1:263020/1‑49 (MQ=255)
                                                      atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                     >  1:1889572/1‑49 (MQ=255)
                                                       tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                 <  1:4241631/68‑1 (MQ=255)
                                                       tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa              <  1:1601331/71‑1 (MQ=255)
                                                         tatTTGCGTGGATTATGTTCCAGTTTTTCGATCCAGATcccc                                         >  1:182441/1‑42 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCa                                                >  1:3672/1‑35 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCa                                                >  1:2948184/1‑35 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:713891/36‑1 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:1151885/36‑1 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                               <  1:3442904/36‑1 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                 >  1:685870/1‑66 (MQ=255)
                                                         tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                 >  1:3493630/1‑66 (MQ=255)
                                                          atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                        >  1:691836/1‑42 (MQ=255)
                                                          atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                        >  1:2851571/1‑42 (MQ=255)
                                                           tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGt          <  1:562367/71‑1 (MQ=255)
                                                            ttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGt          >  1:5620087/1‑70 (MQ=255)
                                                            tggCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGt          >  1:3042698/3‑70 (MQ=255)
                                                             tGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTc                          >  1:4645775/1‑53 (MQ=255)
                                                             ggCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTAGCTTTAATGTTACGGt          >  1:5529572/2‑69 (MQ=255)
                                                              gCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTcag       <  1:493611/71‑1 (MQ=255)
                                                              gCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag       >  1:1544333/1‑71 (MQ=255)
                                                               ccTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTcag       <  1:2585167/70‑1 (MQ=255)
                                                                cTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTcagcg     <  1:765686/71‑2 (MQ=255)
                                                                 tGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGtt                                 <  1:2352472/42‑1 (MQ=16)
                                                                  ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag    >  1:5459541/1‑70 (MQ=255)
                                                                  ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca     >  1:3118813/1‑69 (MQ=255)
                                                                  ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca     >  1:2699696/1‑69 (MQ=255)
                                                                  ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTAGGTcagcag    >  1:2488653/1‑69 (MQ=255)
                                                                   gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                   >  1:2269201/1‑54 (MQ=255)
                                                                   gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag  >  1:3866169/1‑71 (MQ=255)
                                                                   gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag  >  1:4189154/1‑71 (MQ=255)
                                                                   gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag  >  1:4853656/1‑71 (MQ=255)
                                                                   gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag  <  1:4837800/71‑1 (MQ=255)
                                                                   gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTAGCTTTGATATTACGGTCagcaga   <  1:489632/70‑1 (MQ=255)
                                                                    aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                   >  1:5114168/1‑53 (MQ=255)
                                                                    aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                   >  1:2264516/1‑53 (MQ=255)
                                                                    |                                                                     
TTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAG  >  minE/1014978‑1015115

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: