Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,541 G→A 17.9% D176D (GAC→GAT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,5410GA17.9% 195.4 / 31.5 89D176D (GAC→GATfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base A (7/9);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATAT  >  minE/1015472‑1015611
                                                                     |                                                                      
ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                                                       >  1:5297930/1‑71 (MQ=255)
ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                                                       >  1:3400818/1‑71 (MQ=255)
ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                                                       >  1:422688/1‑71 (MQ=255)
ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcg                                                                        >  1:1547230/1‑70 (MQ=255)
 gTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                                                       <  1:4363897/70‑1 (MQ=255)
    tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGc                                                                    >  1:2245741/1‑70 (MQ=255)
    tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGc                                                                    <  1:3758264/70‑1 (MQ=255)
    tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgtc                                                                      <  1:3839973/68‑1 (MQ=255)
    tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc                                                                    <  1:4641737/70‑1 (MQ=255)
    tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc                                                                    >  1:2995995/1‑70 (MQ=255)
     ccTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCa                                                                  >  1:3948516/1‑71 (MQ=255)
       tggtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                                                 >  1:1893027/4‑70 (MQ=255)
       tggtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                                                 >  1:3978994/4‑70 (MQ=255)
         gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCTATc                                                                <  1:2453718/68‑1 (MQ=255)
         gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                                                                <  1:1750586/68‑1 (MQ=255)
         gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                                                                <  1:184854/68‑1 (MQ=255)
         aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                                                                <  1:4116209/69‑1 (MQ=255)
         aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                                                                <  1:1020268/69‑1 (MQ=255)
         aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                                              <  1:3230301/71‑1 (MQ=255)
         aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTAGTCGCCATc                                                                <  1:4432689/69‑1 (MQ=255)
          tAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                                                 >  1:1283257/1‑67 (MQ=255)
          tAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc                                                             >  1:1370218/1‑71 (MQ=255)
               aTACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCa                                                         >  1:2072584/1‑70 (MQ=255)
               aTACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCa                                                         >  1:2886802/1‑70 (MQ=255)
                tACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCAATCAACGGCAt                                                        >  1:2007269/1‑70 (MQ=255)
                 aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtcatcg                                                                     >  1:791211/1‑56 (MQ=255)
                 aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtcatcg                                                                     >  1:2748907/1‑56 (MQ=255)
                 aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATcta                                                              >  1:2418078/1‑61 (MQ=255)
                  cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTaca                                                     <  1:4390952/71‑3 (MQ=255)
                    ttttgttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAgg                                                   >  1:2328173/3‑70 (MQ=255)
                      ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       >  1:5291368/1‑65 (MQ=255)
                      ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga                                                   <  1:1196073/69‑1 (MQ=255)
                        gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                                              >  1:4162793/1‑56 (MQ=255)
                        gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                                              >  1:4382572/1‑56 (MQ=255)
                         tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat                                              >  1:360859/1‑71 (MQ=255)
                         tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat                                              >  1:4978529/1‑71 (MQ=255)
                         tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcga                                               >  1:453718/1‑70 (MQ=255)
                              cGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcg                                         <  1:4692434/71‑1 (MQ=255)
                                aaCCACCACCTTTGGCGATACAGAGGAATTTGTATtcg                                                                        >  1:3157211/1‑38 (MQ=255)
                                aaCCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgatc                                             <  1:743428/65‑1 (MQ=255)
                                 accaccaccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTg                                      <  1:725818/71‑1 (MQ=255)
                                  ccaccaccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTgg                                     <  1:4156516/71‑1 (MQ=255)
                                       accTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCt                                            <  1:2720264/59‑1 (MQ=255)
                                        ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGc                                                          >  1:2795932/1‑44 (MQ=255)
                                        ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       <  1:4584255/47‑1 (MQ=255)
                                        ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       <  1:4458752/47‑1 (MQ=255)
                                        ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAtt                               <  1:4524809/71‑1 (MQ=255)
                                        ccTTTCGCGATACAGAGGAATTTGTATTCGTCGCCa                                                                  >  1:1318878/1‑36 (MQ=255)
                                         cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgat                                                  >  1:2622213/1‑51 (MQ=255)
                                         cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgat                                                  >  1:3159724/1‑51 (MQ=255)
                                         cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTg                              <  1:3775936/71‑1 (MQ=255)
                                          tttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTgg                             >  1:2072133/1‑71 (MQ=255)
                                          tttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTgg                             >  1:3325390/1‑71 (MQ=255)
                                           ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                                              >  1:2488213/1‑37 (MQ=255)
                                           ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga                                                   >  1:4013716/1‑48 (MQ=255)
                                           ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAg                                                    <  1:4629740/47‑1 (MQ=255)
                                           ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGt                            >  1:1891320/1‑71 (MQ=255)
                                           ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGt                            >  1:4099552/1‑71 (MQ=255)
                                           ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGt                            >  1:226183/1‑71 (MQ=255)
                                                gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       <  1:2108055/39‑1 (MQ=255)
                                                gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       <  1:4176970/39‑1 (MQ=255)
                                                gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       <  1:4971077/39‑1 (MQ=255)
                                                gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                                                       <  1:5315601/39‑1 (MQ=255)
                                                 aTACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTg                              <  1:1310120/63‑1 (MQ=255)
                                                  tACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcgc                                        <  1:5610401/52‑1 (MQ=255)
                                                  tACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTa                     >  1:4344834/1‑71 (MQ=255)
                                                  tACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTa                     >  1:3207396/1‑71 (MQ=255)
                                                    caaagGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGtt                   <  1:4681343/68‑1 (MQ=255)
                                                    caaagGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGtt                   <  1:1091046/68‑1 (MQ=255)
                                                    caaagGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGtt                   <  1:1514313/68‑1 (MQ=255)
                                                    cAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTAtt                   <  1:5620944/71‑1 (MQ=255)
                                                    cAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTAtt                   <  1:731829/71‑1 (MQ=255)
                                                    cAGAGGAATTTGTATTCATCGCCATCAACGGCATAAAGATCGATCTGCGCTg                                      >  1:249488/1‑52 (MQ=255)
                                                     agagGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTc                  <  1:2360023/71‑1 (MQ=255)
                                                       agGAATTTGTACTCATCGCCATCTACCGCTTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGtt                   <  1:396478/68‑1 (MQ=255)
                                                        ggAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAtt                               >  1:2818779/1‑55 (MQ=255)
                                                        ggAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAtt                               >  1:4419217/1‑55 (MQ=255)
                                                        ggAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAtt                               >  1:4584288/1‑55 (MQ=255)
                                                        ggAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCAct               <  1:3024302/71‑1 (MQ=255)
                                                          aaTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcgc                                        >  1:13705/1‑44 (MQ=255)
                                                          aaTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcgc                                        >  1:2143245/1‑44 (MQ=255)
                                                          aaTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGa                                 <  1:2050158/51‑1 (MQ=255)
                                                          aaTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGa                                 <  1:5270838/51‑1 (MQ=255)
                                                          aaTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGa                                 <  1:5111803/51‑1 (MQ=255)
                                                          aaTTCGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcgc                                        >  1:4555774/1‑44 (MQ=255)
                                                           aTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCa                                   <  1:554684/48‑1 (MQ=255)
                                                           aTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCa                                   <  1:5671698/48‑1 (MQ=255)
                                                           aTTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCa                                   <  1:781674/48‑1 (MQ=255)
                                                            tttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCa                                   <  1:1503289/47‑1 (MQ=255)
                                                            tttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAt                                <  1:1423350/50‑1 (MQ=255)
                                                            tttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCActtctt           <  1:4432141/71‑1 (MQ=255)
                                                             ttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCActtctt           >  1:4243759/1‑70 (MQ=255)
                                                             ttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCActtctt           >  1:4250465/1‑70 (MQ=255)
                                                              tGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTg                              >  1:2239319/1‑50 (MQ=255)
                                                                tATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTtata       >  1:142132/1‑71 (MQ=255)
                                                                tATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTtat        <  1:5629182/70‑1 (MQ=255)
                                                                 aTTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTtata       >  1:4939793/1‑70 (MQ=255)
                                                                   tcgtcgCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCAc                >  1:2745223/1‑59 (MQ=255)
                                                                   tcgtcgCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACa     >  1:2891749/1‑70 (MQ=255)
                                                                     gtcgCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACatat  <  1:1752028/71‑1 (MQ=255)
                                                                     |                                                                      
GGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATAT  >  minE/1015472‑1015611

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: