Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,843 C→T 18.9% G409E (GGG→GAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8430CT18.9% 230.0 / 41.5 107G409E (GGG→GAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (46/40);  new base T (8/12);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 109 total aligned reads displayed.

ACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCG  >  minE/1014779‑1014911
                                                                |                                                                    
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)
      cGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGccc                                                                     <  1:2412323/60‑1 (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:3490139/1‑71 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                   >  1:1345685/1‑71 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                   >  1:2099960/1‑71 (MQ=255)
             ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                   >  1:1742759/1‑71 (MQ=255)
               ggCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGc                                                  <  1:5612886/70‑1 (MQ=255)
                 ccGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                               <  1:4025183/71‑1 (MQ=255)
                 ccGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                               <  1:3958427/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)
                      cccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                          >  1:4683493/1‑71 (MQ=255)
                         gcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                       <  1:4713331/71‑1 (MQ=255)
                            tagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                            <  1:5550914/63‑1 (MQ=255)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCggc                                                            >  1:1577843/1‑43 (MQ=12)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCggc                                                            >  1:1971180/1‑43 (MQ=12)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCggc                                                            >  1:2740048/1‑43 (MQ=12)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                  >  1:591435/1‑71 (MQ=255)
                              gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                  <  1:159008/71‑1 (MQ=255)
                              gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                               <  1:2701980/58‑1 (MQ=255)
                              gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                               <  1:669181/58‑1 (MQ=255)
                              gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                          >  1:4561844/1‑63 (MQ=255)
                               taAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                  <  1:5614760/70‑1 (MQ=255)
                                aaaTCGGATGATCTTTGATGTACTGCGGCAGCCCt                                                                    >  1:3967715/1‑35 (MQ=255)
                                aaaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc                                                                <  1:1106746/39‑1 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                             >  1:4673677/1‑41 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                             >  1:2718549/1‑41 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                  >  1:2767455/1‑68 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                               >  1:5363902/1‑71 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                               >  1:4459167/1‑71 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                               >  1:3063084/1‑71 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                               <  1:3039072/71‑1 (MQ=255)
                                 aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                               <  1:2528025/71‑1 (MQ=255)
                                  aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                     <  1:1891405/64‑1 (MQ=255)
                                  aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                     <  1:571955/64‑1 (MQ=255)
                                  aTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCg                                                              >  1:2889820/1‑39 (MQ=255)
                                  aTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCg                                                              >  1:2794645/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)
                                     gggtgatCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAat                           <  1:333256/68‑1 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                          >  1:4045622/1‑56 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                          >  1:3956383/1‑56 (MQ=255)
                                     ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                            <  1:2865223/70‑1 (MQ=255)
                                       atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                       <  1:1704753/57‑1 (MQ=255)
                                       atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat                         <  1:584630/71‑1 (MQ=255)
                                       atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat                         <  1:401998/71‑1 (MQ=255)
                                        tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                      <  1:845235/57‑1 (MQ=255)
                                        tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                      <  1:793295/57‑1 (MQ=255)
                                        tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATc                        <  1:3127483/71‑1 (MQ=255)
                                        ggatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                          >  1:5177737/2‑69 (MQ=255)
                                         gatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATc                        >  1:2108383/1‑70 (MQ=255)
                                          atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                              <  1:1608258/47‑1 (MQ=255)
                                          atCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAatat                         <  1:1186710/68‑1 (MQ=255)
                                           tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                       >  1:1742229/1‑53 (MQ=255)
                                           tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                       >  1:1096361/1‑53 (MQ=255)
                                           tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                           <  1:5089026/65‑1 (MQ=255)
                                           tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                     >  1:4376413/1‑71 (MQ=255)
                                           tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                     >  1:3013531/1‑71 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                            >  1:1976161/1‑47 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                            >  1:1768194/1‑47 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                            >  1:1393894/1‑47 (MQ=255)
                                            cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                            >  1:3973381/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)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                            <  1:3857005/53‑1 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcga                    >  1:1522583/1‑61 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATaa          >  1:4373141/1‑71 (MQ=255)
                                                      cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATa           <  1:2024512/70‑1 (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)
                                                              gCCCTTCACCGTAATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt            <  1:5290922/61‑1 (MQ=255)
                                                               cccTTCACCGTTATCCATCCGCTCTTTCAGTTTgg                                     <  1:1365618/35‑1 (MQ=255)
                                                                |                                                                    
ACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCG  >  minE/1014779‑1014911

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: