Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 965,844 A→T 15.5% V17D (GTT→GAT)  ydeM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE965,8440AT15.5% 191.3 / 25.6 77V17D (GTT→GAT) ydeMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/24);  new base T (4/8);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 6.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGG  >  minE/965776‑965913
                                                                    |                                                                     
tCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAaca                                                                     >  1:3367171/1‑71 (MQ=255)
 cAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAtca                                                                     <  1:2407403/70‑1 (MQ=255)
 cAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAaca                                                                     >  1:425990/1‑70 (MQ=255)
        cTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTgag                                                             <  1:1215504/71‑1 (MQ=255)
        cTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTga                                                              <  1:2305268/70‑1 (MQ=255)
        cTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTga                                                              <  1:4082731/70‑1 (MQ=255)
          gTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATcaca                                                                   <  1:765196/63‑1 (MQ=255)
          gTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACAttt                                                                <  1:622198/66‑1 (MQ=255)
          gTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTg                                                               >  1:1520676/1‑67 (MQ=255)
           tCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTgaga                                                            >  1:1950539/1‑69 (MQ=255)
           tCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGAtt                                                          <  1:2128243/71‑1 (MQ=255)
           tCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGAtt                                                          <  1:2133038/71‑1 (MQ=255)
             atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACAttt                                                                >  1:932266/1‑63 (MQ=255)
             atccatccATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTTAGATTAc                                                        <  1:3873208/71‑1 (MQ=255)
                    catTTTTCATGAGTAAACTGCGACTCTTTTTGAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  >  1:321938/1‑70 (MQ=255)
                    catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaa                                                 <  1:1600048/71‑1 (MQ=255)
                    catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  <  1:3968506/70‑1 (MQ=255)
                    catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  >  1:101026/1‑70 (MQ=255)
                    catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  <  1:1865227/70‑1 (MQ=255)
                    catTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  <  1:88770/70‑1 (MQ=255)
                     atTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTa                                                         <  1:412593/62‑1 (MQ=255)
                     atTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTa                                                         <  1:1754669/62‑1 (MQ=255)
                     atTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                <  1:3925720/71‑1 (MQ=255)
                      tttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGAt                                                           >  1:1751548/1‑59 (MQ=255)
                      tttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  <  1:3006077/68‑1 (MQ=255)
                      tttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAAc                                               >  1:4084009/1‑71 (MQ=255)
                      tttttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAAc                                               <  1:1428919/71‑1 (MQ=255)
                        tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:2757771/1‑68 (MQ=255)
                        tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:2282056/1‑68 (MQ=255)
                        tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:3840430/1‑68 (MQ=255)
                        tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaa                                                 >  1:2286931/1‑67 (MQ=255)
                        tttCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTg                                             >  1:3457659/1‑71 (MQ=255)
                          tCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTga                                                              <  1:1606441/52‑1 (MQ=255)
                            aTGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaa                                                 <  1:1899852/63‑1 (MQ=255)
                              gAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGc                                       >  1:242664/1‑71 (MQ=255)
                               aGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTaaa                                                                       >  1:2854150/1‑38 (MQ=255)
                               aGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTaaa                                                                       >  1:69519/1‑38 (MQ=255)
                               aGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAaca                                                                     <  1:963551/40‑1 (MQ=255)
                               aGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCt                                      >  1:1453691/1‑71 (MQ=255)
                                 tAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTg                                    <  1:897553/71‑1 (MQ=255)
                                 tAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTg                                    >  1:3159773/1‑71 (MQ=255)
                                   aaCTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGa                                           >  1:780374/1‑62 (MQ=255)
                                   aaCTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGa                                           >  1:214944/1‑62 (MQ=255)
                                   aaCTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGc                                  >  1:53371/1‑71 (MQ=255)
                                     cTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTg                                >  1:1767555/1‑71 (MQ=255)
                                      tGCGACTCTTTTTCAAGGTAAAAACAGTAAACACAtt                                                                 >  1:1122556/1‑37 (MQ=255)
                                      tGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTg                                >  1:2617831/1‑70 (MQ=255)
                                        cGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:1577470/1‑52 (MQ=255)
                                        cGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTa                              <  1:2587788/70‑1 (MQ=255)
                                        cGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTa                              <  1:879020/70‑1 (MQ=255)
                                        cGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTa                              >  1:147067/1‑70 (MQ=255)
                                           ctctTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:1234498/1‑49 (MQ=255)
                                           ctctTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:2212611/1‑49 (MQ=255)
                                           ctctTTTTCAAGGTAAAAACAGTAAACACATTTGAGATACATTGAAAACTGGAGGGCTTGGCTGTAACGTg                         <  1:3053779/71‑1 (MQ=255)
                                             ctTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTg                         >  1:1871060/1‑70 (MQ=255)
                                             ctTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACa                                                       >  1:2574246/1‑40 (MQ=255)
                                              tttttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaa                                                  <  1:3583908/44‑1 (MQ=255)
                                               ttttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGaaaa                                                >  1:1586746/1‑45 (MQ=255)
                                                tttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATc                     >  1:3436399/1‑71 (MQ=255)
                                                 ttCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCg                    >  1:3969067/1‑71 (MQ=255)
                                                   cAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGc                   <  1:1293309/70‑1 (MQ=255)
                                                    aaGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCt                  <  1:1504517/70‑1 (MQ=255)
                                                      ggTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTg                         >  1:1380069/1‑61 (MQ=255)
                                                      ggTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtc               >  1:478600/1‑71 (MQ=255)
                                                      ggTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtc               >  1:68863/1‑71 (MQ=255)
                                                      ggTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtc               >  1:2028062/1‑71 (MQ=255)
                                                       gTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTg                         <  1:814012/60‑1 (MQ=255)
                                                       gTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtct              >  1:715040/1‑71 (MQ=255)
                                                        tAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAtc               <  1:551551/69‑1 (MQ=255)
                                                        tAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCtt                                     <  1:1775155/47‑1 (MQ=255)
                                                         aaaaaCAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATc                     <  1:2099121/62‑1 (MQ=255)
                                                          aaaaCAGTAATCACATTTGAGATTACATTGAAAACTGGa                                           >  1:28418/1‑39 (MQ=255)
                                                          aaaaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGc           >  1:995289/1‑71 (MQ=255)
                                                            aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCtt                                     >  1:1477166/1‑43 (MQ=255)
                                                            aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCt          >  1:2626487/1‑70 (MQ=255)
                                                            aaCAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCt          >  1:2115602/1‑70 (MQ=255)
                                                              cAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCaa       >  1:4014143/1‑71 (MQ=255)
                                                                 tAAACACATTTGAGATTACATTGAAAACTGGAGGGCtt                                     <  1:3443342/38‑1 (MQ=255)
                                                                 tAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCa                       >  1:2686511/1‑52 (MQ=255)
                                                                 tAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAt                >  1:1825941/1‑59 (MQ=255)
                                                                 tAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAt                >  1:1188108/1‑59 (MQ=255)
                                                                 tAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAAt      >  1:695369/1‑69 (MQ=255)
                                                                 tAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATaa    <  1:511632/71‑1 (MQ=255)
                                                                   aaCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAgg  >  1:908820/1‑71 (MQ=255)
                                                                    acacATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTAt                >  1:3874997/1‑56 (MQ=255)
                                                                    acacATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAAt      <  1:142486/66‑1 (MQ=255)
                                                                    acacATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAAt      <  1:1296279/66‑1 (MQ=255)
                                                                    |                                                                     
TCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAAACACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGG  >  minE/965776‑965913

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: