Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,701 T→G 18.3% G456G (GGA→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7010TG18.3% 145.9 / 38.5 82G456G (GGA→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base G (9/6);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCG  >  minE/1014634‑1014770
                                                                   |                                                                     
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:2412482/70‑1 (MQ=255)
            tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                                                         >  1:118167/1‑70 (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)
                       gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                                    <  1:5680458/48‑1 (MQ=255)
                       gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                                                                    <  1:5604358/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:4813594/44‑1 (MQ=255)
                                         ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                      <  1:469002/44‑1 (MQ=255)
                                         ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                      <  1:4838559/44‑1 (MQ=255)
                                         ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                      <  1:5237625/44‑1 (MQ=255)
                                         ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAg                                                      <  1:1075610/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:3270859/2‑69 (MQ=255)
                                           accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                         >  1:5384115/2‑69 (MQ=255)
                                           accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                         >  1:3405713/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:4895398/40‑1 (MQ=255)
                                                         tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                          <  1:2032453/40‑1 (MQ=255)
                                                         tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                          <  1:469015/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:2968097/1‑49 (MQ=255)
                                                           aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                               >  1:3933305/1‑49 (MQ=255)
                                                           aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                     <  1:430467/43‑1 (MQ=255)
                                                           aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                     <  1:3905659/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)
                                                                   tCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg  >  1:3303521/1‑70 (MQ=255)
                                                                   |                                                                     
CGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCG  >  minE/1014634‑1014770

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: