Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,949 C→T 20.8% A435A (GCG→GCA)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9490CT20.8% 216.2 / 47.3 101A435A (GCG→GCA) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base T (12/9);  total (52/49)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCG  >  minE/608879‑609018
                                                                      |                                                                     
ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                                                                       <  1:2996212/71‑1 (MQ=255)
 cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                                      <  1:3665939/71‑1 (MQ=255)
    cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                                    <  1:2835929/70‑1 (MQ=255)
    cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                                                                   >  1:1781102/1‑71 (MQ=255)
      ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                 <  1:921355/71‑1 (MQ=255)
          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                              <  1:675343/70‑1 (MQ=255)
          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                              >  1:2215457/1‑70 (MQ=255)
           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                            >  1:922123/1‑71 (MQ=255)
           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                            >  1:3825802/1‑71 (MQ=255)
              ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                         >  1:3185564/1‑71 (MQ=255)
               tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                         >  1:3422021/1‑70 (MQ=255)
                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                                                                >  1:1259392/1‑62 (MQ=255)
                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                       >  1:3152497/1‑71 (MQ=255)
                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                       >  1:3789051/1‑71 (MQ=255)
                 tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                      >  1:311948/1‑71 (MQ=255)
                  cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                      <  1:973447/70‑1 (MQ=255)
                  cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                      <  1:1852187/70‑1 (MQ=255)
                   gTCCACGCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                     <  1:1468722/70‑1 (MQ=255)
                   gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                     <  1:2778101/70‑1 (MQ=255)
                     cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                  >  1:206006/1‑71 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                 <  1:2912248/71‑1 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                 >  1:2708089/1‑71 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                 >  1:3542483/1‑71 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                  <  1:2048483/70‑1 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                  <  1:3233016/70‑1 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                  <  1:2327653/70‑1 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGG‑TTGCACGCGCAcc                                                                  >  1:489912/1‑53 (MQ=255)
                       aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                                <  1:2576364/71‑1 (MQ=255)
                        cgccgcCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                                             <  1:2662450/57‑1 (MQ=255)
                         cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGAGCGTGCACCAAAGAACTGCATATgc                                                  >  1:682230/2‑66 (MQ=255)
                           cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                             >  1:1643088/1‑70 (MQ=255)
                            gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                                           >  1:2542308/1‑71 (MQ=255)
                             ccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                              <  1:3901684/51‑1 (MQ=255)
                             ccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                          <  1:3946145/71‑1 (MQ=255)
                               gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                      <  1:2242866/57‑1 (MQ=255)
                               gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                                        >  1:2423550/1‑71 (MQ=255)
                               gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                                        >  1:3502839/1‑71 (MQ=255)
                               gTTGATTGCGTACACCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                      <  1:2564288/57‑1 (MQ=255)
                                  gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                       >  1:1525207/1‑69 (MQ=255)
                                   aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                 <  1:823223/42‑1 (MQ=255)
                                   aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                 <  1:2077600/42‑1 (MQ=255)
                                   aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                    <  1:3354044/71‑1 (MQ=255)
                                   aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     <  1:3315875/70‑1 (MQ=255)
                                   aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     <  1:2226859/70‑1 (MQ=255)
                                   aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     <  1:3074143/70‑1 (MQ=255)
                                      gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                 <  1:2344491/39‑1 (MQ=255)
                                      gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCg                                  <  1:3931339/70‑1 (MQ=255)
                                       cGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                     >  1:148422/1‑50 (MQ=255)
                                       cGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                     >  1:3137387/1‑50 (MQ=255)
                                       cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCggg                                >  1:1808979/1‑71 (MQ=255)
                                        gTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                           <  1:1013124/43‑1 (MQ=255)
                                        gTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGc                               >  1:1335388/1‑71 (MQ=255)
                                         tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                              <  1:1903957/71‑1 (MQ=255)
                                           cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                            <  1:3093794/71‑1 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                 <  1:2632525/48‑1 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCggg                                <  1:1882474/65‑1 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                 >  1:3197307/1‑64 (MQ=255)
                                                tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                       >  1:2136562/1‑71 (MQ=255)
                                                gcgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                       <  1:574363/69‑1 (MQ=255)
                                                 cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                       <  1:1296101/69‑1 (MQ=255)
                                                 cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                       >  1:3598461/2‑70 (MQ=255)
                                                  gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAcgccgc                      <  1:1695153/70‑1 (MQ=255)
                                                  gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                     <  1:256287/71‑1 (MQ=255)
                                                  gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                      <  1:2876279/70‑1 (MQ=255)
                                                  gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                      >  1:3293049/1‑70 (MQ=255)
                                                   ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                     <  1:453039/70‑1 (MQ=255)
                                                   ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                        <  1:300387/35‑1 (MQ=16)
                                                   ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                        <  1:3319818/35‑1 (MQ=16)
                                                   ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcag                    >  1:3433081/1‑70 (MQ=255)
                                                    tttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                   >  1:3199830/1‑69 (MQ=255)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAGGATCATCGGGCTTACGCagcaag                   >  1:1778673/1‑70 (MQ=255)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                 >  1:3470676/1‑56 (MQ=255)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                   >  1:527550/1‑70 (MQ=255)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                   >  1:1305649/1‑70 (MQ=255)
                                                     ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                    >  1:1313612/1‑37 (MQ=21)
                                                      tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                  >  1:3820084/1‑70 (MQ=255)
                                                       cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                       <  1:351852/64‑1 (MQ=255)
                                                         ccAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa              <  1:2004439/71‑1 (MQ=255)
                                                         ccAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc              <  1:2743812/71‑2 (MQ=255)
                                                          cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     >  1:2448292/1‑47 (MQ=255)
                                                          cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     >  1:2690320/1‑47 (MQ=255)
                                                          cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     >  1:434905/1‑47 (MQ=255)
                                                          cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtcg             <  1:2406376/71‑3 (MQ=255)
                                                           aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTa                                               >  1:3778134/1‑36 (MQ=255)
                                                           aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc            >  1:746943/1‑71 (MQ=255)
                                                           aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc            >  1:1523287/1‑71 (MQ=255)
                                                            ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCATCAAGCAATAGCg           >  1:3409611/1‑71 (MQ=255)
                                                            ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg           >  1:584689/1‑71 (MQ=255)
                                                             gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                          <  1:2842939/39‑1 (MQ=255)
                                                              tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg           <  1:3466562/69‑1 (MQ=255)
                                                               ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg        >  1:3045171/1‑71 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg        <  1:1345190/70‑1 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt       >  1:1450037/1‑71 (MQ=255)
                                                                  cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa     >  1:1403558/1‑71 (MQ=255)
                                                                  cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa     >  1:1887477/1‑71 (MQ=255)
                                                                   aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                <  1:430575/59‑1 (MQ=255)
                                                                   aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                <  1:435288/59‑1 (MQ=255)
                                                                    cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc   <  1:3810445/71‑1 (MQ=255)
                                                                     gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg  <  1:2264441/71‑1 (MQ=255)
                                                                      cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                              >  1:238756/1‑42 (MQ=255)
                                                                      |                                                                     
CCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCG  >  minE/608879‑609018

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: ‑