Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,640 A→G 8.9% N143N (AAT→AAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6400AG 8.9% 185.5 / 8.9 67N143N (AAT→AACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base G (3/3);  total (32/35)
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.74e-01

CGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCA  >  minE/1015571‑1015705
                                                                     |                                                                 
cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                                   >  1:1831867/1‑70 (MQ=255)
 gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                                 <  1:2949303/71‑1 (MQ=255)
 gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                                 <  1:1557821/71‑1 (MQ=255)
    ggCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATc                                                              >  1:3230298/1‑71 (MQ=255)
             gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                      <  1:592169/70‑1 (MQ=255)
             gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                      <  1:1561873/70‑1 (MQ=255)
             gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                                      <  1:2029589/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)
                ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                               >  1:1658889/1‑58 (MQ=255)
                ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                                  <  1:3889355/71‑1 (MQ=255)
                ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                                  <  1:3399309/71‑1 (MQ=255)
                ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                                  >  1:4129317/1‑71 (MQ=255)
                ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                                  <  1:3687280/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)
                      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)
                        cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                              >  1:2127447/1‑67 (MQ=255)
                        cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                              >  1:2090403/1‑67 (MQ=255)
                        cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                          >  1:2286823/1‑71 (MQ=255)
                         aCTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGac                                         <  1:2048024/71‑1 (MQ=255)
                          cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCTATATAAgt                                                  <  1:3722959/61‑1 (MQ=255)
                          cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                              <  1:1927646/65‑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)
                                                         cGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccaccac         <  1:490723/71‑1 (MQ=255)
                                                           aGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCg       <  1:369095/71‑1 (MQ=255)
                                                            gTAGCGCAGATTATCTTCGATATAAGTGTTATAGacac                                       >  1:1402853/1‑38 (MQ=255)
                                                            gTAGCGCAGATTATCTTCGATATAAGTGTTATAGaca                                        <  1:2867547/37‑1 (MQ=255)
                                                             tAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGt     <  1:3475443/71‑1 (MQ=255)
                                                               gcgcAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTcc   <  1:1394865/71‑1 (MQ=255)
                                                                 gcAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCa  >  1:163965/1‑70 (MQ=255)
                                                                     |                                                                 
CGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCA  >  minE/1015571‑1015705

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: