Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,076 C→T 36.6% Q331Q (CAG→CAA) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0760CT36.6% 97.1 / 93.3 101Q331Q (CAG→CAA) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base T (14/23);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

Only 100 of 109 total aligned reads displayed.

CTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATG  >  minE/1015006‑1015146
                                                                      |                                                                      
ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcccc                                                                        >  1:2551854/1‑71 (MQ=255)
   ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                                                       <  1:323718/69‑1 (MQ=255)
      cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                  >  1:1005824/1‑71 (MQ=255)
      cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                  <  1:3582069/71‑1 (MQ=255)
      cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                  <  1:139382/71‑1 (MQ=255)
      cTTTGCGCAGCTCTTCAGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                  >  1:5533592/1‑71 (MQ=255)
        tggcgcAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGtt                                                                >  1:3669571/3‑71 (MQ=255)
         tGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGtt                                                                >  1:1342141/1‑70 (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:1180815/1‑71 (MQ=255)
                  ctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCt                                                      >  1:3193210/1‑71 (MQ=255)
                     ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATACCCTGACGGTTGATCTTCGCt                                                      <  1:5277373/68‑1 (MQ=255)
                      cGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                  <  1:190477/71‑1 (MQ=255)
                          atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                                                    >  1:263020/1‑49 (MQ=255)
                          atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                                                    >  1:1889572/1‑49 (MQ=255)
                           tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                <  1:4241631/68‑1 (MQ=255)
                             tatTTGCGTGGATTATGTTCCAGTTTTTCGATCCAGATcccc                                                                        >  1:182441/1‑42 (MQ=255)
                             tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                >  1:3493630/1‑66 (MQ=255)
                             tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                >  1:685870/1‑66 (MQ=255)
                              atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                                                       >  1:691836/1‑42 (MQ=255)
                              atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCt                                                                       >  1:2851571/1‑42 (MQ=255)
                                ttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGt                                         >  1:5620087/1‑70 (MQ=255)
                                tggCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGt                                         >  1:3042698/3‑70 (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:2699696/1‑69 (MQ=255)
                                      ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                    >  1:3118813/1‑69 (MQ=255)
                                      ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTAGGTcagcag                                   >  1:2488653/1‑69 (MQ=255)
                                       gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                  >  1:2269201/1‑54 (MQ=255)
                                       gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                 >  1:4853656/1‑71 (MQ=255)
                                       gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                 <  1:4837800/71‑1 (MQ=255)
                                       gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                 >  1:3866169/1‑71 (MQ=255)
                                       gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                 >  1:4189154/1‑71 (MQ=255)
                                       gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTAGCTTTGATATTACGGTCagcaga                                  <  1:489632/70‑1 (MQ=255)
                                        aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                  >  1:2264516/1‑53 (MQ=255)
                                         ttATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagc                                <  1:2001343/70‑1 (MQ=255)
                                         ttATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                 >  1:943428/1‑69 (MQ=255)
                                          tgtGTTCCAGTTTTTCGATGCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagc                                <  1:4692495/67‑1 (MQ=255)
                                           gtGTTCCAGTTTTTCGATCCATATACCTTCGCGGTTGATTTTCGCTTTAATGTTACCGTCagcggagcag                              <  1:3984026/69‑1 (MQ=255)
                                           gtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                              <  1:812238/69‑1 (MQ=255)
                                            tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                            <  1:1212523/71‑1 (MQ=255)
                                            tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                            <  1:3690239/71‑1 (MQ=255)
                                            tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                            <  1:4598355/71‑1 (MQ=255)
                                            tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                            <  1:5354398/71‑1 (MQ=255)
                                            tGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCtt                                                     <  1:2805948/46‑1 (MQ=25)
                                            tGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCtt                                                     <  1:4131354/46‑1 (MQ=25)
                                            tGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                            <  1:19175/71‑1 (MQ=255)
                                             gTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                            <  1:1146097/70‑1 (MQ=255)
                                              ttCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                          >  1:1604107/1‑71 (MQ=255)
                                              ttCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                          >  1:3603061/1‑71 (MQ=255)
                                              ttCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                          <  1:1525592/71‑1 (MQ=255)
                                              ttCCAGTTTTTCGATCCAGATACCTTAGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                          <  1:4348262/71‑1 (MQ=255)
                                               tCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                         >  1:3339198/1‑71 (MQ=255)
                                                ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                         <  1:5284528/70‑1 (MQ=255)
                                                ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACg                        >  1:55464/1‑71 (MQ=255)
                                                ccAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTa                                             >  1:5324530/1‑50 (MQ=255)
                                                ccAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACg                        >  1:1805971/1‑71 (MQ=255)
                                                 cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                              <  1:425225/64‑1 (MQ=255)
                                                 cAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGGCAGCGGAGCAGGAGACGc                       <  1:4591599/71‑1 (MQ=255)
                                                  aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagca                               >  1:2883009/1‑62 (MQ=255)
                                                  aGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCttt                                                    <  1:465666/41‑1 (MQ=19)
                                                  aGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCttt                                                    <  1:5193442/41‑1 (MQ=19)
                                                   gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                              <  1:988283/62‑1 (MQ=255)
                                                   gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                              <  1:3605047/62‑1 (MQ=255)
                                                   gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                              <  1:1152191/62‑1 (MQ=255)
                                                    tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                   <  1:1178992/56‑1 (MQ=255)
                                                    tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                      <  1:1659209/53‑1 (MQ=255)
                                                    tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                      <  1:4554022/53‑1 (MQ=255)
                                                    tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                      <  1:939032/53‑1 (MQ=255)
                                                     ttttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAatg                                                >  1:610987/1‑41 (MQ=16)
                                                     ttttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAatg                                                >  1:3770352/1‑41 (MQ=16)
                                                      tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                  >  1:4828660/1‑39 (MQ=255)
                                                      tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagca                               >  1:195876/1‑58 (MQ=255)
                                                      tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                          >  1:1768927/1‑63 (MQ=255)
                                                      tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCAtg                  >  1:5106137/1‑70 (MQ=255)
                                                         cGATCCAGATACCTTCGGGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCTGGAGACGCCCATgcc                <  1:959164/70‑1 (MQ=255)
                                                         cGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgcc                >  1:631225/1‑70 (MQ=255)
                                                          gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccga              <  1:2196846/71‑1 (MQ=255)
                                                             ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                    >  1:1112260/1‑46 (MQ=255)
                                                             ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                    >  1:1884346/1‑46 (MQ=255)
                                                             ccAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccgaccg           >  1:2889007/1‑71 (MQ=255)
                                                                gATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                    <  1:4953286/43‑1 (MQ=255)
                                                                   ccccTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg      >  1:2236504/1‑70 (MQ=255)
                                                                   ccccTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCACATACCGACCGGGCAg      >  1:2294138/1‑70 (MQ=255)
                                                                    cccTCACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGa    <  1:3037613/71‑1 (MQ=255)
                                                                      cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGCCCGGGCAGGAt   >  1:5521129/1‑70 (MQ=255)
                                                                      cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg  >  1:3316979/1‑71 (MQ=255)
                                                                      cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg  >  1:4567470/1‑71 (MQ=255)
                                                                      |                                                                      
CTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATG  >  minE/1015006‑1015146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: