Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,847 C→T 17.1% E408K (GAA→AAA)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8470CT17.1% 269.6 / 41.0 117E408K (GAA→AAA) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (51/46);  new base T (8/12);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 124 total aligned reads displayed.

ACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGC  >  minE/1014779‑1014917
                                                                    |                                                                      
aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa                                                                       >  1:5414016/1‑70 (MQ=255)
aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa                                                                       >  1:3814459/1‑70 (MQ=255)
aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc                                                                      >  1:101591/1‑71 (MQ=255)
 ccTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc                                                                     >  1:4527079/1‑71 (MQ=255)
    tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc                                                                     >  1:2011655/1‑68 (MQ=255)
    tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                   >  1:3026398/1‑70 (MQ=255)
           ttttGGCCGGACGCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatccatc                                                           >  1:4626157/1‑71 (MQ=255)
           ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                   <  1:1013865/63‑1 (MQ=255)
           ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatccatc                                                           >  1:3414559/1‑71 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatcca                                                             >  1:4590974/1‑67 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                         >  1:1742759/1‑71 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                         >  1:2099960/1‑71 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                         >  1:1345685/1‑71 (MQ=255)
                 ccGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                     <  1:3958427/71‑1 (MQ=255)
                 ccGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                     <  1:4025183/71‑1 (MQ=255)
                  cGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCtt                                                    <  1:386899/71‑1 (MQ=255)
                  cGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                    <  1:5068988/71‑1 (MQ=255)
                   ggACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                   <  1:4819662/71‑1 (MQ=255)
                      cccGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAg                                                >  1:488525/1‑71 (MQ=255)
                      cccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                >  1:4683493/1‑71 (MQ=255)
                            tagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                  <  1:5550914/63‑1 (MQ=255)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCggc                                                                  >  1:2740048/1‑43 (MQ=12)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCggc                                                                  >  1:1971180/1‑43 (MQ=12)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                        <  1:159008/71‑1 (MQ=255)
                              gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                     <  1:2701980/58‑1 (MQ=255)
                              gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                     <  1:669181/58‑1 (MQ=255)
                               taAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                        <  1:5614760/70‑1 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                   >  1:2009132/1‑41 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                   >  1:2718549/1‑41 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                   >  1:4673677/1‑41 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                     >  1:4459167/1‑71 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                     <  1:3039072/71‑1 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                     >  1:3063084/1‑71 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                     <  1:2528025/71‑1 (MQ=255)
                                  aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                           <  1:571955/64‑1 (MQ=255)
                                  aTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCg                                                                    >  1:2794645/1‑39 (MQ=255)
                                  aTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCg                                                                    >  1:2889820/1‑39 (MQ=255)
                                   tCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                   <  1:2257829/71‑1 (MQ=255)
                                   tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                    <  1:5657746/54‑1 (MQ=255)
                                   tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTgg                                           <  1:5341876/63‑1 (MQ=255)
                                    cGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCaa                                  >  1:76916/1‑71 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                >  1:3956383/1‑56 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                >  1:4045622/1‑56 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                  <  1:2865223/70‑1 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                 <  1:5541178/71‑1 (MQ=255)
                                       atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                             <  1:1704753/57‑1 (MQ=255)
                                       atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat                               <  1:584630/71‑1 (MQ=255)
                                        tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                            <  1:845235/57‑1 (MQ=255)
                                        tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                            <  1:793295/57‑1 (MQ=255)
                                         gatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATc                              >  1:2108383/1‑70 (MQ=255)
                                         gatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCa                             >  1:304630/1‑71 (MQ=255)
                                          atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                    <  1:3620500/47‑1 (MQ=255)
                                          atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                    <  1:1608258/47‑1 (MQ=255)
                                          atCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAatat                               <  1:1186710/68‑1 (MQ=255)
                                           tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                             >  1:1096361/1‑53 (MQ=255)
                                           tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                 <  1:5089026/65‑1 (MQ=255)
                                           tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                           >  1:4376413/1‑71 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                  >  1:3973381/1‑47 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                  >  1:1393894/1‑47 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                  >  1:1768194/1‑47 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                  >  1:1976161/1‑47 (MQ=255)
                                              ttGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgacc                        >  1:1243314/1‑71 (MQ=255)
                                               tGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                           <  1:1014956/67‑1 (MQ=255)
                                                gATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcggccc                       <  1:3442358/70‑2 (MQ=255)
                                                 aTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcggccca                      <  1:2628690/70‑3 (MQ=255)
                                                    tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACcgacg                    >  1:2855539/1‑69 (MQ=255)
                                                    tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                  <  1:1044178/71‑1 (MQ=255)
                                                    tACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                  <  1:2416040/71‑1 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                             <  1:2640018/42‑1 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                             <  1:4970889/42‑1 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                             <  1:3486280/42‑1 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                  <  1:1117566/53‑1 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcga                          >  1:1522583/1‑61 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATaa                >  1:4373141/1‑71 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATa                 >  1:4221716/1‑70 (MQ=255)
                                                        gCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGAGTGAGCAat                                 >  1:5339659/1‑52 (MQ=255)
                                                         cGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCg             >  1:4918468/1‑71 (MQ=255)
                                                          ggCAGCGCTTCACCGTTATCCATCCGCTCTTTCAGt                                               <  1:3237745/36‑1 (MQ=255)
                                                          ggCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGt                                               <  1:3230830/36‑1 (MQ=255)
                                                          ggCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCg             <  1:3571505/70‑1 (MQ=255)
                                                            cAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                         >  1:5531401/1‑56 (MQ=255)
                                                              gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                  <  1:2042683/61‑1 (MQ=255)
                                                              gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg        >  1:359039/1‑71 (MQ=255)
                                                              gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg        >  1:1479363/1‑71 (MQ=255)
                                                              gCCCTTCACCGTAATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                  <  1:5290922/61‑1 (MQ=255)
                                                                 cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAGTCGTGc          >  1:1366287/1‑66 (MQ=255)
                                                                 cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc          >  1:3192557/1‑66 (MQ=255)
                                                                 cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt      >  1:1772105/1‑70 (MQ=255)
                                                                 cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa     <  1:5256898/71‑1 (MQ=255)
                                                                 cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa     <  1:348849/71‑1 (MQ=255)
                                                                  ttCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgacc                        <  1:16583/51‑1 (MQ=255)
                                                                   tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa     <  1:5148124/69‑1 (MQ=255)
                                                                   tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg   >  1:4027284/1‑71 (MQ=255)
                                                                   tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCAACGATAATCGTGCCGTTAAg   >  1:174298/1‑71 (MQ=255)
                                                                    cACCGTTATCCATCCGCTCTTTCAGTTTGGCTTGAGCAATATCAcgaga                        >  1:3351006/1‑47 (MQ=39)
                                                                    cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgaccga                      >  1:2061168/1‑51 (MQ=255)
                                                                    cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgaccga                      >  1:511593/1‑51 (MQ=255)
                                                                    cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAc                            >  1:2737524/1‑45 (MQ=255)
                                                                    cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg   <  1:321464/70‑1 (MQ=255)
                                                                    cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGc  >  1:4912231/1‑71 (MQ=255)
                                                                    |                                                                      
ACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGC  >  minE/1014779‑1014917

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: