Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,698 A→G 19.3% S457S (AGT→AGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6980AG19.3% 167.5 / 33.7 83S457S (AGT→AGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (30/37);  new base G (9/7);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAG  >  minE/1014630‑1014765
                                                                    |                                                                   
aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCt                                                                   <  1:530231/71‑1 (MQ=255)
    cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                               >  1:4858543/1‑71 (MQ=255)
    cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                               <  1:4932291/71‑1 (MQ=255)
    cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                               <  1:5149872/71‑1 (MQ=255)
            tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcgc                                                       >  1:3348219/1‑71 (MQ=255)
                tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                                                    >  1:118167/1‑70 (MQ=255)
                tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                                                    <  1:2412482/70‑1 (MQ=255)
                  cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                   <  1:3500963/69‑1 (MQ=255)
                  cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCa                                                  <  1:1257073/70‑1 (MQ=255)
                      cATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                             >  1:1866316/1‑71 (MQ=255)
                       aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                                            >  1:785037/1‑71 (MQ=255)
                          cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATAc                                                                    >  1:3708461/1‑44 (MQ=255)
                          cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATAc                                                                    >  1:3604764/1‑44 (MQ=255)
                           gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                               <  1:5604358/48‑1 (MQ=255)
                           gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                               <  1:5680458/48‑1 (MQ=255)
                           gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCa                                                  >  1:72032/1‑61 (MQ=255)
                              aCCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACAt                                      >  1:2767763/1‑70 (MQ=255)
                                  gctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa                                 >  1:3249287/1‑71 (MQ=255)
                                  gctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGag                                 >  1:472056/1‑70 (MQ=255)
                                   ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTg                                              <  1:349343/57‑1 (MQ=255)
                                   ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGa                                  <  1:3328560/69‑1 (MQ=255)
                                   ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCAACATAAGa                                  <  1:2458362/69‑1 (MQ=255)
                                   ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTa                                     >  1:5578970/1‑66 (MQ=255)
                                     ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                   >  1:538187/1‑50 (MQ=255)
                                     gactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtc                               >  1:2740101/3‑70 (MQ=255)
                                       ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccat                            >  1:919753/1‑71 (MQ=255)
                                       ctgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccat                            <  1:635382/71‑1 (MQ=255)
                                        tgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCAcgt                                      <  1:5649444/60‑3 (MQ=255)
                                        tgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccata                           >  1:33357/1‑70 (MQ=255)
                                         gCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                          >  1:4892227/1‑71 (MQ=255)
                                            gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcca                             <  1:3370554/65‑1 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                 <  1:1075610/44‑1 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                 <  1:5237625/44‑1 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                 <  1:4838559/44‑1 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                 <  1:4813594/44‑1 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                 <  1:469002/44‑1 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtcc                      >  1:2943055/1‑71 (MQ=255)
                                             ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCATGCAg                                                 <  1:4093900/44‑1 (MQ=255)
                                              taccTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcct                     >  1:2604965/3‑70 (MQ=255)
                                              tGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtccg                     >  1:2437893/1‑71 (MQ=255)
                                               gCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtccg                     >  1:153873/1‑70 (MQ=255)
                                               gCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgt                        <  1:4345819/67‑1 (MQ=255)
                                               accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                    >  1:3405713/2‑69 (MQ=255)
                                               accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                    >  1:5384115/2‑69 (MQ=255)
                                               accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                    >  1:3270859/2‑69 (MQ=255)
                                               accTTTCGCCAGCATGATCATGCCGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcct                     >  1:5142767/2‑69 (MQ=255)
                                                 cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggcg                  >  1:2273606/1‑71 (MQ=255)
                                                  tttCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtcca                             <  1:4624087/59‑1 (MQ=255)
                                                   ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                >  1:3373890/1‑71 (MQ=255)
                                                    tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa                                       >  1:1007613/1‑47 (MQ=255)
                                                    tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa                                       >  1:2347324/1‑47 (MQ=255)
                                                    tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                <  1:1482870/70‑1 (MQ=255)
                                                    tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc               >  1:4617686/1‑71 (MQ=255)
                                                      gCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCggtggt             <  1:208187/71‑1 (MQ=255)
                                                      gCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGtg               <  1:2335407/69‑2 (MQ=255)
                                                      gCCAGCATGATCATACTTCCGCCCTGCGCTAGCAGTTGATCGACATAAGaa                                 >  1:179943/1‑51 (MQ=255)
                                                       ccAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGtt            <  1:5185377/71‑1 (MQ=255)
                                                       ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatc                               <  1:3114174/52‑1 (MQ=255)
                                                       ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccg                         <  1:5585939/58‑1 (MQ=255)
                                                       ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt             <  1:3198675/70‑1 (MQ=255)
                                                         aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg          >  1:1936462/1‑71 (MQ=255)
                                                         aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg          <  1:3163240/71‑1 (MQ=255)
                                                           cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt             <  1:5507567/66‑1 (MQ=255)
                                                           cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc        >  1:2295277/1‑71 (MQ=255)
                                                             tgttCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCa      <  1:3217070/68‑1 (MQ=255)
                                                             tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                     <  1:469015/40‑1 (MQ=255)
                                                             tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                     <  1:2032453/40‑1 (MQ=255)
                                                             tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                     <  1:4895398/40‑1 (MQ=255)
                                                             tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg           >  1:1444583/1‑66 (MQ=255)
                                                             tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc       <  1:4603298/70‑1 (MQ=255)
                                                              gATCATACTTCCGCCCTGTGCTTGCAGTTGATCGACATAAGaa                                 <  1:1828306/43‑1 (MQ=255)
                                                              gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt            <  1:2871881/64‑1 (MQ=255)
                                                               aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                          >  1:3933305/1‑49 (MQ=255)
                                                               aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                          >  1:2968097/1‑49 (MQ=255)
                                                               aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                <  1:3905659/43‑1 (MQ=255)
                                                               aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                <  1:430467/43‑1 (MQ=255)
                                                               aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                <  1:4334084/43‑1 (MQ=255)
                                                               aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc       >  1:520481/1‑68 (MQ=255)
                                                                tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                   <  1:5372804/39‑1 (MQ=255)
                                                                tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                   <  1:1640576/39‑1 (MQ=255)
                                                                tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                   <  1:1201293/39‑1 (MQ=255)
                                                                tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg              >  1:4616898/1‑60 (MQ=255)
                                                                  aTACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAagag  >  1:5288592/1‑70 (MQ=255)
                                                                    aCTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                   >  1:276670/1‑35 (MQ=255)
                                                                    aCTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                   >  1:917123/1‑35 (MQ=255)
                                                                    |                                                                   
AAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAG  >  minE/1014630‑1014765

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: