Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,854 A→G 14.1% D405D (GAT→GAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8540AG14.1% 289.0 / 33.9 120D405D (GAT→GACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (52/51);  new base G (5/12);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 123 total aligned reads displayed.

TTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAG  >  minE/1014790‑1014923
                                                                |                                                                     
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:1742759/1‑71 (MQ=255)
  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg                                                               >  1:2099960/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:855344/63‑1 (MQ=255)
                 tagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                        <  1:5550914/63‑1 (MQ=255)
                   gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                              <  1:159008/71‑1 (MQ=255)
                   gtaGATCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGt                                              >  1:591435/1‑71 (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)
                      aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                              >  1:2767455/1‑68 (MQ=255)
                      aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                           <  1:2528025/71‑1 (MQ=255)
                      aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                           >  1:5363902/1‑71 (MQ=255)
                      aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                           >  1:3063084/1‑71 (MQ=255)
                      aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                           >  1:4459167/1‑71 (MQ=255)
                      aaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAg                                           <  1:3039072/71‑1 (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)
                          ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                      >  1:3956383/1‑56 (MQ=255)
                          ggATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAg                                                      >  1:4045622/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: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)
                             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)
                               atCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCtt                                                          <  1:1608258/47‑1 (MQ=255)
                               atCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAatat                                     <  1:1186710/68‑1 (MQ=255)
                                tCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                   >  1:1096361/1‑53 (MQ=255)
                                tCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                                 >  1:3013531/1‑71 (MQ=255)
                                 cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                        >  1:1976161/1‑47 (MQ=255)
                                 cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                        >  1:1393894/1‑47 (MQ=255)
                                 cTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTc                                                        >  1:1768194/1‑47 (MQ=255)
                                    tGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                 <  1:1014956/67‑1 (MQ=255)
                                     gATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcggccc                             <  1:3442358/70‑2 (MQ=255)
                                      aTGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcggccca                            <  1:2628690/70‑3 (MQ=255)
                                         tACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGAt                        <  1:1044178/71‑1 (MQ=255)
                                         tACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTCGGCCCACGAt                        <  1:4260998/71‑1 (MQ=255)
                                           cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                   <  1:3486280/42‑1 (MQ=255)
                                           cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                   <  1:4970889/42‑1 (MQ=255)
                                           cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                   <  1:2640018/42‑1 (MQ=255)
                                           cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                        <  1:3857005/53‑1 (MQ=255)
                                           cTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                        <  1:1117566/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)
                                                   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)
                                                      cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAGTCGTGc                >  1:1366287/1‑66 (MQ=255)
                                                      cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc                >  1:3192557/1‑66 (MQ=255)
                                                      cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt            >  1:5052641/1‑70 (MQ=255)
                                                      cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt            >  1:1772105/1‑70 (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)
                                                         cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAc                                  >  1:2737524/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:883362/1‑38 (MQ=255)
                                                           ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                       >  1:4881699/1‑38 (MQ=255)
                                                           ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                 <  1:194354/44‑1 (MQ=255)
                                                            cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcg                                 <  1:3509224/43‑1 (MQ=255)
                                                             gTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGAt     >  1:2827710/1‑70 (MQ=255)
                                                              ttATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                               >  1:4105199/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)
                                                                |                                                                     
TTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAG  >  minE/1014790‑1014923

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: