Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,946 A→G 20.2% R436R (CGT→CGC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9460AG20.2% 212.8 / 48.8 104R436R (CGT→CGC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (42/41);  new base G (12/9);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

Only 100 of 108 total aligned reads displayed.

ACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCA  >  minE/608878‑609015
                                                                    |                                                                     
aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                     >  1:1475953/1‑71 (MQ=255)
 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)
         ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                                                             <  1:3563169/70‑1 (MQ=255)
         ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                              >  1:3240407/1‑69 (MQ=255)
           cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                           >  1:2215457/1‑70 (MQ=255)
           cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                           <  1:675343/70‑1 (MQ=255)
            aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                         >  1:922123/1‑71 (MQ=255)
            aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                         >  1:3825802/1‑71 (MQ=255)
              tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                      >  1:570867/1‑56 (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:3789051/1‑71 (MQ=255)
                 ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                    >  1:3152497/1‑71 (MQ=255)
                  tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                   >  1:311948/1‑71 (MQ=255)
                   cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                   <  1:1852187/70‑1 (MQ=255)
                    gTCCACGCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                  <  1:1468722/70‑1 (MQ=255)
                    gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                  <  1:2778101/70‑1 (MQ=255)
                     tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                                 <  1:2557423/54‑1 (MQ=255)
                      cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                                      <  1:1149445/48‑1 (MQ=255)
                      cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                               >  1:206006/1‑71 (MQ=255)
                       aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                              <  1:2912248/71‑1 (MQ=255)
                       aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                              >  1:2708089/1‑71 (MQ=255)
                       aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                               <  1:2048483/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)
                            cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgt                                                                    >  1:2889951/1‑43 (MQ=16)
                            cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                          >  1:1643088/1‑70 (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:3502839/1‑71 (MQ=255)
                                gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                                     >  1:2423550/1‑71 (MQ=255)
                                gTTGATTGCGTACACCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                   <  1:2564288/57‑1 (MQ=255)
                                   gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                    >  1:1525207/1‑69 (MQ=255)
                                    aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                              <  1:2077600/42‑1 (MQ=255)
                                    aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                              <  1:823223/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: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:1305649/1‑70 (MQ=255)
                                                      ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                >  1:527550/1‑70 (MQ=255)
                                                      ttCGCCAGGTTTGCACGCGCACCGAACAACTGCAt                                                   >  1:2062714/1‑35 (MQ=255)
                                                      ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                 >  1:1313612/1‑37 (MQ=21)
                                                       tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc               >  1:3820084/1‑70 (MQ=255)
                                                        cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                    <  1:351852/64‑1 (MQ=255)
                                                          ccAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc           <  1:2743812/71‑2 (MQ=255)
                                                           cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                  >  1:2690320/1‑47 (MQ=255)
                                                           cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                  >  1:2448292/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)
                                                             ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg        >  1:1617089/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:1887477/1‑71 (MQ=255)
                                                                   cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa  >  1:1403558/1‑71 (MQ=255)
                                                                    aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa             <  1:430575/59‑1 (MQ=255)
                                                                    aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa             <  1:435288/59‑1 (MQ=255)
                                                                    |                                                                     
ACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCA  >  minE/608878‑609015

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