Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,747 C→T 21.2% L169L (CTG→CTA) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,7470CT21.2% 206.7 / 48.2 99L169L (CTG→CTA) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (46/32);  new base T (11/10);  total (57/42)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCC  >  minE/609678‑609813
                                                                     |                                                                  
gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCa                                                                   <  1:2268395/71‑1 (MQ=255)
  aGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                                                 >  1:878512/1‑71 (MQ=255)
   gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                                                 <  1:1835672/70‑1 (MQ=255)
   gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                                                >  1:1429495/1‑71 (MQ=255)
   gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                                                >  1:2999190/1‑71 (MQ=255)
      gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                             >  1:772263/1‑71 (MQ=255)
      gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                             <  1:651234/71‑1 (MQ=255)
      gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                             >  1:1975028/1‑71 (MQ=255)
      gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCg                                                              <  1:3947610/70‑1 (MQ=255)
        tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                                                >  1:1150230/1‑66 (MQ=255)
         aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                                             >  1:3559219/1‑68 (MQ=255)
         aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                                                            >  1:2029271/1‑69 (MQ=255)
          ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGAATCTGGCAGACCGGTCa                                                          <  1:2138738/70‑1 (MQ=255)
           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)
                                    aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCa                                                          >  1:2631320/1‑44 (MQ=255)
                                    aCCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAgc                                                        <  1:3941157/46‑2 (MQ=255)
                                     ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                              <  1:1036060/71‑1 (MQ=255)
                                     ccGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGt                                                            <  1:2146723/41‑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)
                                        atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                                           >  1:1744646/1‑39 (MQ=255)
                                        atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                                           >  1:2591318/1‑39 (MQ=255)
                                        atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                   >  1:1713932/1‑63 (MQ=255)
                                        atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                   >  1:563614/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:403082/1‑69 (MQ=255)
                                            tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                         >  1:1819363/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:2393813/1‑70 (MQ=255)
                                                    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                >  1:2683014/1‑70 (MQ=255)
                                                    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                >  1:2389073/1‑70 (MQ=255)
                                                    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                >  1:3015164/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:3667396/1‑56 (MQ=255)
                                                           atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                       >  1:2681569/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:622401/1‑65 (MQ=255)
                                                               aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa          >  1:3472631/1‑65 (MQ=255)
                                                                tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                 >  1:2598261/1‑41 (MQ=19)
                                                                tCTGGCAGACCGGTCAGAACACCAGATTTACGGCa                                       <  1:398088/35‑1 (MQ=255)
                                                                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                >  1:430826/1‑57 (MQ=255)
                                                                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                >  1:75521/1‑57 (MQ=255)
                                                                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc  >  1:683712/1‑71 (MQ=255)
                                                                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc  >  1:3124160/1‑71 (MQ=255)
                                                                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc  >  1:1200076/1‑71 (MQ=255)
                                                                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc  >  1:3202735/1‑71 (MQ=255)
                                                                   ggCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa          <  1:3232434/61‑1 (MQ=255)
                                                                    gCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa          <  1:2694772/60‑1 (MQ=255)
                                                                     |                                                                  
GTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCC  >  minE/609678‑609813

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