Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,807 C→T 9.7% V149V (GTG→GTA pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,8070CT 9.7% 233.9 / 11.2 82V149V (GTG→GTApflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base T (6/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTG  >  minE/609738‑609876
                                                                     |                                                                     
tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      <  1:3713503/71‑1 (MQ=255)
tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      <  1:3863764/71‑1 (MQ=255)
tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      <  1:1557305/71‑1 (MQ=255)
  cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACg                                                                    <  1:3294142/71‑1 (MQ=255)
   aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGc                                                                   >  1:934271/1‑71 (MQ=255)
     cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:3124160/1‑71 (MQ=255)
     cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:683712/1‑71 (MQ=255)
     cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:1200076/1‑71 (MQ=255)
     cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:3202735/1‑71 (MQ=255)
           aaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGtta                                                           >  1:3413941/2‑70 (MQ=255)
            aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                          >  1:1962760/1‑71 (MQ=255)
            aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                          <  1:794064/71‑1 (MQ=255)
             ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                                                         <  1:1080095/71‑1 (MQ=255)
              cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                        <  1:472274/71‑1 (MQ=255)
               ggTCAGAACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCACTGGTtatggg                                                       >  1:2884392/1‑69 (MQ=255)
                gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                        >  1:2894459/1‑69 (MQ=255)
                 tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGGCGAACACGCCCTGGTTGTGAGtt                                                     <  1:1754112/71‑1 (MQ=255)
                  cAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                    >  1:10591/1‑71 (MQ=255)
                   aGAACTCCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                                   <  1:2585645/71‑1 (MQ=255)
                   aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTTTAAACGTCGAACACGCCCTGGGTGTGAGtttt                                                   <  1:3339867/71‑1 (MQ=255)
                   aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                    <  1:858471/70‑1 (MQ=255)
                   aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                    <  1:2490538/70‑1 (MQ=255)
                   aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                                    <  1:186373/70‑1 (MQ=255)
                     aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                                  <  1:816445/70‑1 (MQ=255)
                      acgccAGATTTACGGCATCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                                >  1:1815777/4‑71 (MQ=255)
                       cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                               >  1:437384/1‑71 (MQ=255)
                       cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                               >  1:1129451/1‑71 (MQ=255)
                        gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                                                    <  1:1631026/64‑1 (MQ=255)
                        gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                                >  1:2061190/2‑69 (MQ=255)
                          cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                                            >  1:1531336/1‑71 (MQ=255)
                           aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                                            >  1:1452631/1‑54 (MQ=255)
                            gATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                                >  1:1326924/1‑65 (MQ=255)
                            gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                                                    >  1:1028313/1‑45 (MQ=255)
                            gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGt                                              >  1:1099139/1‑67 (MQ=255)
                            gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                                           >  1:593872/1‑70 (MQ=255)
                            gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                                          <  1:1668607/71‑1 (MQ=255)
                               ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                                               >  1:2056657/1‑47 (MQ=255)
                               ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                       >  1:3031529/1‑71 (MQ=255)
                               ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                       >  1:1708068/1‑71 (MQ=255)
                                tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                       <  1:2080340/70‑1 (MQ=255)
                                 aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                      <  1:2314118/70‑1 (MQ=255)
                                 aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                      <  1:3194530/70‑1 (MQ=255)
                                  cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGa                                     <  1:1607697/70‑1 (MQ=255)
                                      aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt                                >  1:1050678/1‑71 (MQ=255)
                                       aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:649819/1‑37 (MQ=255)
                                       aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:1502410/1‑37 (MQ=255)
                                       aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGa                                  >  1:3117557/1‑68 (MQ=255)
                                         gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAg                                                       <  1:82200/45‑1 (MQ=255)
                                         gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                             <  1:1930890/71‑1 (MQ=255)
                                          cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                            >  1:1831887/1‑71 (MQ=255)
                                           aGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt                           >  1:1797344/1‑71 (MQ=255)
                                            ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                          <  1:1372368/71‑1 (MQ=255)
                                             gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                         <  1:2074140/71‑1 (MQ=255)
                                             gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                         <  1:2008984/71‑1 (MQ=255)
                                             gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                         <  1:2937464/71‑1 (MQ=255)
                                              aTGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                        >  1:3690858/1‑71 (MQ=255)
                                               tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatc                       <  1:1569361/71‑1 (MQ=255)
                                               tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                        >  1:2696638/1‑70 (MQ=255)
                                                  ccGGAGTGTAAACGTCTAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                       <  1:1123984/52‑1 (MQ=255)
                                                   cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc                    <  1:3007321/70‑1 (MQ=255)
                                                   cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc                    <  1:392250/70‑1 (MQ=255)
                                                   cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                   <  1:418031/71‑1 (MQ=255)
                                                   cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                   >  1:622611/1‑71 (MQ=255)
                                                    ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg                  <  1:2932642/71‑1 (MQ=255)
                                                    ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg                  <  1:351944/71‑1 (MQ=255)
                                                    ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg                  <  1:3602155/71‑1 (MQ=255)
                                                     gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                                        <  1:1069834/48‑1 (MQ=255)
                                                     gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa                 >  1:985011/1‑71 (MQ=255)
                                                        tgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATcc              <  1:2847058/71‑1 (MQ=255)
                                                          tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                              <  1:3734233/53‑1 (MQ=255)
                                                          tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg            >  1:3697829/1‑71 (MQ=255)
                                                          tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg            >  1:2603706/1‑71 (MQ=255)
                                                            aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGtt          >  1:3641481/1‑71 (MQ=255)
                                                             aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa                 >  1:3677054/1‑63 (MQ=255)
                                                             aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc         >  1:1798891/1‑71 (MQ=255)
                                                             aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc         >  1:3148315/1‑71 (MQ=255)
                                                              cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                             <  1:1078831/50‑1 (MQ=255)
                                                              cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg        <  1:3790389/71‑1 (MQ=255)
                                                              cGTCGAACACGCCCTGGTTGTGAGATTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg        >  1:1820147/1‑71 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTAGGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcgc       >  1:1263605/1‑71 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcgc       >  1:253230/1‑71 (MQ=255)
                                                               gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg        >  1:3811134/1‑70 (MQ=255)
                                                                  gAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                          <  1:664459/49‑1 (MQ=255)
                                                                  gAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                          <  1:1816992/49‑1 (MQ=255)
                                                                  gAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                          <  1:1402262/49‑1 (MQ=255)
                                                                    acacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTg  <  1:2620820/71‑1 (MQ=255)
                                                                     cacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCgg     >  1:1398152/1‑67 (MQ=255)
                                                                     |                                                                     
TGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTG  >  minE/609738‑609876

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: