Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 982,608 A→C 19.4% P18P (CCA→CCC)  ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE982,6080AC19.4% 281.8 / 84.3 160P18P (CCA→CCC) ydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (61/68);  new base C (11/20);  total (72/88)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 161 total aligned reads displayed.

GGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCC  >  minE/982539‑982677
                                                                     |                                                                     
gggAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                                                      <  1:4421474/71‑1 (MQ=255)
   aaaTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                                                   <  1:139504/71‑1 (MQ=255)
      tGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAgg                                                                <  1:4847719/71‑1 (MQ=255)
       gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                                               >  1:3614632/1‑71 (MQ=255)
       gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                                               >  1:3253683/1‑71 (MQ=255)
       gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGAGTGTTTCCAAACCAGGa                                                               >  1:3639645/1‑71 (MQ=255)
        aGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                              >  1:560113/1‑71 (MQ=255)
         ggTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACg                                                             <  1:254718/71‑1 (MQ=255)
             aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                         >  1:4293902/1‑71 (MQ=255)
             aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                         >  1:4859633/1‑71 (MQ=255)
              gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                                        <  1:1306891/71‑1 (MQ=255)
              gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTAAGGCGTGTTTCCAAACCAGGACGGCACg                                                        <  1:3760797/71‑1 (MQ=255)
               cATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGt                                                       >  1:1875911/1‑71 (MQ=255)
                aTGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                                                      <  1:1904906/71‑1 (MQ=255)
                 tGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                                                   <  1:2114/57‑1 (MQ=255)
                 tGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTc                                                     <  1:3036222/71‑1 (MQ=255)
                   gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                   <  1:2359212/71‑1 (MQ=255)
                   gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                   <  1:2339448/71‑1 (MQ=255)
                    aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                              >  1:767283/1‑59 (MQ=255)
                    aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                                                  >  1:93754/1‑71 (MQ=255)
                     tACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                 <  1:1743708/71‑1 (MQ=255)
                       cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                <  1:3261743/70‑1 (MQ=255)
                       cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                <  1:165066/70‑1 (MQ=255)
                       cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                <  1:4439219/70‑1 (MQ=255)
                        tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACgg                                                            <  1:115781/57‑1 (MQ=255)
                        tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                    >  1:5086101/1‑65 (MQ=255)
                        aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGTAACGTTCACTGcc                                               >  1:2488450/2‑69 (MQ=255)
                        aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGcca                                              >  1:5656746/2‑69 (MQ=255)
                          tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                                                      >  1:103855/1‑61 (MQ=255)
                          tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                             >  1:2985504/1‑70 (MQ=255)
                          tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                             >  1:5238324/1‑70 (MQ=255)
                           cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg                                            >  1:2933886/1‑70 (MQ=255)
                           cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                           >  1:3132958/1‑71 (MQ=255)
                           cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                           >  1:2927463/1‑71 (MQ=255)
                            gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                          >  1:1627117/1‑71 (MQ=255)
                               ctcggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACg                                         <  1:2015484/66‑1 (MQ=255)
                               cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                       >  1:2658389/1‑71 (MQ=255)
                               cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                       >  1:4194718/1‑71 (MQ=255)
                                 cggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGa                                           >  1:2261816/2‑65 (MQ=255)
                                  ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAGGTTCACCGCAATGACGTATAc                                    >  1:5165757/1‑71 (MQ=255)
                                  ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                         <  1:1106693/66‑1 (MQ=255)
                                   gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                   <  1:5574020/71‑1 (MQ=255)
                                   gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                   <  1:4958940/71‑1 (MQ=255)
                                    caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                <  1:342520/57‑1 (MQ=255)
                                    caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                   >  1:4368369/1‑70 (MQ=255)
                                     aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                 <  1:587474/71‑1 (MQ=255)
                                     aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                 <  1:4636744/71‑1 (MQ=255)
                                     aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                 <  1:5425847/71‑1 (MQ=255)
                                      acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa                                <  1:1357286/71‑1 (MQ=255)
                                      acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa                                <  1:3543836/71‑1 (MQ=255)
                                       caacGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                                                   >  1:2182536/1‑35 (MQ=255)
                                        aacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                               <  1:2332322/54‑1 (MQ=255)
                                           gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                           <  1:411584/71‑1 (MQ=255)
                                           gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                           <  1:3691746/71‑1 (MQ=255)
                                           gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                           <  1:348569/71‑1 (MQ=255)
                                           gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                           <  1:147057/71‑1 (MQ=255)
                                            aaTCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                           >  1:5408693/1‑70 (MQ=255)
                                             aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaaaa                         <  1:5534185/71‑1 (MQ=255)
                                             aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaaaa                         <  1:4875973/71‑1 (MQ=255)
                                              tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                                                           <  1:672586/36‑1 (MQ=255)
                                              tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc                             <  1:4885688/66‑1 (MQ=255)
                                               ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                          >  1:1740954/1‑68 (MQ=255)
                                               ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                          >  1:2412398/1‑68 (MQ=255)
                                               ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                          >  1:2656630/1‑68 (MQ=255)
                                                tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata                                     >  1:1999018/1‑56 (MQ=255)
                                                 ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTAAACCAAAAGCAAAACGtt                     >  1:476508/1‑71 (MQ=255)
                                                  tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                       <  1:3554955/68‑1 (MQ=255)
                                                    gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                  <  1:5040189/71‑1 (MQ=255)
                                                    gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTa                   <  1:312039/70‑1 (MQ=255)
                                                    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                  <  1:3155242/71‑1 (MQ=255)
                                                    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                  <  1:5059234/71‑1 (MQ=255)
                                                    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                  <  1:3045327/71‑1 (MQ=255)
                                                    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                  <  1:4165288/71‑1 (MQ=255)
                                                    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                   >  1:5248812/1‑70 (MQ=255)
                                                    gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                   <  1:5164881/70‑1 (MQ=255)
                                                     tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                 <  1:3289837/71‑1 (MQ=255)
                                                     tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                 <  1:718535/71‑1 (MQ=255)
                                                     tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                 <  1:4568437/71‑1 (MQ=255)
                                                     tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                  >  1:1653029/1‑70 (MQ=255)
                                                     tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                  >  1:2696590/1‑70 (MQ=255)
                                                     tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                  >  1:157976/1‑70 (MQ=255)
                                                     tgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaa                 <  1:5295608/71‑1 (MQ=255)
                                                         cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                              >  1:477412/1‑38 (MQ=255)
                                                         cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                              >  1:1703812/1‑38 (MQ=255)
                                                         cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                              >  1:3903751/1‑38 (MQ=255)
                                                         cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCg             <  1:2442740/71‑1 (MQ=255)
                                                         cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCg             <  1:2175992/71‑1 (MQ=255)
                                                         cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCg             <  1:1420245/71‑1 (MQ=255)
                                                          ggCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                         <  1:3771331/42‑1 (MQ=255)
                                                          ggCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGa            <  1:2237662/71‑1 (MQ=255)
                                                          ggCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGa            <  1:3707637/71‑1 (MQ=255)
                                                            cGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaa          <  1:1957919/71‑1 (MQ=255)
                                                              tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                      <  1:3591019/57‑1 (MQ=255)
                                                              tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaaa                >  1:2095174/1‑63 (MQ=255)
                                                                tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGttt                    >  1:4496144/1‑57 (MQ=255)
                                                                  tCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                       >  1:5614570/1‑36 (MQ=255)
                                                                  tCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                       >  1:1068631/1‑36 (MQ=255)
                                                                    cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc  >  1:578358/1‑71 (MQ=255)
                                                                    cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc  >  1:4718072/1‑71 (MQ=255)
                                                                    cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc   <  1:1170512/70‑1 (MQ=255)
                                                                     |                                                                     
GGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCC  >  minE/982539‑982677

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