Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,783 G→C 29.8% I157M (ATC→ATG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,7830GC29.8% 167.9 / 87.5 114I157M (ATC→ATGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (55/25);  new base C (17/17);  total (72/42)
Fisher's exact test for biased strand distribution p-value = 8.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 121 total aligned reads displayed.

ACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTG  >  minE/609714‑609852
                                                                     |                                                                     
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:1634714/1‑71 (MQ=255)
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:830147/1‑71 (MQ=255)
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:1910296/1‑71 (MQ=255)
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:5506485/1‑71 (MQ=255)
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:2137121/1‑71 (MQ=255)
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:3452840/1‑71 (MQ=255)
aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAgg                                                                       <  1:438047/70‑1 (MQ=255)
  cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                                                    <  1:269486/71‑1 (MQ=255)
   gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                                                                   >  1:840490/1‑69 (MQ=255)
    atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                                  >  1:417548/1‑71 (MQ=255)
    atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                                  >  1:3351202/1‑71 (MQ=255)
       atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                <  1:1230781/70‑1 (MQ=255)
       atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                                >  1:5603563/1‑70 (MQ=255)
          cgcccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAt                                                                     >  1:291098/4‑62 (MQ=255)
          cgcccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAt                                                                     >  1:2696991/4‑62 (MQ=255)
             ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGTCAACGCAGGATGTCCGGAGTGTa                                                         <  1:4426013/71‑1 (MQ=255)
             ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:119365/71‑1 (MQ=255)
              cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                        <  1:609488/71‑1 (MQ=255)
              cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:472761/70‑1 (MQ=255)
              cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                        <  1:5234997/71‑1 (MQ=255)
               acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                         <  1:826443/69‑1 (MQ=255)
               acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                         <  1:5222851/69‑1 (MQ=255)
                cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:1699264/1‑70 (MQ=255)
                cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                      >  1:5312540/1‑71 (MQ=255)
                cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                      >  1:2896784/1‑71 (MQ=255)
                 ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                     <  1:114108/71‑1 (MQ=255)
                 ggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAca                                                     <  1:5066423/71‑2 (MQ=255)
                   ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATc                                                   <  1:4805222/71‑1 (MQ=255)
                    caTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                  <  1:1392433/71‑1 (MQ=255)
                      tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                                     <  1:1389149/50‑1 (MQ=255)
                      tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:2781513/62‑1 (MQ=255)
                      tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                <  1:1112119/71‑1 (MQ=255)
                      tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                <  1:4442888/71‑1 (MQ=255)
                       agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                >  1:2467907/4‑70 (MQ=255)
                        tGCATCTGGCAGACCGGGCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:3822322/1‑71 (MQ=255)
                        gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                       >  1:35780/3‑62 (MQ=255)
                        gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaat                                               <  1:4228793/68‑2 (MQ=255)
                            tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATcc                                                                 <  1:2465098/48‑1 (MQ=16)
                            tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATcc                                                                 <  1:4288415/48‑1 (MQ=16)
                            tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggtg                                                             <  1:3678607/52‑3 (MQ=16)
                            tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                         <  1:685436/56‑1 (MQ=255)
                            tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                         <  1:145809/56‑1 (MQ=255)
                             cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                         >  1:4247352/1‑71 (MQ=255)
                              tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:1846544/54‑1 (MQ=255)
                              tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                        >  1:5075076/1‑71 (MQ=255)
                              tGGCAGACCGGTCAGAACACCAGATTTACGGAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                       >  1:2475998/1‑71 (MQ=255)
                                 cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                >  1:3868989/1‑60 (MQ=255)
                                 cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                >  1:4702682/1‑60 (MQ=255)
                                  aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                     >  1:1024750/3‑70 (MQ=255)
                                   gACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                                       <  1:1134571/35‑1 (MQ=255)
                                   gACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                                       >  1:1265131/1‑35 (MQ=255)
                                     ccGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                                     >  1:2787592/1‑35 (MQ=255)
                                      cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatgg                                >  1:3710304/1‑70 (MQ=255)
                                       ggTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:3069347/45‑1 (MQ=255)
                                        gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                         >  1:4575174/1‑60 (MQ=12)
                                        gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                         >  1:4581481/1‑60 (MQ=12)
                                        gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                         >  1:4106300/1‑60 (MQ=12)
                                        gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                         >  1:5571530/1‑60 (MQ=12)
                                        gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:1149188/1‑55 (MQ=255)
                                        gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAac                                               >  1:1333273/1‑54 (MQ=255)
                                         tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:2025101/1‑45 (MQ=255)
                                         tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:3094220/1‑45 (MQ=255)
                                          cAGAACACCAGATTTACGGCAACGCATGATGTCCGGAGTGTAAACGTCg                                                  <  1:417847/49‑1 (MQ=255)
                                           aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                            >  1:4540961/1‑70 (MQ=255)
                                            gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTTAACGTCGaa                                                >  1:4909086/1‑49 (MQ=255)
                                            gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                >  1:2810451/1‑49 (MQ=255)
                                            gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                >  1:1354427/1‑49 (MQ=255)
                                             aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:3661589/1‑50 (MQ=255)
                                             aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                          >  1:1366394/1‑70 (MQ=255)
                                                 ccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGtt                             <  1:1475341/63‑1 (MQ=16)
                                                 ccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACgcag                     <  1:5526446/71‑4 (MQ=255)
                                                  cAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                        >  1:4448415/1‑67 (MQ=255)
                                                   aGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                        >  1:5020873/1‑66 (MQ=255)
                                                   aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                   >  1:4220959/1‑71 (MQ=255)
                                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:3553148/1‑43 (MQ=255)
                                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:425310/1‑43 (MQ=255)
                                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                   >  1:5335648/1‑70 (MQ=255)
                                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                  <  1:2987674/71‑1 (MQ=255)
                                                     aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGt                              <  1:1816085/58‑1 (MQ=12)
                                                      tttAAGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                >  1:5282951/1‑39 (MQ=255)
                                                        tACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                        >  1:624521/1‑61 (MQ=16)
                                                        tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                <  1:718323/53‑1 (MQ=255)
                                                        tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                <  1:5670167/53‑1 (MQ=255)
                                                         aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg              >  1:4481659/1‑70 (MQ=255)
                                                            gCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                <  1:257130/49‑1 (MQ=255)
                                                             cAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                           >  1:1238475/1‑53 (MQ=255)
                                                              aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt        >  1:3581036/1‑71 (MQ=255)
                                                               aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg              >  1:2324190/1‑64 (MQ=255)
                                                                cGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt      <  1:1350574/71‑1 (MQ=255)
                                                                 gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                          >  1:1039979/1‑50 (MQ=255)
                                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                        >  1:1739280/1‑51 (MQ=255)
                                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                        >  1:2460056/1‑51 (MQ=255)
                                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                        >  1:683568/1‑51 (MQ=255)
                                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                        >  1:2670501/1‑51 (MQ=255)
                                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt    >  1:1047864/1‑71 (MQ=255)
                                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt     <  1:5433830/70‑1 (MQ=255)
                                                                    ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAg                               <  1:3102581/42‑1 (MQ=255)
                                                                    ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCGGTGAAGATTTTTTTg  <  1:1513943/71‑1 (MQ=255)
                                                                    ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg  <  1:3189165/71‑1 (MQ=255)
                                                                     gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                          <  1:3490918/46‑1 (MQ=255)
                                                                     |                                                                     
ACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTG  >  minE/609714‑609852

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