Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,279,253 G→C 15.3% V344V (GTG→GTC yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,279,2530GC15.3% 208.9 / 26.6 85V344V (GTG→GTCyeeJadhesin
Reads supporting (aligned to +/- strand):  ref base G (44/28);  new base C (9/4);  total (53/32)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

AGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCC  >  minE/1279185‑1279322
                                                                    |                                                                     
aGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGc                                                                     >  1:215755/1‑71 (MQ=255)
 gCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGcc                                                                    >  1:4672137/1‑71 (MQ=255)
   tGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCt                                                                  >  1:5298725/1‑71 (MQ=255)
       tACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtt                                                               <  1:5211454/70‑1 (MQ=255)
       tACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:789441/1‑71 (MQ=255)
        aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:276759/1‑70 (MQ=255)
        aCCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:3480601/1‑70 (MQ=255)
          ccGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttg                                                             >  1:930567/1‑67 (MQ=255)
           cGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtttgg                                                            >  1:582751/1‑66 (MQ=255)
           cGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGAt                                                           >  1:399332/1‑70 (MQ=255)
            gCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaa                                                         <  1:4646674/71‑1 (MQ=255)
             ccTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATaaa                                                        <  1:3471370/71‑1 (MQ=255)
                ggCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGAc                                                     >  1:622087/1‑71 (MQ=255)
                 gCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaacat                                                     <  1:2447591/70‑4 (MQ=255)
                 gCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                      <  1:5044865/69‑1 (MQ=255)
                  ccGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   <  1:212857/71‑1 (MQ=255)
                    gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc                                                 >  1:5469915/1‑71 (MQ=255)
                    gCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATc                                                 <  1:1670750/71‑1 (MQ=255)
                      aCCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg                                               >  1:3774442/1‑71 (MQ=255)
                       ccTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCgg                                               >  1:1661023/1‑70 (MQ=255)
                        cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGccc                                                                   >  1:2790648/1‑49 (MQ=255)
                        cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                                                          <  1:4725793/58‑1 (MQ=255)
                        cTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATa                                                          <  1:5462014/58‑1 (MQ=255)
                            gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg                                          <  1:5163512/70‑1 (MQ=255)
                            gCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAg                                          <  1:4415364/70‑1 (MQ=255)
                             cGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa                                        <  1:4158708/71‑1 (MQ=255)
                               gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                                                                <  1:74015/45‑1 (MQ=255)
                               gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGt                                                                <  1:3053862/45‑1 (MQ=255)
                               gTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAAt                                      >  1:4577140/1‑71 (MQ=255)
                                tAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTg                                                                 <  1:4469042/43‑1 (MQ=255)
                                 aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTc                                                              >  1:4641011/1‑45 (MQ=255)
                                 aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:3612501/1‑64 (MQ=255)
                                 aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc                                    >  1:4000099/1‑71 (MQ=255)
                                 aaaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATcc                                    >  1:3678873/1‑71 (MQ=255)
                                  aaCTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGtt                                                               >  1:1547429/1‑43 (MQ=255)
                                      ggTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATg                               >  1:5237960/1‑71 (MQ=255)
                                       gTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc                              <  1:713950/71‑1 (MQ=255)
                                        tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCa                                             <  1:5048629/55‑1 (MQ=255)
                                        tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCa                                             <  1:4801790/55‑1 (MQ=255)
                                         cTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGc                              >  1:2481538/1‑69 (MQ=255)
                                          tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   <  1:4473728/47‑1 (MQ=255)
                                          tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   <  1:1009565/47‑1 (MQ=255)
                                          tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTa                                        <  1:1886378/58‑1 (MQ=255)
                                            tGAACAGTATTATGGCTATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:947842/1‑40 (MQ=16)
                                            tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGc                                              >  1:132388/1‑50 (MQ=255)
                                            tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATaa                         >  1:1902824/1‑71 (MQ=255)
                                            tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATa                          <  1:4411527/70‑1 (MQ=255)
                                            tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:644916/1‑40 (MQ=16)
                                            tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:5562060/1‑40 (MQ=16)
                                            tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:4522852/1‑40 (MQ=16)
                                            tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTaaa                                                        >  1:247319/1‑40 (MQ=16)
                                            tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                      >  1:547210/1‑42 (MQ=19)
                                            tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGa                                                      >  1:3504138/1‑42 (MQ=19)
                                             gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACgct                                                  >  1:1865255/1‑43 (MQ=255)
                                             gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:4982896/1‑52 (MQ=255)
                                             gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:5495788/1‑52 (MQ=255)
                                             gAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCaaa                                           >  1:1861234/1‑52 (MQ=255)
                                              aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGa                                                   >  1:1065708/1‑43 (MQ=255)
                                              aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGAt                                                  >  1:1696726/1‑44 (MQ=255)
                                              aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGAt                                                  >  1:390284/1‑44 (MQ=255)
                                              aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAAcc                       >  1:3649832/1‑71 (MQ=255)
                                               aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCg                      <  1:2886515/71‑1 (MQ=255)
                                                  gTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCa                                 <  1:3520479/57‑1 (MQ=255)
                                                    attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCg                                                >  1:3319100/1‑40 (MQ=255)
                                                    attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCg                                                >  1:1046219/1‑40 (MQ=255)
                                                     ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCa                                 >  1:1897872/1‑54 (MQ=255)
                                                     ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACGGCTgg                  <  1:2634890/69‑1 (MQ=255)
                                                     ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTgg                  <  1:690727/69‑1 (MQ=255)
                                                     ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTgg                  <  1:3714117/69‑1 (MQ=255)
                                                     ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGAc                >  1:3419483/1‑71 (MQ=255)
                                                     ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGAc                <  1:409899/71‑1 (MQ=255)
                                                      tatGGCCATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATaa                         >  1:4468066/1‑61 (MQ=255)
                                                          gCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc           <  1:795407/71‑1 (MQ=255)
                                                           cGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAAc           >  1:3659153/1‑70 (MQ=255)
                                                            gATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACta         >  1:260368/1‑71 (MQ=255)
                                                             aTGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTaa                                       <  1:4151459/40‑1 (MQ=255)
                                                              tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACtata       <  1:1159844/71‑1 (MQ=255)
                                                               gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                            >  1:1669836/1‑49 (MQ=255)
                                                               gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                            >  1:5288081/1‑49 (MQ=255)
                                                               gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCa                            >  1:773287/1‑49 (MQ=255)
                                                               gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAc      >  1:4278948/1‑71 (MQ=255)
                                                                  gTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAcccc   >  1:3522912/1‑71 (MQ=255)
                                                                   tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTc                                  >  1:1573677/1‑39 (MQ=255)
                                                                   tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTg                   <  1:1941665/54‑1 (MQ=255)
                                                                   tGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATAccccc  >  1:5304496/1‑71 (MQ=255)
                                                                    |                                                                     
AGCTGGCTACCCGCCTGGCCGCACCTTGGCGGTAAACTGGTCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATACCCCC  >  minE/1279185‑1279322

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: