Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,904 G→C 11.9% G450G (GGC→GGG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9040GC11.9% 243.4 / 20.0 92G450G (GGC→GGG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (48/33);  new base C (3/8);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTT  >  minE/608835‑608972
                                                                     |                                                                    
atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                                     >  1:1906577/1‑71 (MQ=255)
        ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                              >  1:3316256/1‑70 (MQ=255)
           cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                          <  1:3893848/71‑1 (MQ=255)
           cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                          <  1:3944075/71‑1 (MQ=255)
           cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                          >  1:266883/1‑71 (MQ=255)
               gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                       >  1:3545287/1‑70 (MQ=255)
               gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                       >  1:3283948/1‑70 (MQ=255)
               gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                       >  1:685629/1‑70 (MQ=255)
                 cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                    <  1:2993798/71‑1 (MQ=255)
                    tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                  >  1:401538/1‑70 (MQ=255)
                    tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                  >  1:3720852/1‑70 (MQ=255)
                     gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGc                                                 <  1:3535708/70‑1 (MQ=255)
                       tCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                              >  1:934303/1‑71 (MQ=255)
                        cGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                             >  1:3949552/1‑71 (MQ=255)
                         ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                            >  1:2842195/1‑71 (MQ=255)
                         ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                            <  1:798967/71‑1 (MQ=255)
                           ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                                     <  1:111933/44‑1 (MQ=255)
                           ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                          <  1:3929044/71‑1 (MQ=255)
                             cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                         <  1:2184315/70‑1 (MQ=255)
                             cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                        >  1:743499/1‑71 (MQ=255)
                             cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                        >  1:479043/1‑71 (MQ=255)
                             cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                        >  1:941882/1‑71 (MQ=255)
                               gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                           <  1:2513308/50‑1 (MQ=255)
                               gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                      <  1:2220937/71‑1 (MQ=255)
                                 cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                          >  1:3858894/1‑49 (MQ=255)
                                 cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                             >  1:2574848/1‑62 (MQ=255)
                                   tCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                   <  1:2889250/70‑1 (MQ=255)
                                    cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                 >  1:2065075/1‑71 (MQ=255)
                                    cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                 >  1:1660998/1‑71 (MQ=255)
                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                 >  1:2574866/1‑71 (MQ=255)
                                    cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                 >  1:2579454/1‑71 (MQ=255)
                                     ggACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCATGtt                                <  1:3209922/71‑1 (MQ=255)
                                     ggACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                 <  1:1811441/70‑1 (MQ=255)
                                      gACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                          >  1:3177881/1‑44 (MQ=255)
                                        ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGc                             <  1:485992/68‑1 (MQ=255)
                                          gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                           <  1:2869311/70‑3 (MQ=255)
                                          aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                             <  1:3071286/37‑1 (MQ=255)
                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                           >  1:2564407/1‑38 (MQ=255)
                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                  >  1:1137796/1‑63 (MQ=255)
                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                 >  1:2764278/1‑64 (MQ=255)
                                           aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                          >  1:1475953/1‑71 (MQ=255)
                                           aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                   <  1:2463701/62‑1 (MQ=255)
                                           aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                          <  1:1521300/71‑1 (MQ=255)
                                            ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                         <  1:2996212/71‑1 (MQ=255)
                                            ccTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                 <  1:3059050/63‑1 (MQ=255)
                                             cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                        <  1:3665939/71‑1 (MQ=255)
                                                cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                      <  1:2835929/70‑1 (MQ=255)
                                                cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                     >  1:1781102/1‑71 (MQ=255)
                                                  ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                         <  1:833394/49‑1 (MQ=255)
                                                  ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                   <  1:921355/71‑1 (MQ=255)
                                                   cttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                 <  1:3800219/55‑1 (MQ=255)
                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                             <  1:3159143/43‑1 (MQ=255)
                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                     >  1:3053695/1‑51 (MQ=255)
                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                     >  1:3646875/1‑51 (MQ=255)
                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                     >  1:1858069/1‑51 (MQ=255)
                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                  <  1:3563169/70‑1 (MQ=255)
                                                    ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                   >  1:3240407/1‑69 (MQ=255)
                                                     tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                          >  1:569910/1‑45 (MQ=255)
                                                     tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                          >  1:3214221/1‑45 (MQ=255)
                                                      cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                                   >  1:2858225/1‑35 (MQ=255)
                                                      cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                                   >  1:659043/1‑35 (MQ=255)
                                                      cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                >  1:2215457/1‑70 (MQ=255)
                                                      cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                <  1:675343/70‑1 (MQ=255)
                                                       aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                                <  1:3144932/53‑1 (MQ=255)
                                                       aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa              >  1:922123/1‑71 (MQ=255)
                                                       aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa              >  1:3825802/1‑71 (MQ=255)
                                                        gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGttt                               <  1:2342304/53‑1 (MQ=255)
                                                         tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                              >  1:3416380/1‑53 (MQ=255)
                                                         tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                              >  1:2451095/1‑53 (MQ=255)
                                                         tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                           >  1:570867/1‑56 (MQ=255)
                                                          ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg           >  1:3185564/1‑71 (MQ=255)
                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCa                                    >  1:1163429/1‑45 (MQ=255)
                                                           tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg           >  1:3422021/1‑70 (MQ=255)
                                                            ttCGTCAACGCCGCCGTTGATTGCGTACATCATGGTTTTCGcc                                     <  1:1492905/43‑1 (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)
                                                                tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                      <  1:2557423/54‑1 (MQ=255)
                                                                 cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                           <  1:1149445/48‑1 (MQ=255)
                                                                 cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt    >  1:206006/1‑71 (MQ=255)
                                                                  aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCa                                    >  1:2346010/1‑38 (MQ=255)
                                                                  aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt   <  1:2912248/71‑1 (MQ=255)
                                                                  aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt   >  1:3542483/1‑71 (MQ=255)
                                                                  aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt   >  1:2708089/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)
                                                                     |                                                                    
ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTT  >  minE/608835‑608972

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