Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,708,913 G→A 22.2% E54E (GAG→GAA yhfZ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,708,9130GA22.2% 207.6 / 54.2 99E54E (GAG→GAAyhfZconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/43);  new base A (12/10);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ATAGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATATACGAA  >  minE/2708845‑2708981
                                                                    |                                                                    
aTAGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGt                                                                     >  1:146153/1‑70 (MQ=255)
 tAGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGt                                                                     >  1:2950888/1‑69 (MQ=255)
 tAGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTc                                                                    >  1:4109908/1‑70 (MQ=255)
 tAGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCa                                                                   <  1:1860611/71‑1 (MQ=255)
       aGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCgg                                                              <  1:205267/70‑1 (MQ=255)
        gCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAg                                                                 <  1:105933/66‑1 (MQ=255)
         cTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGc                                                           <  1:2570202/71‑1 (MQ=255)
            gCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCg                                                          >  1:3155246/1‑69 (MQ=255)
            gCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGa                                                         >  1:1515502/1‑70 (MQ=255)
              aaaTGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATAc                                                      >  1:2706802/1‑71 (MQ=255)
               aatgaatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACg                                                     <  1:475255/71‑1 (MQ=255)
               aatgaatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACg                                                     <  1:2739228/71‑1 (MQ=255)
                 tgaatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGa                                                   >  1:3283824/1‑71 (MQ=255)
                  gaatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGAAACGGAt                                                  >  1:966581/1‑71 (MQ=255)
                   aatgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCaa                                                                  >  1:2382129/1‑54 (MQ=255)
                    atgCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTg                                                <  1:3878530/71‑1 (MQ=255)
                         gTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCgg                                                              <  1:3833675/52‑1 (MQ=255)
                          ttCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGcc                                          <  1:1176151/71‑1 (MQ=255)
                           tCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGcc                                          >  1:250537/1‑70 (MQ=255)
                            cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGGGGAGCGATACGGATTGAACGCCGt                                        >  1:1756292/1‑71 (MQ=255)
                            cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGt                                        >  1:3724457/1‑71 (MQ=255)
                            cATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGa                                                   >  1:1375940/1‑60 (MQ=255)
                             aTCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGAt                                                        <  1:1284677/54‑1 (MQ=255)
                             aTCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGAt                                                        <  1:788844/54‑1 (MQ=255)
                              tCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGAtt                                                 <  1:1509512/60‑1 (MQ=255)
                               ccGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTg                                                <  1:2214552/60‑1 (MQ=255)
                                cGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTgg                                      >  1:2672779/1‑69 (MQ=255)
                                cGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgc                                    >  1:31357/1‑71 (MQ=255)
                                cGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgc                                    >  1:2903755/1‑71 (MQ=255)
                                cGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgc                                    >  1:3725484/1‑71 (MQ=255)
                                cGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGaa                                              >  1:2909924/1‑61 (MQ=255)
                                  ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                                  <  1:2040263/71‑1 (MQ=255)
                                  ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                                  <  1:3154134/71‑1 (MQ=255)
                                  ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                                  <  1:3677519/71‑1 (MQ=255)
                                  ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                                  <  1:2259726/71‑1 (MQ=255)
                                  ttGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGGGGgcgc                                  <  1:36450/71‑1 (MQ=255)
                                    ggCCTGACGCAGGCGGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCaa                                >  1:1504774/1‑71 (MQ=255)
                                    ggCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTc                                                                    >  1:557081/1‑35 (MQ=255)
                                    ggCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGa                                                   >  1:1943040/1‑52 (MQ=255)
                                    ggCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCaa                                >  1:1253655/1‑71 (MQ=255)
                                     gCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACg                                                     <  1:3026154/49‑1 (MQ=255)
                                     gCCTGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAAt                               <  1:3407387/71‑1 (MQ=255)
                                      ccTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATg                              <  1:2462936/71‑1 (MQ=255)
                                       cTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGAtt                                                 >  1:923096/1‑51 (MQ=255)
                                        tGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATgg                             >  1:1916551/1‑70 (MQ=255)
                                        tGACGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATgg                             >  1:2629406/1‑70 (MQ=255)
                                         gACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATgg                             >  1:4048491/1‑69 (MQ=255)
                                          aCGCTGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                                  <  1:1392994/63‑1 (MQ=255)
                                          aCGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCa                           <  1:444106/70‑1 (MQ=255)
                                          aCGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAg                          <  1:1094105/71‑1 (MQ=255)
                                          aCGCAGGCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCa                           <  1:3040757/70‑1 (MQ=255)
                                           cGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATg                              >  1:1904813/1‑66 (MQ=255)
                                           cGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                         >  1:1601438/1‑71 (MQ=255)
                                           cGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                         <  1:273320/71‑1 (MQ=255)
                                            gCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGt                         <  1:3598109/70‑1 (MQ=255)
                                              aGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTAt                      <  1:2904745/71‑1 (MQ=255)
                                               ggCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATc                     <  1:809831/71‑1 (MQ=255)
                                               ggCCGCGTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATc                     <  1:3741221/71‑1 (MQ=255)
                                                gCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATc                     >  1:888524/1‑70 (MQ=255)
                                                 ccGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCaa                                >  1:1421141/1‑58 (MQ=255)
                                                 ccGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTg                   >  1:1558413/1‑71 (MQ=255)
                                                  cgcgTTGAAAACGCTGGAGTCAGGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTgg                  >  1:3666824/1‑71 (MQ=255)
                                                  cgcgTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTAt                      >  1:386572/1‑67 (MQ=255)
                                                   gcgTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGt                 >  1:2201822/1‑71 (MQ=255)
                                                   gcgTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATc                     <  1:2943590/67‑1 (MQ=255)
                                                     gTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCg               <  1:2222932/71‑1 (MQ=255)
                                                     gTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTc                >  1:1990978/1‑70 (MQ=255)
                                                     gTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTc                >  1:437269/1‑70 (MQ=255)
                                                     gTTGAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTc                >  1:36353/1‑70 (MQ=255)
                                                       tGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGtt                        <  1:2294338/60‑1 (MQ=255)
                                                        gAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGgcgc                                  <  1:3785943/49‑1 (MQ=255)
                                                        gAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATc                     >  1:2404091/1‑62 (MQ=255)
                                                        gAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg             >  1:4018181/1‑70 (MQ=255)
                                                        gAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg             >  1:3002208/1‑70 (MQ=255)
                                                        gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTgg                  >  1:3916998/1‑65 (MQ=255)
                                                        gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg             <  1:2620358/70‑1 (MQ=255)
                                                        gAAAACGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAg             <  1:621377/70‑1 (MQ=255)
                                                        gAAAACGATGGAGTCAAGCGGAGCGATACGGATTGAACg                                            >  1:190533/1‑39 (MQ=255)
                                                         aaaaCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGa           <  1:1524085/71‑1 (MQ=255)
                                                         aaaaCGCTCGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGa           <  1:1311578/71‑1 (MQ=255)
                                                          aaaCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTAt                      >  1:1104349/1‑59 (MQ=255)
                                                          aaaCGCTGGAATCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAgat           <  1:637390/70‑3 (MQ=255)
                                                           aaCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGaa          >  1:4142287/1‑70 (MQ=255)
                                                           aaCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAat         <  1:3416209/71‑1 (MQ=255)
                                                            aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata        <  1:2944265/71‑1 (MQ=255)
                                                            aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata        <  1:628937/71‑1 (MQ=255)
                                                            aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata        <  1:1681288/71‑1 (MQ=255)
                                                            aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata        <  1:2505037/71‑1 (MQ=255)
                                                            aCGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAata        <  1:1073749/71‑1 (MQ=255)
                                                             cGCTGGAGTCAAGCGTAGCGATACGGATTGAACACCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAatat       <  1:3900053/71‑1 (MQ=255)
                                                             cGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAatat       <  1:3762785/71‑1 (MQ=255)
                                                             cGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAatat       <  1:455362/71‑1 (MQ=255)
                                                             cGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAatat       <  1:3644871/71‑1 (MQ=255)
                                                             cGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAatat       <  1:3087029/71‑1 (MQ=255)
                                                              gCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTgg                                      <  1:169084/39‑1 (MQ=255)
                                                                tGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATATAc     >  1:62294/1‑70 (MQ=255)
                                                                  gAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATATACg    >  1:2372440/1‑69 (MQ=255)
                                                                  gAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATATACGaa  <  1:1140455/71‑1 (MQ=255)
                                                                  gAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATATACGaa  <  1:289195/71‑1 (MQ=255)
                                                                    |                                                                    
ATAGATGAGCTGGCAAATGAATGCCGTTCATCCGTTGGCCTGACGCAGGCCGCGTTGAAAACGCTGGAGTCAAGCGGAGCGATACGGATTGAACGCCGTGGGCGCAATGGCAGTTATCTGGTCGAGGAATATACGAA  >  minE/2708845‑2708981

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: