Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,955 G→A 25.2% F433F (TTC→TTT)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9550GA25.2% 178.1 / 60.8 99F433F (TTC→TTT) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base A (13/12);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 114 total aligned reads displayed.

TTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGA  >  minE/608885‑609025
                                                                      |                                                                      
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)
     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:3789051/1‑71 (MQ=255)
          ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                              >  1:3152497/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:2708089/1‑71 (MQ=255)
                aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                        >  1:3542483/1‑71 (MQ=255)
                aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                         <  1:3233016/70‑1 (MQ=255)
                aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                         <  1:2048483/70‑1 (MQ=255)
                aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                                                         <  1:2327653/70‑1 (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)
                         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)
                                          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)
                                            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:3319818/35‑1 (MQ=16)
                                             ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                               <  1:300387/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)
                                               ttCGCCAGGTTTGCACGCGCACCGAACAACTGCAt                                                             >  1:2062714/1‑35 (MQ=255)
                                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                           >  1:1313612/1‑37 (MQ=21)
                                                tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                         >  1:3820084/1‑70 (MQ=255)
                                                   ccAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                     <  1:2004439/71‑1 (MQ=255)
                                                   ccAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                     <  1:2743812/71‑2 (MQ=255)
                                                    cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                            >  1:437652/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)
                                                            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)
                                                              cgtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTc          <  1:3689070/68‑1 (MQ=255)
                                                              cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc          <  1:2905324/71‑1 (MQ=255)
                                                               gtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTaa               >  1:3461263/3‑64 (MQ=255)
                                                               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg         <  1:2264441/71‑1 (MQ=255)
                                                                cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                     >  1:238756/1‑42 (MQ=255)
                                                                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt       <  1:2275509/71‑1 (MQ=255)
                                                                  cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg      <  1:444757/71‑1 (MQ=255)
                                                                   aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgt     <  1:1484372/71‑1 (MQ=255)
                                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                              >  1:3632553/1‑44 (MQ=255)
                                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg   <  1:404777/71‑1 (MQ=255)
                                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg   <  1:76445/71‑1 (MQ=255)
                                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg   >  1:2659473/1‑70 (MQ=255)
                                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa  <  1:2835979/71‑1 (MQ=255)
                                                                      |                                                                      
TTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGA  >  minE/608885‑609025

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