Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,568,329 A→G 5.5% L333L (TTG→CTG)  mntH ← manganese/divalent cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,568,3290AG 5.5% 312.4 / 7.1 114L333L (TTG→CTG) mntHmanganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (47/60);  new base G (3/3);  total (51/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 123 total aligned reads displayed.

AACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACC  >  minE/1568259‑1568393
                                                                      |                                                                
aaCAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCaa                                                                  >  1:1949485/1‑71 (MQ=255)
aaCAGTACCTGACTCATAACCAGAATCCGTGTCCGATCTAATCCCATCAGAATGACAATAAATGACGGCaa                                                                  >  1:2185113/1‑71 (MQ=255)
 aCAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCaa                                                                  <  1:1405311/70‑1 (MQ=255)
 aCAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAAc                                                                 <  1:2669298/71‑1 (MQ=255)
    gcaCCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACGATAAATGACGGCAGCATg                                                              >  1:3467979/3‑71 (MQ=255)
    gTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATg                                                              <  1:68188/71‑1 (MQ=255)
    gTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATg                                                              <  1:3608310/71‑1 (MQ=255)
         tGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGAc                                                         >  1:2295000/1‑71 (MQ=255)
          gACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAAGATGGTGACt                                                        >  1:4266855/1‑71 (MQ=255)
          gACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACt                                                        >  1:4631640/1‑71 (MQ=255)
           aCTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACt                                                        <  1:78638/70‑1 (MQ=255)
           aCTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAAAATGGTGACTg                                                       <  1:4068730/71‑1 (MQ=255)
            cTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTg                                                       <  1:1072022/70‑1 (MQ=255)
            cTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTg                                                       >  1:1735957/1‑70 (MQ=255)
             tCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTg                                                       <  1:4237062/69‑1 (MQ=255)
               aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACGATAAATGACGGCAGCATGGTGACTGTacg                                                   <  1:4145472/71‑1 (MQ=255)
               aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacg                                                   <  1:1972489/71‑1 (MQ=255)
               aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTac                                                    <  1:102029/70‑1 (MQ=255)
                tAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCaa                                                                  <  1:1000414/55‑1 (MQ=255)
                tAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCaa                                                                  <  1:1425607/55‑1 (MQ=255)
                tAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacga                                                  >  1:3142008/1‑71 (MQ=255)
                 aaCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacgac                                                 >  1:1268643/1‑71 (MQ=255)
                 aaCCAGAATACGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacga                                                  >  1:1817418/1‑70 (MQ=255)
                  aCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacgacg                                                <  1:488269/71‑1 (MQ=255)
                    cAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACAt                                                               >  1:1691415/1‑54 (MQ=255)
                     aGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCATCATGGTGACTGTACGACGCAc                                             >  1:3198360/1‑71 (MQ=255)
                      gAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTg                                                       <  1:4608138/60‑1 (MQ=255)
                       aaTCCGTGTCGGATCTAATCCCATCAGAATGACGATAAATGACGGCAGCATGGTGACTGTacg                                                   <  1:4179255/63‑1 (MQ=255)
                       aaTCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACa                                                                <  1:276329/50‑1 (MQ=255)
                       aaTCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTac                                                    <  1:1954519/62‑1 (MQ=255)
                       aaTCCGTGTCGGATCGAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTac                                                    <  1:2938959/62‑1 (MQ=255)
                          ccGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGa                                                          >  1:3159372/1‑53 (MQ=255)
                           cGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAg                                       >  1:1297703/1‑71 (MQ=255)
                            gtgtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGt                                                            <  1:1952108/49‑1 (MQ=255)
                            gtgtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGc                                      <  1:186277/71‑1 (MQ=255)
                             tgtCGTATCTAATCCAATCAGAATGACAATAAATGACGGCAACATg                                                              <  1:3478555/46‑1 (MQ=37)
                             tgtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCcac                                         <  1:617441/67‑1 (MQ=255)
                             tgtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCg                                     <  1:2392013/71‑1 (MQ=255)
                             tgtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCg                                     <  1:2448889/71‑1 (MQ=255)
                               tCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGt                                                            >  1:1588870/1‑46 (MQ=255)
                               tCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacga                                                  <  1:1945141/56‑1 (MQ=255)
                               tCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCggg                                   <  1:2901957/71‑1 (MQ=255)
                                cGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGa                                  >  1:2268712/1‑71 (MQ=255)
                                 ggATCTAATCCCATCAGAATGACAATAAATGACGTCAACAt                                                               >  1:3112294/1‑41 (MQ=255)
                                 ggATCTAATCCCATCAGAATGACAATAAATGACGGCAAc                                                                 >  1:1687538/1‑39 (MQ=255)
                                 ggATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGa                                  <  1:875551/70‑1 (MQ=255)
                                 ggATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGa                                  <  1:2652847/70‑1 (MQ=255)
                                 ggATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGa                                  >  1:2303363/1‑70 (MQ=255)
                                   aTCTAATCCCATCAGAATGACAATAAATGACGGCAAc                                                                 <  1:872701/37‑1 (MQ=255)
                                     cTAATCCCATCAGAATGACGATAAATGACGGCAGCATGGTGACTGTACGACGCAccc                                           >  1:3900890/1‑57 (MQ=38)
                                     cTAATCCCATCAGAATGACGATAAATGACGGCAGCATGGTGACTGTACGACGCAccc                                           >  1:2884876/1‑57 (MQ=38)
                                     cTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTa                                                     >  1:1826545/1‑47 (MQ=255)
                                     cTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCa                                              >  1:635352/1‑54 (MQ=255)
                                      tAATCCCATCAGAATGACAATAAATGACGGCAACa                                                                <  1:2170702/35‑1 (MQ=255)
                                      tAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACt                                                        <  1:633236/43‑1 (MQ=255)
                                      tAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACt                                                        <  1:56063/43‑1 (MQ=255)
                                      tAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACt                                                        <  1:4195118/43‑1 (MQ=255)
                                         tCCCATCATAATGACAATAAATGACGGCAACATGGt                                                            <  1:333423/36‑1 (MQ=255)
                                         tCCCATCAGAATGACAATAAATGACGGCAACATGGt                                                            <  1:1267697/36‑1 (MQ=255)
                                         tCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTacga                                                  >  1:436541/1‑46 (MQ=255)
                                          cccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCcac                                         >  1:4335068/1‑54 (MQ=255)
                                          cccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGc                         >  1:2181775/1‑70 (MQ=255)
                                          cccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGc                         >  1:3358239/1‑70 (MQ=255)
                                           ccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCg                        <  1:3555313/70‑1 (MQ=255)
                                           ccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCg                        <  1:2043110/70‑1 (MQ=255)
                                           ccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGa                       <  1:1738546/71‑1 (MQ=255)
                                           ccATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGa                       <  1:1083618/71‑1 (MQ=255)
                                            cATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGaa                      <  1:1290262/71‑1 (MQ=255)
                                             aTCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCgg                                    <  1:1771259/56‑1 (MQ=255)
                                             aTCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAAt                     <  1:3743097/71‑1 (MQ=255)
                                              tCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGa                             <  1:2677028/62‑1 (MQ=255)
                                              tCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGaa                      <  1:1203026/69‑1 (MQ=255)
                                              tCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAAt                     >  1:1611417/1‑70 (MQ=255)
                                              tCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATg                    <  1:3774612/71‑1 (MQ=255)
                                                aGAATGACAATAAATGACGGCAACATGGTGACTGTa                                                     >  1:1848527/1‑36 (MQ=255)
                                                aGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGaaa                           >  1:486820/1‑62 (MQ=255)
                                                aGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATg                    >  1:2835131/1‑69 (MQ=255)
                                                aGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATg                    >  1:2125275/1‑69 (MQ=255)
                                                aGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGaa                  <  1:3911814/71‑1 (MQ=255)
                                                 gAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAg                 >  1:2334537/1‑71 (MQ=255)
                                                    tGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGccc              >  1:388289/1‑71 (MQ=255)
                                                      aCAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTg            <  1:3490651/71‑1 (MQ=255)
                                                       cAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTg            >  1:2884009/1‑70 (MQ=255)
                                                        aaTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCg                        <  1:832072/57‑1 (MQ=255)
                                                        aaTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCg                        <  1:2297455/57‑1 (MQ=255)
                                                        aaTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCa          <  1:3278205/71‑1 (MQ=255)
                                                        aaTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCa          <  1:4308934/71‑1 (MQ=255)
                                                         aTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGa                   <  1:2981091/61‑1 (MQ=255)
                                                          tAAATGACGGCAACATGGTGACTTTACGACGCACCCACAGCGGGATATGAAAGCg                        <  1:1450543/55‑1 (MQ=255)
                                                          tAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCAt         >  1:3410902/1‑70 (MQ=255)
                                                           aaaTGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGaaa                           >  1:1375191/1‑51 (MQ=255)
                                                            aaTGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGaaa                           >  1:2184314/1‑50 (MQ=255)
                                                             aTGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAAt                     <  1:258015/55‑1 (MQ=255)
                                                             aTGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATcacc     <  1:2583095/71‑1 (MQ=255)
                                                              tGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGa                                  >  1:179459/1‑41 (MQ=255)
                                                              tGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATcacca    <  1:2144218/71‑1 (MQ=255)
                                                                     aaCATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCaccacc  >  1:1883033/1‑66 (MQ=255)
                                                                      aCATGGTGACTGTACGACGCACCCACAGCGGGATATg                              <  1:3593759/37‑1 (MQ=255)
                                                                      aCATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATcacca    <  1:1279169/63‑1 (MQ=255)
                                                                      aCATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATc        >  1:3726994/1‑59 (MQ=255)
                                                                      |                                                                
AACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATGAAAGCGAATGAAGCCCTGCATCACCACC  >  minE/1568259‑1568393

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: