Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,852 T→G 14.5% N406T (AAC→ACC) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8520TG14.5% 242.3 / 40.7 117N406T (AAC→ACC) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (51/49);  new base G (5/12);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 127 total aligned reads displayed.

TCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAA  >  minE/1014783‑1014922
                                                                     |                                                                      
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:3490139/1‑71 (MQ=255)
         ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                              >  1:1345685/1‑71 (MQ=255)
         ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                              >  1:2099960/1‑71 (MQ=255)
           ggCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGc                                                             <  1:5612886/70‑1 (MQ=255)
             ccGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGctct                                                          <  1:3958427/71‑1 (MQ=255)
              cGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                         <  1:5068988/71‑1 (MQ=255)
               ggACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCttt                                                        <  1:4819662/71‑1 (MQ=255)
                  cccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                     <  1:2653573/71‑1 (MQ=255)
                  cccGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                     >  1:4683493/1‑71 (MQ=255)
                     gcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                  <  1:4713331/71‑1 (MQ=255)
                        tagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                       <  1:855344/63‑1 (MQ=255)
                          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)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                        >  1:2009132/1‑41 (MQ=255)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt                                                                        >  1:2718549/1‑41 (MQ=255)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                             >  1:2767455/1‑68 (MQ=255)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                          >  1:4459167/1‑71 (MQ=255)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                          <  1:3039072/71‑1 (MQ=255)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                          >  1:5363902/1‑71 (MQ=255)
                             aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                          <  1:2528025/71‑1 (MQ=255)
                              aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                <  1:1891405/64‑1 (MQ=255)
                              aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                <  1:571955/64‑1 (MQ=255)
                               tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                         <  1:5657746/54‑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)
                                 ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                      <  1:5541178/71‑1 (MQ=255)
                                   atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                  <  1:1704753/57‑1 (MQ=255)
                                   atgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat                                    <  1:401998/71‑1 (MQ=255)
                                    tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                                 <  1:845235/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)
                                     gatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCa                                  >  1:304630/1‑71 (MQ=255)
                                      atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                         <  1:3620500/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)
                                        cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                       >  1:1393894/1‑47 (MQ=255)
                                        cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                       >  1:1768194/1‑47 (MQ=255)
                                        cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                       >  1:1976161/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)
                                                tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACcgacg                         >  1:2855539/1‑69 (MQ=255)
                                                tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                       <  1:1044178/71‑1 (MQ=255)
                                                tACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                       <  1:2416040/71‑1 (MQ=255)
                                                tACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                       <  1:4260998/71‑1 (MQ=255)
                                                  cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                  <  1:2640018/42‑1 (MQ=255)
                                                  cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                  <  1:4970889/42‑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:4221716/1‑70 (MQ=255)
                                                  cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATa                      <  1:2024512/70‑1 (MQ=255)
                                                    gCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGAGTGAGCAat                                      >  1:5339659/1‑52 (MQ=255)
                                                     cGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCg                  >  1:4918468/1‑71 (MQ=255)
                                                      ggCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGt                                                    <  1:3230830/36‑1 (MQ=255)
                                                        cAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                              >  1:5531401/1‑56 (MQ=255)
                                                          gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                       <  1:2042683/61‑1 (MQ=255)
                                                          gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg             >  1:1479363/1‑71 (MQ=255)
                                                          gCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg             >  1:359039/1‑71 (MQ=255)
                                                          gCCCTTCACCGTAATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                       <  1:5290922/61‑1 (MQ=255)
                                                           cccTTCACCGTTATCCATCCGCTCTTTCAGTTTgg                                                <  1:1365618/35‑1 (MQ=255)
                                                             cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc               >  1:3192557/1‑66 (MQ=255)
                                                             cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt           >  1:1772105/1‑70 (MQ=255)
                                                             cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa          <  1:348849/71‑1 (MQ=255)
                                                             cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa          <  1:5256898/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:511593/1‑51 (MQ=255)
                                                                cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAc                                 >  1:2372162/1‑45 (MQ=255)
                                                                cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg        <  1:321464/70‑1 (MQ=255)
                                                                cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGc       >  1:4912231/1‑71 (MQ=255)
                                                                 aCCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCg      <  1:2110683/71‑1 (MQ=255)
                                                                  ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                      >  1:781696/1‑38 (MQ=255)
                                                                  ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                      >  1:883362/1‑38 (MQ=255)
                                                                  ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                      >  1:4881699/1‑38 (MQ=255)
                                                                   cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                <  1:3509224/43‑1 (MQ=255)
                                                                    gTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGAt    >  1:2827710/1‑70 (MQ=255)
                                                                     ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                              >  1:4098501/1‑43 (MQ=255)
                                                                     ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                              >  1:3658776/1‑43 (MQ=255)
                                                                     ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATaa  <  1:2487399/71‑1 (MQ=255)
                                                                     |                                                                      
TCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAA  >  minE/1014783‑1014922

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