Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,037 T→C 26.1% K344K (AAA→AAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0370TC26.1% 98.7 / 37.8 65K344K (AAA→AAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base C (7/10);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGT  >  minE/1014968‑1015107
                                                                     |                                                                      
cATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatatat                                                                        <  1:5345646/70‑1 (MQ=255)
      gcGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTg                                                                 <  1:1535543/70‑1 (MQ=255)
        ggTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTg                                                                 <  1:3973303/69‑1 (MQ=255)
          ttAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTgggt                                                              >  1:3086546/1‑68 (MQ=255)
          ttAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGc                                                                    >  1:994112/1‑64 (MQ=255)
          ttAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGc                                                                    >  1:4762380/1‑64 (MQ=255)
          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:3939654/69‑1 (MQ=255)
                       tcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                                                 <  1:2587755/69‑1 (MQ=255)
                       gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                <  1:1686862/71‑1 (MQ=255)
                        cACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                <  1:2858984/70‑1 (MQ=255)
                          ccGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATAttt                                                                      <  1:3619840/46‑1 (MQ=255)
                           cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                            >  1:2906151/1‑71 (MQ=255)
                           cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                            >  1:5163860/1‑71 (MQ=255)
                           cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                            >  1:189977/1‑71 (MQ=255)
                           cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                            >  1:788805/1‑71 (MQ=255)
                            gCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                                           >  1:5046768/1‑71 (MQ=255)
                            gCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGAt                                           <  1:3122600/71‑1 (MQ=255)
                             cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                          <  1:4133638/71‑1 (MQ=255)
                             cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                          <  1:3000213/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:5063043/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)
                                  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:5297305/40‑1 (MQ=255)
                                            cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATg                                                          <  1:3716480/40‑1 (MQ=255)
                                            cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                           >  1:1005824/1‑71 (MQ=255)
                                            cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                           <  1:139382/71‑1 (MQ=255)
                                            cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                           <  1:3582069/71‑1 (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:1889572/1‑49 (MQ=255)
                                                                atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                             >  1:263020/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:1151885/36‑1 (MQ=255)
                                                                   tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAg                                       <  1:713891/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)
                                                                     |                                                                      
CATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGT  >  minE/1014968‑1015107

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: