Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,759 A→C 20.8% V165V (GTT→GTG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,7590AC20.8% 149.3 / 55.5 96V165V (GTT→GTG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base C (12/8);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTA  >  minE/609689‑609828
                                                                      |                                                                     
cGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGa                                                                       <  1:614699/71‑1 (MQ=255)
 gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                                      <  1:2199452/71‑1 (MQ=255)
    cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                                      <  1:1353960/68‑1 (MQ=255)
     aGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGcc                                                                  <  1:3794751/71‑1 (MQ=255)
        gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGa                                                               <  1:2344681/68‑1 (MQ=255)
           aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACAcc                                                                  <  1:3513508/65‑1 (MQ=255)
           aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAttt                                                            <  1:3367157/71‑1 (MQ=255)
            cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                                           <  1:3746778/71‑1 (MQ=255)
            cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                                                           >  1:1437964/1‑71 (MQ=255)
                aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg                                                        >  1:3381941/1‑70 (MQ=255)
                 cGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                                       >  1:2347202/1‑70 (MQ=255)
                  ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                                                     >  1:277157/1‑71 (MQ=255)
                    tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                                                   <  1:3363619/71‑1 (MQ=255)
                      gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCg                                                   >  1:1254477/1‑69 (MQ=255)
                       tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                             <  1:2711657/58‑1 (MQ=255)
                       tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                          <  1:3535560/61‑1 (MQ=255)
                       tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                                <  1:2118094/71‑1 (MQ=255)
                          ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                             <  1:1036060/71‑1 (MQ=255)
                           cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                >  1:2603756/1‑67 (MQ=255)
                           cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                >  1:762822/1‑67 (MQ=255)
                           cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                            >  1:1591625/1‑71 (MQ=255)
                           cGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGAtt                                                             <  1:1316290/54‑1 (MQ=255)
                            gatgatACGGCCACGGCCATATGGATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                           >  1:3618208/1‑71 (MQ=255)
                            gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                           >  1:2490979/1‑71 (MQ=255)
                            gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                                           <  1:3795941/71‑3 (MQ=255)
                             atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                  >  1:563614/1‑63 (MQ=255)
                             atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                  >  1:1713932/1‑63 (MQ=255)
                             atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                          >  1:8738/1‑71 (MQ=255)
                             atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAAGAtag                                           >  1:1777735/1‑68 (MQ=255)
                                atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                       <  1:2110042/71‑1 (MQ=255)
                                atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                        <  1:264934/70‑1 (MQ=255)
                                 tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                        >  1:1819363/1‑69 (MQ=255)
                                 tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                        >  1:403082/1‑69 (MQ=255)
                                    ggccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                     <  1:3154099/69‑1 (MQ=255)
                                     gccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCa                                                 <  1:681140/56‑1 (MQ=255)
                                     cccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAgc                                               >  1:2368655/2‑57 (MQ=255)
                                      ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgtgt                                  >  1:2180042/1‑70 (MQ=255)
                                      ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgt                                    >  1:2812195/1‑68 (MQ=255)
                                       cacggccaTATGCATGTGGCAGACCGGTCAGAACACCAGAtt                                                             <  1:1626566/42‑1 (MQ=255)
                                       cacggccaTATGCATCTGTCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                <  1:5535/71‑1 (MQ=255)
                                       cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                             <  1:3645305/42‑1 (MQ=255)
                                       cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                <  1:3006292/71‑1 (MQ=255)
                                        acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                               >  1:1205796/1‑71 (MQ=255)
                                         cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                 >  1:2006983/1‑68 (MQ=255)
                                         cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                 <  1:1759126/68‑1 (MQ=255)
                                         cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                               >  1:3015164/1‑70 (MQ=255)
                                         cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                               >  1:2389073/1‑70 (MQ=255)
                                         cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                               >  1:2683014/1‑70 (MQ=255)
                                         cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                               >  1:2393813/1‑70 (MQ=255)
                                          tgccaGATGCATATGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                             >  1:2894211/2‑71 (MQ=255)
                                          ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                             >  1:96175/1‑71 (MQ=255)
                                          ggccaTATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                             <  1:1154638/71‑1 (MQ=255)
                                            ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAcat                            >  1:3254109/1‑68 (MQ=255)
                                               tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg                                                        <  1:735503/39‑1 (MQ=255)
                                               tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                <  1:1244351/47‑1 (MQ=255)
                                                atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                      >  1:2681569/1‑56 (MQ=255)
                                                atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                      >  1:3667396/1‑56 (MQ=255)
                                                 tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                               >  1:89121/1‑62 (MQ=255)
                                                 tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                      <  1:3713503/71‑1 (MQ=255)
                                                 tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                      <  1:1557305/71‑1 (MQ=255)
                                                 tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                      <  1:3863764/71‑1 (MQ=255)
                                                  gCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                      >  1:2082483/1‑54 (MQ=255)
                                                  ccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggt                                      >  1:1971649/2‑53 (MQ=16)
                                                   cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACg                    <  1:3294142/71‑1 (MQ=255)
                                                   cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                >  1:2722572/1‑59 (MQ=255)
                                                   cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                >  1:3389602/1‑59 (MQ=255)
                                                    aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATcaa                         <  1:2573065/65‑3 (MQ=255)
                                                    aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGc                   >  1:934271/1‑71 (MQ=255)
                                                    aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                     >  1:84688/1‑53 (MQ=255)
                                                    aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                         >  1:3472631/1‑65 (MQ=255)
                                                    aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                         >  1:622401/1‑65 (MQ=255)
                                                     tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                                >  1:2598261/1‑41 (MQ=19)
                                                     tCTGGCAGACCGGTCAGAACACCAGATTTACGGCa                                                      <  1:398088/35‑1 (MQ=255)
                                                      cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                               >  1:75521/1‑57 (MQ=255)
                                                      cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                               >  1:430826/1‑57 (MQ=255)
                                                      cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                 >  1:3124160/1‑71 (MQ=255)
                                                      cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                 >  1:3202735/1‑71 (MQ=255)
                                                      cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                 >  1:1200076/1‑71 (MQ=255)
                                                      cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                 >  1:683712/1‑71 (MQ=255)
                                                        ggCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                         <  1:3232434/61‑1 (MQ=255)
                                                         gCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                         <  1:2694772/60‑1 (MQ=255)
                                                           aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                 <  1:1952353/50‑1 (MQ=255)
                                                            aaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGtta           >  1:3413941/2‑70 (MQ=255)
                                                             aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                        <  1:540759/57‑1 (MQ=255)
                                                             aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt          >  1:1962760/1‑71 (MQ=255)
                                                             aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt          <  1:794064/71‑1 (MQ=255)
                                                              ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                              >  1:3332771/1‑50 (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)
                                                                    aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                          <  1:3287562/48‑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)
                                                                      |                                                                     
CGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTA  >  minE/609689‑609828

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