Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,631 G→A 6.8% Y146Y (TAC→TAT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6310GA 6.8% 235.9 / 5.4 73Y146Y (TAC→TATfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base A (3/2);  total (38/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 = 8.79e-01

GATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCG  >  minE/1015561‑1015700
                                                                      |                                                                     
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: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:329373/71‑1 (MQ=255)
                       gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                                <  1:4114088/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:891835/1‑66 (MQ=255)
                          ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                  >  1:2358300/1‑66 (MQ=255)
                          ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                          >  1:1658889/1‑58 (MQ=255)
                          ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                             <  1:3687280/71‑1 (MQ=255)
                          ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                             <  1:3399309/71‑1 (MQ=255)
                          ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                             <  1:3889355/71‑1 (MQ=255)
                          ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                             >  1:4129317/1‑71 (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)
                                  cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAg                                                                    >  1:1727413/1‑40 (MQ=255)
                                  cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAg                                                                    >  1:3101627/1‑40 (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)
                                    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:3850327/1‑42 (MQ=255)
                                       cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                             >  1:653114/1‑42 (MQ=255)
                                       cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                             >  1:2928814/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)
                                                                      |                                                                     
GATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCG  >  minE/1015561‑1015700

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: