Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,941 C→T 5.6% Q376Q (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,014,9410CT 5.6% 325.4 / 5.4 106Q376Q (CAG→CAAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (48/52);  new base T (2/4);  total (50/56)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 110 total aligned reads displayed.

GTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCT  >  minE/1014871‑1015007
                                                                      |                                                                  
gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc                                                                    <  1:5365638/71‑1 (MQ=255)
gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc                                                                    <  1:4520888/71‑1 (MQ=255)
   tGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGc                                                                 >  1:5207530/1‑71 (MQ=255)
   tGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGc                                                                 >  1:528083/1‑71 (MQ=255)
     gCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGa                                                               >  1:2852516/1‑71 (MQ=255)
       gTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGAc                                                              >  1:5142449/1‑70 (MQ=255)
          aGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                          >  1:102225/1‑71 (MQ=255)
          aGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                          >  1:1396740/1‑71 (MQ=255)
          aGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACt                                                          >  1:501252/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)
                 tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAAc                                                           >  1:3412687/1‑63 (MQ=255)
                 tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:3168626/1‑71 (MQ=255)
                 tCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:3647118/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)
                        cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               >  1:3748881/1‑68 (MQ=255)
                        cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                            <  1:703900/71‑1 (MQ=255)
                        cgacgaTAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATctct                                            <  1:4903658/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:4535664/1‑50 (MQ=255)
                              tAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAAc                                                           >  1:1669128/1‑50 (MQ=255)
                               aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTg                                                                  >  1:2832831/1‑42 (MQ=255)
                               aaTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGc                                                                 >  1:5418795/1‑43 (MQ=255)
                                aTCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCg                                     <  1:2195968/70‑1 (MQ=255)
                                 tggTGCCGGTGAGCGACAAACGAGTGGATACCGGGTATTGCGAAAGCTGGGCGAGGATCTCttt                                          <  1:4478725/62‑1 (MQ=255)
                                 tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACaa                                                            <  1:3826976/46‑1 (MQ=255)
                                 tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:2716435/1‑55 (MQ=255)
                                 tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAg                                                   >  1:2490004/1‑55 (MQ=255)
                                 tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               <  1:1818387/59‑1 (MQ=255)
                                 tCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAtc                                               <  1:2305385/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: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:5584027/71‑1 (MQ=255)
                                        gTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTaa                            <  1:1318515/71‑1 (MQ=255)
                                         ttAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:1960752/1‑50 (MQ=255)
                                           aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:589074/1‑40 (MQ=255)
                                           aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:1777248/1‑40 (MQ=255)
                                           aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:2066077/1‑40 (MQ=255)
                                           aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACtgt                                                        >  1:2123600/1‑40 (MQ=255)
                                           aaGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGt                         >  1:4422320/1‑71 (MQ=255)
                                            aGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                        <  1:5274508/71‑1 (MQ=255)
                                            aGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                        <  1:1479785/71‑1 (MQ=255)
                                            aGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTc                        <  1:2106476/71‑1 (MQ=255)
                                             gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:1115304/1‑46 (MQ=255)
                                             gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:5489895/1‑46 (MQ=255)
                                             gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:3878176/1‑46 (MQ=255)
                                             gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:1593455/1‑46 (MQ=255)
                                             gCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGAt                                                >  1:3349422/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)
                                               gATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAc                     <  1:5200294/71‑1 (MQ=255)
                                                aTAAGCGTGTAGAAACGGGATACTGCGACACCTGTGCGAGGATCTCTTTCAt                                       <  1:2803054/52‑1 (MQ=255)
                                                aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcg                    <  1:2166107/71‑1 (MQ=255)
                                                aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcg                    <  1:4589979/71‑1 (MQ=255)
                                                aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGCTCAAcg                    <  1:6717/71‑1 (MQ=255)
                                                aTAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTACATCGGACGGTTAAGGTCAAcg                    <  1:5308426/71‑1 (MQ=255)
                                                 tAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAcgc                   <  1:2421370/71‑1 (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:5422825/71‑1 (MQ=255)
                                                         tAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCt           <  1:734071/71‑1 (MQ=255)
                                                           gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTc         <  1:3522649/71‑1 (MQ=255)
                                                           gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTc         <  1:4334335/71‑1 (MQ=255)
                                                           gAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTc         <  1:1180063/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)
                                                                    tACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACg                                 >  1:2859702/1‑38 (MQ=255)
                                                                    tACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACg                                 >  1:1055246/1‑38 (MQ=255)
                                                                      |                                                                  
GTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCT  >  minE/1014871‑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: