Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,895 T→C 12.0% E453E (GAA→GAG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,8950TC12.0% 194.8 / 15.7 75E453E (GAA→GAG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (38/28);  new base C (2/7);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 7.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCA  >  minE/608825‑608965
                                                                      |                                                                      
tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                        <  1:2916003/71‑1 (MQ=255)
tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                        <  1:1473380/71‑1 (MQ=255)
tCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGtttt                                                                        >  1:3446003/1‑71 (MQ=255)
  aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                      <  1:38364/71‑1 (MQ=255)
  aCTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                      <  1:2151196/71‑1 (MQ=255)
      catcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTc                                                                   >  1:2737071/1‑70 (MQ=255)
       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                                    <  1:119104/68‑1 (MQ=255)
       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                  <  1:2478310/70‑1 (MQ=255)
       atcatAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCa                                                                  <  1:1936985/70‑1 (MQ=255)
          atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                              >  1:1906577/1‑71 (MQ=255)
                  ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTc                                                                      >  1:2616954/1‑55 (MQ=255)
                  ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                       >  1:3316256/1‑70 (MQ=255)
                     cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                   <  1:3893848/71‑1 (MQ=255)
                     cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                   >  1:266883/1‑71 (MQ=255)
                     cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                   <  1:3944075/71‑1 (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)
                              tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGt                                                                    >  1:3864525/1‑45 (MQ=255)
                              tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                           >  1:3720852/1‑70 (MQ=255)
                              tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                           >  1:401538/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:941882/1‑71 (MQ=255)
                                       cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                 >  1:479043/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:1660998/1‑71 (MQ=255)
                                              cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                          >  1:2065075/1‑71 (MQ=255)
                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                          >  1:2579454/1‑71 (MQ=255)
                                              cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                          >  1:2574866/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:3646875/1‑51 (MQ=255)
                                                              ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                              >  1:3053695/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:3214221/1‑45 (MQ=255)
                                                               tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                   >  1:569910/1‑45 (MQ=255)
                                                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                            >  1:2858225/1‑35 (MQ=255)
                                                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                            >  1:659043/1‑35 (MQ=255)
                                                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag         <  1:675343/70‑1 (MQ=255)
                                                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag         >  1:2215457/1‑70 (MQ=255)
                                                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                         <  1:3144932/53‑1 (MQ=255)
                                                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa       >  1:3825802/1‑71 (MQ=255)
                                                                 aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa       >  1:922123/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:3152497/1‑71 (MQ=255)
                                                                      ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa  >  1:3789051/1‑71 (MQ=255)
                                                                      |                                                                      
TCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCA  >  minE/608825‑608965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑