Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,938 A→G 5.5% Y377Y (TAT→TAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,9380AG 5.5% 318.2 / 6.1 108Y377Y (TAT→TACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (47/55);  new base G (2/4);  total (49/59)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 118 total aligned reads displayed.

TCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCT  >  minE/1014868‑1015007
                                                                      |                                                                     
tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       >  1:1602650/1‑71 (MQ=255)
 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       >  1:1814396/1‑70 (MQ=255)
 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       <  1:3911465/70‑1 (MQ=255)
 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       <  1:4176942/70‑1 (MQ=255)
 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       <  1:5000338/70‑1 (MQ=255)
   gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc                                                                    <  1:1345450/71‑1 (MQ=255)
   gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc                                                                    <  1:4520888/71‑1 (MQ=255)
      tGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGc                                                                 >  1:528083/1‑71 (MQ=255)
      tGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGc                                                                 >  1:5207530/1‑71 (MQ=255)
        gCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGa                                                               >  1:2852516/1‑71 (MQ=255)
        gCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGa                                                               >  1:5280783/1‑71 (MQ=255)
          gTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGAc                                                              >  1:5142449/1‑70 (MQ=255)
           tGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       <  1:2126123/60‑1 (MQ=255)
             aGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                          >  1:1396740/1‑71 (MQ=255)
             aGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                          >  1:102225/1‑71 (MQ=255)
                 atatCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGAc                                                              <  1:1835641/63‑1 (MQ=255)
                  tatCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGc                                                      >  1:5199266/1‑70 (MQ=255)
                  tatCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCg                                                     <  1:2578893/71‑1 (MQ=255)
                    tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAAc                                                           >  1:4159673/1‑63 (MQ=255)
                    tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAAc                                                           >  1:1075172/1‑63 (MQ=255)
                    tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:3168626/1‑71 (MQ=255)
                    tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:5663141/1‑71 (MQ=255)
                       cgaccgacGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                <  1:3468341/71‑1 (MQ=255)
                        gaccgacGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               <  1:3212204/71‑1 (MQ=255)
                        gaccgacGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               <  1:2712482/71‑1 (MQ=255)
                          ccgacGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       <  1:4707559/45‑1 (MQ=255)
                           cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               >  1:3748881/1‑68 (MQ=255)
                           cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                            <  1:4903658/71‑1 (MQ=255)
                           cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                            <  1:703900/71‑1 (MQ=255)
                           cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                            <  1:2169117/71‑1 (MQ=255)
                           ccacgaTAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATctct                                            >  1:5406507/3‑71 (MQ=255)
                            gacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCtt                                           <  1:5125769/71‑1 (MQ=255)
                            cacgaTAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCtt                                           <  1:4133051/70‑1 (MQ=255)
                              cgaTAATGGTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTTCGAAAGCTGGGCGAGGATCTCTTTc                                         <  1:1733960/71‑1 (MQ=255)
                              cgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCttt                                          >  1:1726081/1‑70 (MQ=255)
                              cgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCttt                                          >  1:4447970/1‑70 (MQ=255)
                               gaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGAc                                                              <  1:1081887/49‑1 (MQ=255)
                               gaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGAc                                                              <  1:752029/49‑1 (MQ=255)
                               gaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGAc                                                              <  1:5636689/49‑1 (MQ=255)
                                aTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACt                                                                   <  1:268706/43‑1 (MQ=255)
                                aTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCAt                                       <  1:2848773/71‑1 (MQ=255)
                                 tAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAAc                                                           >  1:1669128/1‑50 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       >  1:149437/1‑37 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       >  1:3903844/1‑37 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa                                                                       >  1:3427710/1‑37 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGAt                                                                      >  1:2736351/1‑38 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTg                                                                  >  1:2832831/1‑42 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGc                                                                 >  1:5418795/1‑43 (MQ=255)
                                  aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCg                                     >  1:419350/1‑71 (MQ=255)
                                   aTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCg                                     <  1:2195968/70‑1 (MQ=255)
                                    tggTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCttt                                          <  1:4478725/62‑1 (MQ=255)
                                    tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACaa                                                            <  1:3826976/46‑1 (MQ=255)
                                    tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGa                                                    >  1:50036/1‑54 (MQ=255)
                                    tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:2490004/1‑55 (MQ=255)
                                    tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               <  1:1818387/59‑1 (MQ=255)
                                    tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGa                                   >  1:4508695/1‑71 (MQ=255)
                                      gTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                          >  1:1251165/1‑46 (MQ=255)
                                       tGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCg                                 >  1:4963361/1‑70 (MQ=255)
                                         ccGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGtt                              <  1:3604627/71‑1 (MQ=255)
                                         ccGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGtt                              <  1:876619/71‑1 (MQ=255)
                                         ccGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGtt                              <  1:3337592/71‑1 (MQ=255)
                                          cGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGttt                             <  1:30886/71‑2 (MQ=255)
                                          cGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTa                             <  1:3792735/71‑1 (MQ=255)
                                          cGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTa                             <  1:5377619/71‑1 (MQ=255)
                                          cGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTa                             <  1:1923949/71‑1 (MQ=255)
                                           gTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTaa                            <  1:1318515/71‑1 (MQ=255)
                                           gTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTaa                            <  1:5584027/71‑1 (MQ=255)
                                            ttAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:1960752/1‑50 (MQ=255)
                                              aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:589074/1‑40 (MQ=255)
                                              aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:2123600/1‑40 (MQ=255)
                                              aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:2066077/1‑40 (MQ=255)
                                              aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGt                         >  1:4422320/1‑71 (MQ=255)
                                               aGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                        <  1:1479785/71‑1 (MQ=255)
                                               aGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                        <  1:5274508/71‑1 (MQ=255)
                                                gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:3349422/1‑46 (MQ=255)
                                                gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:1115304/1‑46 (MQ=255)
                                                gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCa                       <  1:3796929/71‑1 (MQ=255)
                                                gCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCa                       <  1:3093759/71‑1 (MQ=255)
                                                 cGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCaa                      <  1:163267/71‑1 (MQ=255)
                                                   aTAAGCGTGTAGAAACGGGATACTGCGACACCTGTGCGAGGATCTCTTTCAt                                       <  1:2803054/52‑1 (MQ=255)
                                                   aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcg                    <  1:4589979/71‑1 (MQ=255)
                                                   aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcg                    <  1:2166107/71‑1 (MQ=255)
                                                   aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTACATCGGACGGTTAAGGTCAAcg                    <  1:5308426/71‑1 (MQ=255)
                                                    tAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcgc                   <  1:2421370/71‑1 (MQ=255)
                                                    tAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAc                     >  1:2453114/1‑69 (MQ=255)
                                                     aaGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcgcg                  >  1:5622352/1‑71 (MQ=255)
                                                      aGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACgcgc                 <  1:700668/71‑1 (MQ=255)
                                                          tgtAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCg             <  1:5153970/71‑1 (MQ=255)
                                                           gtAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCAt                                       <  1:5254472/44‑1 (MQ=255)
                                                            tAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCt           <  1:734071/71‑1 (MQ=255)
                                                            tAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCt           <  1:5422825/71‑1 (MQ=255)
                                                              gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTc         <  1:5587656/71‑1 (MQ=255)
                                                              gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTc         <  1:3522649/71‑1 (MQ=255)
                                                               aaaCGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCg        <  1:1062587/71‑1 (MQ=255)
                                                                aaCGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCg        <  1:3743063/70‑1 (MQ=255)
                                                                aaCGGGATACTGCGACAACTGTGCGAGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGcc      <  1:5172977/71‑1 (MQ=255)
                                                                 aCGGGATACTGTGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGcc      <  1:1118155/71‑1 (MQ=255)
                                                                 aCGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCa                       >  1:3647560/1‑54 (MQ=255)
                                                                 aCGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGcc      >  1:5463478/1‑71 (MQ=255)
                                                                     gATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTtct  >  1:3546580/1‑71 (MQ=255)
                                                                      |                                                                     
TCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCT  >  minE/1014868‑1015007

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: