Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,872 T→C 15.5% K399K (AAA→AAG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8720TC15.5% 281.5 / 38.2 122K399K (AAA→AAGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (53/50);  new base C (4/15);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAC  >  minE/1014804‑1014941
                                                                    |                                                                     
gcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                                     <  1:4713331/71‑1 (MQ=255)
     gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                                                <  1:159008/71‑1 (MQ=255)
     gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                                                >  1:591435/1‑71 (MQ=255)
      taAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                <  1:5614760/70‑1 (MQ=255)
        aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                             <  1:2528025/71‑1 (MQ=255)
        aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                             >  1:5363902/1‑71 (MQ=255)
        aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                             <  1:3039072/71‑1 (MQ=255)
        aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                             >  1:4459167/1‑71 (MQ=255)
        aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                             >  1:3063084/1‑71 (MQ=255)
         aTCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTgg                                                                   <  1:1891405/64‑1 (MQ=255)
          tCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCa                                                           <  1:2257829/71‑1 (MQ=255)
          tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTgg                                                                   <  1:5341876/63‑1 (MQ=255)
           cGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCaa                                                          >  1:76916/1‑71 (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)
              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)
                  tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                                     >  1:1742229/1‑53 (MQ=255)
                  tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                                     >  1:1096361/1‑53 (MQ=255)
                  tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                                         <  1:5089026/65‑1 (MQ=255)
                  tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                                                   >  1:3013531/1‑71 (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)
                             cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                                     <  1:4970889/42‑1 (MQ=255)
                             cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                                     <  1:2640018/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)
                                 ggCAGCGCTTCACCGTTATCCATCCGCTCTTTCAGt                                                                       <  1:3237745/36‑1 (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:359039/1‑71 (MQ=255)
                                        cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAGTCGTGc                                  >  1:1366287/1‑66 (MQ=255)
                                        cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa                             <  1:5256898/71‑1 (MQ=255)
                                        cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa                             <  1:348849/71‑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)
                                           cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAc                                                    >  1:2372162/1‑45 (MQ=255)
                                           cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg                           <  1:321464/70‑1 (MQ=255)
                                           cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGc                          >  1:4912231/1‑71 (MQ=255)
                                            aCCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                                             >  1:5164418/1‑35 (MQ=255)
                                            aCCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCg                         <  1:2110683/71‑1 (MQ=255)
                                             ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                                         >  1:4881699/1‑38 (MQ=255)
                                             ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                                         >  1:781696/1‑38 (MQ=255)
                                             ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                                         >  1:883362/1‑38 (MQ=255)
                                              cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                                   <  1:3509224/43‑1 (MQ=255)
                                               gTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGAt                       >  1:2827710/1‑70 (MQ=255)
                                                ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                                 >  1:4105199/1‑43 (MQ=255)
                                                ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                                 >  1:3658776/1‑43 (MQ=255)
                                                ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATaa                     <  1:2487399/71‑1 (MQ=255)
                                                ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATaa                     <  1:4427120/71‑1 (MQ=255)
                                                 tATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAg                    >  1:4278265/1‑71 (MQ=255)
                                                   tcaatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg                  <  1:3695869/68‑1 (MQ=255)
                                                   tcaatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg                  <  1:792600/68‑1 (MQ=255)
                                                   tcaatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCGGTGAGCGACAAACg                  <  1:1676986/68‑1 (MQ=255)
                                                     cacccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgt                 <  1:5103627/67‑1 (MQ=255)
                                                     aatcaGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGATAATGGTGCCg                                <  1:430306/54‑1 (MQ=255)
                                                      atccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                                          <  1:2976672/44‑1 (MQ=255)
                                                       tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGcc                                 >  1:4689288/1‑52 (MQ=255)
                                                       tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGcc                                 >  1:1322693/1‑52 (MQ=255)
                                                       tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                >  1:3549813/1‑69 (MQ=255)
                                                         cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAg             >  1:1428937/1‑70 (MQ=255)
                                                         cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAg             >  1:933886/1‑70 (MQ=255)
                                                         cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGa            <  1:3308024/71‑1 (MQ=255)
                                                         cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGa            <  1:1470266/71‑1 (MQ=255)
                                                          gCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGaa           <  1:3616907/71‑1 (MQ=255)
                                                           ctctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGaaa          >  1:3602321/1‑71 (MQ=255)
                                                           ctctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGaaa          <  1:4951145/71‑1 (MQ=255)
                                                              tttCAGTTTGGCGTGAGCAATATCACGACCGACGAt                                          <  1:3673054/36‑1 (MQ=255)
                                                              tttCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAAc         <  1:4071831/69‑1 (MQ=255)
                                                               ttCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGt                                    <  1:3355889/41‑1 (MQ=255)
                                                                tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc                                  >  1:177860/1‑42 (MQ=255)
                                                                tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGa                        >  1:591413/1‑52 (MQ=255)
                                                                tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGa                        >  1:3061272/1‑52 (MQ=255)
                                                                tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGa            <  1:2489260/64‑1 (MQ=255)
                                                                tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa     >  1:1602650/1‑71 (MQ=255)
                                                                 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa     >  1:1814396/1‑70 (MQ=255)
                                                                 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa     <  1:5000338/70‑1 (MQ=255)
                                                                 cAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGa     <  1:4176942/70‑1 (MQ=255)
                                                                  aGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGaaa          >  1:244145/1‑64 (MQ=255)
                                                                   gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc  <  1:5365638/71‑1 (MQ=255)
                                                                   gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc  <  1:4520888/71‑1 (MQ=255)
                                                                   gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAc  <  1:1345450/71‑1 (MQ=255)
                                                                    |                                                                     
GCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATAC  >  minE/1014804‑1014941

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: