Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,016,338 G→A 18.9% D6N (GAT→AAT)  yghB → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,016,3380GA18.9% 223.5 / 41.1 95D6N (GAT→AAT) yghBconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base A (13/5);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

ATTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCGTTGTTT  >  minE/2016268‑2016407
                                                                      |                                                                     
aTTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAg                                                                       <  1:813050/71‑1 (MQ=255)
  tCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                                      >  1:2906379/1‑70 (MQ=255)
        gTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcatc                                                               <  1:3797466/71‑1 (MQ=255)
           cAATAGGCAGATAAAGGCTTAAACGCTGTTCCATAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCt                                                            >  1:3686200/1‑71 (MQ=255)
           cAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGatat                                                                   >  1:2020821/1‑64 (MQ=255)
           cAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCt                                                            >  1:1416762/1‑71 (MQ=255)
             aTAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTgc                                                          >  1:3609402/1‑71 (MQ=255)
              tAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcg                                                         >  1:1143948/1‑71 (MQ=255)
              tAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcg                                                         >  1:1290784/1‑71 (MQ=255)
                ggCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGct                                                       <  1:3141166/71‑1 (MQ=255)
                ggCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGct                                                       <  1:2769212/71‑1 (MQ=255)
                 gCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcgc                                                        <  1:333681/69‑1 (MQ=255)
                   aGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCt                                                            >  1:3054020/1‑63 (MQ=255)
                    gATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                                      <  1:2435904/52‑1 (MQ=255)
                    gATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                                             <  1:2236188/61‑1 (MQ=255)
                     aTAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGGTCTGGc                                                  <  1:3500601/71‑1 (MQ=255)
                     aTAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                                  <  1:2970574/71‑1 (MQ=255)
                      tAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                                             >  1:1783815/1‑59 (MQ=255)
                      tAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                                             >  1:3869457/1‑59 (MQ=255)
                      tAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                                  >  1:2666488/1‑70 (MQ=255)
                      tAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGc                                                  >  1:2326041/1‑70 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGctct                                                     >  1:771736/1‑66 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGctct                                                     >  1:999354/1‑66 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:952287/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:1986867/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:518059/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:1884048/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:3816433/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                >  1:3081285/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                                >  1:1897425/1‑71 (MQ=255)
                       aaaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                                >  1:3699759/1‑71 (MQ=255)
                        aaGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCaa                                                <  1:189128/70‑1 (MQ=255)
                         aGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCaa                                                >  1:2865322/1‑69 (MQ=255)
                           gCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                                             >  1:1406251/1‑70 (MQ=255)
                            cTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                                  <  1:1252639/48‑1 (MQ=255)
                            cTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                                  <  1:2879832/48‑1 (MQ=255)
                            cTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGa                                           <  1:1206360/71‑1 (MQ=255)
                             ttAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGAc                                          >  1:3597890/1‑71 (MQ=255)
                             ttAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGAc                                          >  1:325231/1‑71 (MQ=255)
                              tAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcat                                                                >  1:2029687/1‑48 (MQ=255)
                              tAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGc                                                             >  1:2869517/1‑51 (MQ=255)
                               aaaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACtt                                        >  1:1001323/1‑71 (MQ=255)
                               aaaCGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACt                                         <  1:3787397/70‑1 (MQ=255)
                                  cGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAata                                                                    >  1:2123492/1‑40 (MQ=255)
                                   gCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                                  <  1:2348940/57‑1 (MQ=255)
                                    cTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcg                                   >  1:2302415/1‑71 (MQ=255)
                                       ttCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcg                                                         <  1:2332726/46‑1 (MQ=255)
                                        tCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               >  1:3501824/1‑71 (MQ=255)
                                        tCCACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGc                                                  <  1:2513742/52‑1 (MQ=255)
                                        tCCACACGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               >  1:278209/1‑71 (MQ=255)
                                         ccACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGctc                                                      <  1:3689393/47‑1 (MQ=255)
                                         ccACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGctc                                                      <  1:3404273/47‑1 (MQ=255)
                                         ccACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               <  1:4087863/70‑1 (MQ=255)
                                         ccACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               <  1:2299391/70‑1 (MQ=255)
                                         ccACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt                               <  1:621472/70‑1 (MQ=255)
                                         ccACAGGAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTg                              >  1:3932027/1‑71 (MQ=255)
                                          cacaGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg                             >  1:1331817/1‑71 (MQ=255)
                                             aGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCTCGCTGGCgg                          <  1:1491682/71‑1 (MQ=255)
                                              ggAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                                             <  1:1337952/51‑1 (MQ=255)
                                              ggAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                                >  1:3046285/1‑64 (MQ=255)
                                              ggAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa                         <  1:4079030/71‑1 (MQ=255)
                                              ggAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa                         <  1:286294/71‑1 (MQ=255)
                                               gAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGAt                        >  1:1958510/1‑71 (MQ=255)
                                               gAAAGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                                  <  1:1280413/61‑1 (MQ=255)
                                                 aaGTCCATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACACGACtt                                        >  1:954421/1‑53 (MQ=255)
                                                  aGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTg                                      <  1:2293775/54‑1 (MQ=255)
                                                    tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgcg                                 <  1:1153182/57‑1 (MQ=255)
                                                    tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgcg                                 <  1:1539263/57‑1 (MQ=255)
                                                    tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg                             <  1:518132/61‑1 (MQ=255)
                                                    tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg                             <  1:1818467/61‑1 (MQ=255)
                                                    tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa                   >  1:3425914/1‑71 (MQ=255)
                                                    tCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa                   <  1:1672462/71‑1 (MQ=255)
                                                    tCCATGGCTGTTATTCAAAATATCATCGGTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa                   <  1:1892940/71‑1 (MQ=255)
                                                     ccATGGCTGTTATTCAAGATATCATCGCTGCGCTCTgg                                                   >  1:3327133/1‑38 (MQ=255)
                                                     ccATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCa                                                 >  1:2222646/1‑40 (MQ=255)
                                                      cATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCata                 <  1:3121209/71‑1 (MQ=255)
                                                      cATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCata                 <  1:1530913/71‑1 (MQ=255)
                                                      cATGGCTGTTATTCAAAATATCATCGCTGCGCTCTGGCAACGCGACTTTGCCGCGCTGGCGGATCCTCat                  >  1:4042837/1‑70 (MQ=255)
                                                       aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                           <  1:3560448/44‑1 (MQ=255)
                                                       aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                           <  1:2280085/44‑1 (MQ=255)
                                                       aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                           <  1:1335919/44‑1 (MQ=255)
                                                       aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGAc                                          <  1:640691/45‑1 (MQ=255)
                                                       aTGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGAc                                          <  1:1840982/45‑1 (MQ=255)
                                                        tGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATAtt               <  1:332073/71‑1 (MQ=255)
                                                          gCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                                           >  1:521532/1‑41 (MQ=255)
                                                          gCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGt             >  1:2928817/1‑71 (MQ=255)
                                                            tGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACttt                                       >  1:2868032/1‑43 (MQ=255)
                                                            tGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACtt                                        >  1:240032/1‑42 (MQ=255)
                                                               tATTCAAAATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCg        >  1:1437301/1‑71 (MQ=255)
                                                                aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGc                                     >  1:998323/1‑41 (MQ=255)
                                                                aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCata                 >  1:2473215/1‑61 (MQ=255)
                                                                aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATAtt               >  1:864845/1‑63 (MQ=255)
                                                                aTTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgt       >  1:500452/1‑71 (MQ=255)
                                                                    aaGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCgttgt    <  1:3871331/70‑1 (MQ=255)
                                                                     aGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTTATATTGTTAGCGTTGttt  <  1:3943647/71‑1 (MQ=255)
                                                                     aGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCGTTGttt  <  1:1833176/71‑1 (MQ=255)
                                                                      |                                                                     
ATTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCATATTGTTAGCGTTGTTT  >  minE/2016268‑2016407

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: