Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,949 C→T 17.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,9490CT17.8% 220.5 / 35.7 90A435A (GCG→GCA) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (5/11);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCG  >  minE/608879‑609018
                                                                      |                                                                     
ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                                                                       >  1:2043102/1‑71 (MQ=255)
  tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                                      >  1:24167/1‑70 (MQ=255)
    cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACgcgc                                                                     <  1:1895898/69‑1 (MQ=255)
        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                 >  1:3654931/1‑69 (MQ=255)
        ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                                 >  1:2381497/1‑69 (MQ=255)
        ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                                 >  1:1494808/1‑68 (MQ=255)
          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                              <  1:102482/70‑1 (MQ=255)
          cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                              <  1:27440/70‑1 (MQ=255)
          cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                                                <  1:1916967/68‑3 (MQ=255)
          cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                                                                <  1:3865116/68‑3 (MQ=255)
          cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                                 <  1:3618961/67‑2 (MQ=255)
          cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                              <  1:1401060/70‑1 (MQ=255)
          cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                              <  1:1703505/70‑1 (MQ=255)
           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                            >  1:595649/1‑71 (MQ=255)
           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                                             <  1:2728383/70‑1 (MQ=255)
           aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                                             <  1:3160660/70‑1 (MQ=255)
            gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                           <  1:2532648/71‑1 (MQ=255)
             tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGAACGCGCACCGAAGAACt                                                          >  1:2008479/1‑71 (MQ=255)
              ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                         >  1:3852655/1‑71 (MQ=255)
              ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                         >  1:3491701/1‑71 (MQ=255)
               tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                         >  1:3339379/1‑70 (MQ=255)
                ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                              >  1:229808/1‑64 (MQ=255)
                 tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                       <  1:2971132/70‑1 (MQ=255)
                  cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGATCTGCAt                                                      <  1:1270/70‑1 (MQ=255)
                   gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                      <  1:52892/69‑1 (MQ=255)
                    tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                                                   <  1:2163221/71‑1 (MQ=255)
                      aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                 >  1:1856739/1‑71 (MQ=255)
                       aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                                <  1:3815239/71‑1 (MQ=255)
                        ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                               <  1:2023686/69‑1 (MQ=255)
                        ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                               <  1:725150/69‑1 (MQ=255)
                         gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                                              >  1:2708355/1‑71 (MQ=255)
                          ccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaa                                                               >  1:2585727/1‑53 (MQ=255)
                            gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                                           >  1:128065/1‑71 (MQ=255)
                              cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                          >  1:3075564/1‑70 (MQ=255)
                              cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                          >  1:953683/1‑70 (MQ=255)
                              cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                         <  1:677044/71‑1 (MQ=255)
                                 tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                       >  1:657558/1‑70 (MQ=255)
                                  gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     >  1:587614/1‑71 (MQ=255)
                                  gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                     >  1:1852293/1‑71 (MQ=255)
                                  gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                      <  1:4005748/70‑1 (MQ=255)
                                  gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTCGTTTACCAACGatca                                     >  1:685199/1‑71 (MQ=255)
                                    ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                   <  1:910646/71‑1 (MQ=255)
                                    ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                   <  1:1823578/71‑1 (MQ=255)
                                      gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                                            >  1:3238564/1‑44 (MQ=12)
                                      gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                                    >  1:2117140/1‑36 (MQ=255)
                                          aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                    >  1:111947/1‑64 (MQ=255)
                                          aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat                                    >  1:1450022/1‑64 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                              >  1:1288974/1‑67 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                          <  1:2648335/71‑1 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                          <  1:43335/71‑1 (MQ=255)
                                             gCATGGTTTTCGCCAGGTTTGCACCCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                          >  1:3004311/1‑71 (MQ=255)
                                             gCATGGGTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAttt                                                    >  1:538428/1‑45 (MQ=255)
                                                tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                            <  1:1360880/50‑1 (MQ=255)
                                                tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                 <  1:4123172/61‑1 (MQ=255)
                                                tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                 <  1:2636176/61‑1 (MQ=255)
                                                tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                 <  1:1040688/61‑1 (MQ=255)
                                                tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCGTACg                          >  1:3760010/1‑68 (MQ=255)
                                                gcgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                       >  1:1398399/3‑71 (MQ=255)
                                                 ggTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                                       >  1:119737/1‑38 (MQ=255)
                                                  gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                     <  1:3498025/71‑1 (MQ=255)
                                                   ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGc                                                        <  1:408750/35‑1 (MQ=255)
                                                   ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                                <  1:3826659/43‑1 (MQ=255)
                                                   ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                        <  1:2918905/67‑1 (MQ=255)
                                                   ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa                                           <  1:3893971/48‑1 (MQ=25)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                                      >  1:3826947/1‑35 (MQ=255)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                      <  1:1990208/51‑1 (MQ=255)
                                                     ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                      >  1:3624840/1‑51 (MQ=255)
                                                      tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                  <  1:3395187/70‑1 (MQ=255)
                                                       cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                 <  1:11636/70‑1 (MQ=255)
                                                       cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                 <  1:1670833/70‑1 (MQ=255)
                                                       cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                 <  1:1199145/70‑1 (MQ=255)
                                                       cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggca                 <  1:1900468/70‑1 (MQ=255)
                                                          cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg             <  1:1472220/71‑1 (MQ=255)
                                                          cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGctt                                                <  1:2241829/36‑1 (MQ=16)
                                                           aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc            >  1:4119712/1‑71 (MQ=255)
                                                            ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg             >  1:4078426/1‑69 (MQ=255)
                                                            ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg             >  1:2729084/1‑69 (MQ=255)
                                                            ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg           <  1:2076321/71‑1 (MQ=255)
                                                            ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg           >  1:1221440/1‑71 (MQ=255)
                                                             gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGt          <  1:1555201/71‑1 (MQ=255)
                                                              tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTa         <  1:232203/71‑1 (MQ=255)
                                                               ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                >  1:1117979/1‑63 (MQ=255)
                                                               ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa              >  1:2022989/1‑65 (MQ=255)
                                                               ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg        >  1:40416/1‑71 (MQ=255)
                                                               ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg        >  1:593652/1‑71 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg        >  1:488209/1‑70 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt       >  1:3324975/1‑71 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt       >  1:509491/1‑71 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt       >  1:2525102/1‑71 (MQ=255)
                                                                tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt       >  1:2004154/1‑71 (MQ=255)
                                                                  cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                             <  1:1916882/47‑1 (MQ=255)
                                                                   aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                  <  1:3363694/57‑1 (MQ=255)
                                                                   aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAt    <  1:3381657/71‑1 (MQ=255)
                                                                     gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                >  1:2363300/1‑57 (MQ=255)
                                                                     gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                >  1:282008/1‑57 (MQ=255)
                                                                     gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg  <  1:2744558/71‑1 (MQ=255)
                                                                      |                                                                     
CCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCG  >  minE/608879‑609018

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