Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,625 T→C 6.4% Q148Q (CAA→CAG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6250TC 6.4% 224.5 / 6.0 78Q148Q (CAA→CAGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (3/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCA  >  minE/1015557‑1015694
                                                                    |                                                                     
tAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGtttt                                                                      <  1:2381082/70‑1 (MQ=255)
tAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTg                                                                     <  1:3783138/71‑1 (MQ=255)
 aaaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTg                                                                     >  1:1336362/1‑70 (MQ=255)
  aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                    >  1:821474/1‑70 (MQ=255)
  aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                    >  1:3794899/1‑70 (MQ=255)
  aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                    <  1:3462149/70‑1 (MQ=255)
  aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                    >  1:2418447/1‑70 (MQ=255)
    gatcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAg                                                                 <  1:2445750/71‑1 (MQ=255)
      tcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAg                                                                 >  1:2144326/1‑69 (MQ=255)
           cTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAgcg                                                            >  1:1017995/1‑69 (MQ=255)
              cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAgc                                                             <  1:1992101/65‑1 (MQ=255)
              cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCa                                                          <  1:914653/68‑1 (MQ=255)
              cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                        >  1:1831867/1‑70 (MQ=255)
               gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                      <  1:1557821/71‑1 (MQ=255)
               gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                      <  1:2949303/71‑1 (MQ=255)
                  ggCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATc                                                   >  1:3230298/1‑71 (MQ=255)
                           gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                           <  1:2029589/70‑1 (MQ=255)
                           gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                           <  1:592169/70‑1 (MQ=255)
                           gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                           <  1:1561873/70‑1 (MQ=255)
                           gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                          <  1:4114088/71‑1 (MQ=255)
                           gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                          <  1:329373/71‑1 (MQ=255)
                            tGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                       >  1:2106549/1‑57 (MQ=255)
                             gCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                        >  1:2795268/1‑55 (MQ=255)
                             gCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGa                                              >  1:2969842/1‑65 (MQ=255)
                              ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTTTGAATAGCGCAGGTTATCTTCGATATAGgt                                       <  1:3211866/71‑1 (MQ=255)
                              ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                            >  1:2358300/1‑66 (MQ=255)
                              ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                            >  1:891835/1‑66 (MQ=255)
                              ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                    <  1:747230/42‑1 (MQ=255)
                              ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                    >  1:1658889/1‑58 (MQ=255)
                              ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                       <  1:3889355/71‑1 (MQ=255)
                              ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                       >  1:4129317/1‑71 (MQ=255)
                              ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                       <  1:3687280/71‑1 (MQ=255)
                              ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                       <  1:3399309/71‑1 (MQ=255)
                                ggTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                     <  1:3894255/71‑1 (MQ=255)
                                 gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                     >  1:2778410/1‑70 (MQ=255)
                                  tATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                        >  1:1052967/1‑50 (MQ=255)
                                  tATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                   <  1:3613026/71‑1 (MQ=255)
                                    ttCACTTCTTTATACATATCGAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCg                                               <  1:756792/57‑1 (MQ=255)
                                    ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCg                                                                   <  1:33198/37‑1 (MQ=255)
                                    ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                        >  1:2574574/1‑48 (MQ=255)
                                    ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                           <  1:464675/61‑1 (MQ=255)
                                    ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                     <  1:3743346/67‑1 (MQ=255)
                                    ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                                 <  1:1575410/71‑1 (MQ=255)
                                      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAg                                                              >  1:3101627/1‑40 (MQ=255)
                                      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAg                                                              >  1:1727413/1‑40 (MQ=255)
                                      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                   >  1:2090403/1‑67 (MQ=255)
                                      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                   >  1:2127447/1‑67 (MQ=255)
                                      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                               >  1:2286823/1‑71 (MQ=255)
                                       aCTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGac                              <  1:2048024/71‑1 (MQ=255)
                                        cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTa                                                               >  1:2909165/1‑37 (MQ=255)
                                        cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCTATATAAgt                                       <  1:3722959/61‑1 (MQ=255)
                                        cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                   <  1:1927646/65‑1 (MQ=255)
                                         ttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAg                                                         <  1:3751995/42‑1 (MQ=255)
                                         ttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGa                                              <  1:1069832/53‑1 (MQ=255)
                                         ttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                     >  1:63227/1‑62 (MQ=255)
                                           cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                       >  1:2928814/1‑42 (MQ=255)
                                           cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                       >  1:653114/1‑42 (MQ=255)
                                           cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                       >  1:3850327/1‑42 (MQ=255)
                                           cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATATTCGATATAAGTGTTATAg                                >  1:1107220/1‑65 (MQ=255)
                                            tttATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                         >  1:566649/1‑71 (MQ=255)
                                              tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                      >  1:2825859/1‑40 (MQ=255)
                                              tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                    >  1:1710526/1‑42 (MQ=255)
                                              tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgcgc                       >  1:669124/1‑71 (MQ=255)
                                                taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                      <  1:964407/38‑1 (MQ=255)
                                                taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                      <  1:1978949/38‑1 (MQ=255)
                                                taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                     <  1:2887829/71‑1 (MQ=255)
                                                taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                     <  1:1374200/71‑1 (MQ=255)
                                                    tatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                     <  1:3827516/51‑1 (MQ=255)
                                                    tatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgc                 >  1:1916275/1‑71 (MQ=255)
                                                     atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCtt                                                 >  1:570182/1‑38 (MQ=255)
                                                      tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATaa                                         >  1:2666517/1‑45 (MQ=255)
                                                      tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                   <  1:1341398/51‑1 (MQ=255)
                                                       ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc              >  1:4017755/1‑71 (MQ=255)
                                                        cAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                            >  1:3736814/1‑40 (MQ=16)
                                                            ggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                          >  1:3346754/1‑38 (MQ=255)
                                                            ggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAg                  >  1:3534261/1‑62 (MQ=255)
                                                             ccgcaTTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                    <  1:227307/42‑1 (MQ=16)
                                                              cgcgTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                          <  1:2691694/36‑1 (MQ=255)
                                                               gcgTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtcatt      <  1:3193000/71‑2 (MQ=255)
                                                                  ttttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcac    >  1:898395/1‑70 (MQ=255)
                                                                   tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                               >  1:903594/1‑42 (MQ=255)
                                                                   tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca  >  1:3572308/1‑71 (MQ=255)
                                                                    |                                                                     
TAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCA  >  minE/1015557‑1015694

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