Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,010 C→G 11.3% A353A (GCG→GCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0100CG11.3% 192.1 / 13.5 71A353A (GCG→GCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base G (5/3);  total (40/31)
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.85e-01

CTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAC  >  minE/1014941‑1015080
                                                                     |                                                                      
cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCg                                                                       >  1:1089664/1‑71 (MQ=255)
cTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCg                                                                       >  1:1534982/1‑71 (MQ=255)
 tGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGc                                                                      >  1:1858761/1‑71 (MQ=255)
  gCGAAAGCTGGGCGAGGATCTCTTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGcc                                                                     >  1:3034778/1‑70 (MQ=255)
   cGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGc                                                                      <  1:1757436/69‑1 (MQ=255)
    gACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCttt                                                                   >  1:704309/1‑71 (MQ=255)
       aaCTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTgcg                                                                >  1:591585/1‑71 (MQ=255)
         cTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAAGGCGCACCGCTTCGCCTTCTCCCGCTTTgcg                                                                >  1:2612827/1‑69 (MQ=255)
          tgtgCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCtt                                                                    >  1:193981/1‑64 (MQ=255)
               gAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTgcgc                                                               <  1:3066474/64‑1 (MQ=255)
                   aTCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTccc                                                                        <  1:1930734/51‑1 (MQ=255)
                     ctctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAg                                                             <  1:2703865/60‑1 (MQ=255)
                     ctctTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTc                                                       >  1:818005/1‑66 (MQ=255)
                      tctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAg                                                             <  1:1604314/59‑1 (MQ=255)
                       ctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAgtt                                                           >  1:1701736/1‑58 (MQ=255)
                       ctTTCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAgtt                                                           >  1:4069418/1‑58 (MQ=255)
                       ctTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGat                                                 >  1:2833588/1‑70 (MQ=255)
                        tttCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGctc                                                          >  1:3367253/1‑60 (MQ=255)
                        tttCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCgg                                                    <  1:2316345/66‑1 (MQ=255)
                          tCATCGGGCGGTTAAGGTCAACTTTCACCGCTTCGCCTTCACCGGCCTGGCGCAgtt                                                           >  1:1688759/1‑55 (MQ=255)
                          tCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatata                                              <  1:2174847/70‑1 (MQ=255)
                          tCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGata                                                <  1:4284/68‑1 (MQ=255)
                           cATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATAtt                                            >  1:1361354/1‑71 (MQ=255)
                            aTCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATAttt                                           <  1:4089249/71‑1 (MQ=255)
                              cGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGc                                         >  1:1349687/1‑71 (MQ=255)
                               ggACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGcc                                        >  1:3936515/1‑71 (MQ=255)
                                gACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCt                                       >  1:328254/1‑71 (MQ=255)
                                 aCGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTg                                      >  1:612937/1‑71 (MQ=255)
                                  cGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTcc                                                      >  1:2893468/1‑54 (MQ=255)
                                    gTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTgcg                                                                <  1:658069/42‑1 (MQ=255)
                                    gTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                   <  1:1378592/71‑1 (MQ=255)
                                     ttAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAtt                                  <  1:1909983/71‑1 (MQ=255)
                                      tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatat                                               >  1:3177009/1‑57 (MQ=255)
                                      tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                 >  1:899932/1‑71 (MQ=255)
                                      tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                 >  1:2330732/1‑71 (MQ=255)
                                       aaGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                 <  1:141893/70‑1 (MQ=255)
                                        aGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTg                                                                  <  1:2378082/36‑1 (MQ=255)
                                        aGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCa                                                              >  1:237372/1‑40 (MQ=255)
                                          gTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatata                                              >  1:1888592/1‑54 (MQ=255)
                                            cAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTc                            >  1:2154091/1‑70 (MQ=255)
                                             aaCGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGcc                                        <  1:2509229/57‑1 (MQ=255)
                                             aaCGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                   <  1:1811078/62‑1 (MQ=255)
                                               cgcgCACCGCTTCGCCTTCTCCAGCTTTTAGCCGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGt                        >  1:68083/1‑71 (MQ=255)
                                                gcgcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatata                                              <  1:1960902/48‑1 (MQ=255)
                                                 ttcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                      >  1:1674801/3‑71 (MQ=255)
                                                  gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                     >  1:165342/1‑71 (MQ=255)
                                                  gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                     <  1:1504462/71‑1 (MQ=255)
                                                  gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                     <  1:1440185/71‑1 (MQ=255)
                                                  gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                     <  1:2237799/71‑1 (MQ=255)
                                                  gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                     <  1:2405387/71‑1 (MQ=255)
                                                    aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTg                                      <  1:865198/52‑1 (MQ=255)
                                                    aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                   >  1:99078/1‑71 (MQ=255)
                                                    aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                   <  1:869609/71‑1 (MQ=255)
                                                     ccGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGata                                                <  1:3080140/41‑1 (MQ=255)
                                                      cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCg                  >  1:423275/1‑70 (MQ=255)
                                                      cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                 >  1:3842882/1‑71 (MQ=255)
                                                      cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                 >  1:2705091/1‑71 (MQ=255)
                                                        cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatat                                               <  1:3500533/39‑1 (MQ=255)
                                                        cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGatat                                               <  1:4109325/39‑1 (MQ=255)
                                                        cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc               <  1:3447399/71‑1 (MQ=255)
                                                        cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc               <  1:3260101/71‑1 (MQ=255)
                                                           cGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCa             <  1:326395/70‑1 (MQ=255)
                                                            gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                 <  1:1672012/49‑1 (MQ=255)
                                                              cTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCt                                       >  1:2978340/1‑41 (MQ=255)
                                                               ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                   <  1:4140745/44‑1 (MQ=255)
                                                               ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc               <  1:3598271/64‑1 (MQ=255)
                                                               ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcc        >  1:274429/1‑71 (MQ=255)
                                                                tctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATccc       >  1:417230/1‑71 (MQ=255)
                                                                 ctcCCGCTTTGCGGAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                    >  1:3477920/1‑57 (MQ=255)
                                                                 ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                    >  1:727207/1‑57 (MQ=255)
                                                                 ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                    >  1:3289995/1‑57 (MQ=255)
                                                                 ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcccc      >  1:2102682/1‑71 (MQ=255)
                                                                    ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                >  1:1201972/1‑58 (MQ=255)
                                                                    ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc               <  1:1437471/59‑1 (MQ=255)
                                                                    ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc               <  1:1035314/59‑1 (MQ=255)
                                                                    ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTg    <  1:1254966/70‑1 (MQ=255)
                                                                     cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc  >  1:1914875/1‑71 (MQ=255)
                                                                     cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc  >  1:1971624/1‑71 (MQ=255)
                                                                     |                                                                      
CTGCGACAACTGTGCGAGGATCTCTTTCATCGGACGGTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAC  >  minE/1014941‑1015080

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: