Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,749,166 T→G 6.8% I413M (ATT→ATG) ‡ secY → preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,749,1660TG 6.8% 474.8 / 23.7 174I413M (ATT→ATG) ‡secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (80/82);  new base G (4/8);  total (84/90)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

GTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTGAAGAAGGC  >  minE/2749102‑2749234
                                                                |                                                                    
gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                                                <  1:840211/71‑1 (MQ=255)
gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                                                <  1:709671/71‑1 (MQ=255)
gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                                                <  1:623377/71‑1 (MQ=255)
   aTGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                                              >  1:1186307/1‑70 (MQ=255)
   aTGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                                              >  1:3521659/1‑70 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                                                 >  1:990965/1‑65 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                                                 >  1:4875228/1‑65 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                                                 >  1:4153319/1‑65 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                                                 >  1:2684845/1‑65 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:2485480/1‑70 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:4707745/1‑70 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:3750444/1‑70 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:5157511/1‑70 (MQ=255)
     gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:289450/1‑70 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGc                                                        >  1:1920914/1‑70 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                                                       >  1:3650084/1‑71 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                                                       >  1:3991493/1‑71 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                                                       >  1:4630681/1‑71 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                                                       >  1:2069657/1‑71 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                                                       >  1:608828/1‑71 (MQ=255)
         tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTAATGGc                                                        >  1:4399015/1‑70 (MQ=255)
          gAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                                                       >  1:3858092/1‑70 (MQ=255)
           aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTTATTATGGACTTTATGGCTc                                                      <  1:565287/70‑1 (MQ=255)
           aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTc                                                      <  1:3901876/70‑1 (MQ=255)
           aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:5173154/71‑1 (MQ=255)
            aaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCaa                                                    <  1:2597586/71‑1 (MQ=255)
            aaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCaa                                                    <  1:1023758/71‑1 (MQ=255)
            aaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:330823/70‑1 (MQ=255)
              gTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGt                                                  >  1:4843923/1‑71 (MQ=255)
               tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                                              <  1:4667503/58‑1 (MQ=255)
               tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                                              <  1:4515746/58‑1 (MQ=255)
               tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                                              <  1:2179014/58‑1 (MQ=255)
               tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                                              <  1:736380/58‑1 (MQ=255)
                aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTAt                                                                    >  1:1341593/1‑51 (MQ=255)
                aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATgg                                                         >  1:2036085/1‑62 (MQ=255)
                aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATgg                                                         >  1:3445066/1‑62 (MQ=255)
                aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATgg                                                         >  1:3108003/1‑62 (MQ=255)
                   gTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCaaa                                             <  1:1620798/71‑1 (MQ=255)
                   gTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCaaa                                             <  1:2649739/71‑1 (MQ=255)
                   cttCTACTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTGCag                                              >  1:2278205/2‑69 (MQ=255)
                   cttCTACTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTGCag                                              >  1:3950087/2‑69 (MQ=255)
                    ttCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAAc                                            <  1:904888/71‑1 (MQ=255)
                    ttCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAAc                                            >  1:1599686/1‑71 (MQ=255)
                      cTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAActc                                          >  1:5167/1‑71 (MQ=255)
                        aCTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAActct                                         >  1:3354370/1‑70 (MQ=255)
                        aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                                                 <  1:2450033/62‑1 (MQ=255)
                        aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                                                 <  1:3645042/62‑1 (MQ=255)
                            cGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCtgatg                                     <  1:1951655/70‑1 (MQ=255)
                                gggACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTc                                >  1:3145537/1‑71 (MQ=255)
                                gggACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTc                                >  1:1814316/1‑71 (MQ=255)
                                    ccTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCcagt                            <  1:853096/71‑1 (MQ=255)
                                    ccTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCcagt                            <  1:2738942/71‑1 (MQ=255)
                                    ccTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCcagt                            <  1:3253437/71‑1 (MQ=255)
                                      tCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:2952470/44‑1 (MQ=255)
                                      tCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:115012/44‑1 (MQ=255)
                                      tCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:3005357/44‑1 (MQ=255)
                                      tCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:4349559/44‑1 (MQ=255)
                                      tCACTGCTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                                                     <  1:3909292/43‑1 (MQ=39)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:4758041/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:1017749/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:960014/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:306375/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:4957291/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:2374483/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:2566546/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:3351776/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:4089249/1‑35 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGGTCAAGTGCAAACTCTGATGATGTCcagtcag                         <  1:4401192/70‑1 (MQ=255)
                                        aCTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCcagtcagt                        <  1:2225115/71‑1 (MQ=255)
                                        aCTGCCTATCGTTGTTGTCGTGATTATGGACTTTa                                                            >  1:1651664/1‑35 (MQ=255)
                                           gCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCtga                                       >  1:3227004/1‑53 (MQ=255)
                                           gCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCtga                                       >  1:232260/1‑53 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCTGTATGa                    <  1:665052/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:2784193/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:147051/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:1699125/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:5186758/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:4917811/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:1496856/71‑1 (MQ=255)
                                            cTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGa                    <  1:2711893/71‑1 (MQ=255)
                                              tATCGTTGTTGTCTTGATTATGGACTTTATGGCTCat                                                    <  1:3166591/37‑2 (MQ=39)
                                              tATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGt                  >  1:3575880/1‑71 (MQ=255)
                                                tCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCt                >  1:3704028/1‑71 (MQ=255)
                                                tCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCt                >  1:3197322/1‑71 (MQ=255)
                                                     gttgtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCAt            <  1:5052636/70‑1 (MQ=255)
                                                      ttgtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          <  1:4122230/71‑1 (MQ=255)
                                                      ttgtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          <  1:4781405/71‑1 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTga         >  1:1689285/1‑70 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          >  1:1419621/1‑69 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          >  1:3824860/1‑69 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          >  1:3528480/1‑69 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          >  1:1190040/1‑69 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          >  1:4142024/1‑69 (MQ=255)
                                                        gtCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTg          >  1:4276127/1‑69 (MQ=255)
                                                         tCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTgaag       <  1:5018209/71‑1 (MQ=255)
                                                          cGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTgaag       <  1:3293914/70‑1 (MQ=255)
                                                          cGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTgaag       <  1:3412243/70‑1 (MQ=255)
                                                           gTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTgaagaa     <  1:2669898/71‑1 (MQ=255)
                                                           gTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTgaag       >  1:583523/1‑69 (MQ=255)
                                                           gTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTgaag       >  1:887597/1‑69 (MQ=255)
                                                                |                                                                    
GTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTCCAGTCAGTATGAGTCTGCATTGAAGAAGGC  >  minE/2749102‑2749234

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: