Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 93,391 G→T 5.2% G125C (GGC→TGC) ‡ secA → preprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE93,3910GT 5.2% 393.1 / 3.7 114G125C (GGC→TGC) ‡secApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base G (47/61);  new base T (2/4);  total (49/65)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 120 total aligned reads displayed.

ATGCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGA  >  minE/93322‑93461
                                                                     |                                                                      
aTGCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAgg                                                                       <  1:3380934/71‑1 (MQ=255)
aTGCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAgg                                                                       >  1:2717448/1‑71 (MQ=255)
aTGCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAgg                                                                       >  1:2119931/1‑71 (MQ=255)
 tGCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGc                                                                      <  1:2881196/71‑1 (MQ=255)
  gCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGc                                                                      <  1:1218424/70‑1 (MQ=255)
  gCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGc                                                                      <  1:2969752/70‑1 (MQ=255)
    gTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTg                                                                   >  1:1451052/1‑71 (MQ=255)
     tACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTg                                                                   >  1:2957844/1‑70 (MQ=255)
         ggTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGt                                                              <  1:1477433/71‑1 (MQ=255)
          gTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTa                                                             <  1:1175150/71‑1 (MQ=255)
            gAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           >  1:2410822/1‑71 (MQ=255)
            gAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           >  1:213817/1‑71 (MQ=255)
            gAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           >  1:40640/1‑71 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGtt                                                          <  1:3036875/71‑1 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGtt                                                          <  1:1695193/71‑1 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGtt                                                          <  1:723925/71‑1 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGtt                                                          <  1:1397037/71‑1 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           >  1:2987721/1‑70 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           >  1:2511688/1‑70 (MQ=255)
             aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           >  1:687499/1‑70 (MQ=255)
                gAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTAcc                                                       <  1:3572976/71‑1 (MQ=255)
                 aaaaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           <  1:2452216/66‑1 (MQ=255)
                   aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCggt                                                    <  1:295092/71‑3 (MQ=21)
                   aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCggt                                                    <  1:230413/71‑3 (MQ=21)
                   aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCggt                                                    <  1:172975/71‑3 (MQ=21)
                   aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCggt                                                    <  1:2775245/71‑3 (MQ=21)
                     aCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGtaac                                                  <  1:1243325/71‑4 (MQ=255)
                     aCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTc                                                    <  1:1297624/69‑1 (MQ=255)
                      cccTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAAc                                                 <  1:1650267/71‑1 (MQ=255)
                       ccTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAAc                                                 <  1:1678077/70‑1 (MQ=255)
                        cTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCg                                                      >  1:2840826/1‑64 (MQ=21)
                        cTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCg                                                      >  1:152414/1‑64 (MQ=21)
                          gACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTAc                                                        >  1:1495129/1‑60 (MQ=255)
                           aCCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGAc                                              <  1:331779/69‑1 (MQ=255)
                           aCCGCAAAGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                           <  1:2476413/56‑1 (MQ=255)
                             cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                         >  1:210448/1‑56 (MQ=255)
                             cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                         >  1:2195694/1‑56 (MQ=255)
                             cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                         >  1:2256060/1‑56 (MQ=255)
                             cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                         >  1:380900/1‑56 (MQ=255)
                             cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTAcc                                          >  1:775852/1‑71 (MQ=255)
                             cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTAcc                                          >  1:2467446/1‑71 (MQ=255)
                               cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                >  1:3231135/1‑63 (MQ=255)
                               cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                >  1:1776948/1‑63 (MQ=255)
                               cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                >  1:658442/1‑63 (MQ=255)
                               cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTg                                        <  1:3360839/71‑1 (MQ=255)
                               cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTg                                        <  1:750505/71‑1 (MQ=255)
                                aaCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCa                                                                 <  1:2409100/45‑1 (MQ=255)
                                aaCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCa                                                                 <  1:3643017/45‑1 (MQ=255)
                                aaCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCa                                                                 <  1:3217261/45‑1 (MQ=255)
                                aaCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTgg                                       >  1:1835621/1‑71 (MQ=255)
                                 aCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgc                                      <  1:1247164/71‑1 (MQ=255)
                                 aCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgc                                      <  1:1500866/71‑1 (MQ=255)
                                   gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                    <  1:2260362/71‑1 (MQ=255)
                                   gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                    <  1:92963/71‑1 (MQ=255)
                                    ctgcctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAAc                                                 >  1:1626588/1‑57 (MQ=255)
                                       cctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAAc                                 <  1:491967/70‑1 (MQ=255)
                                       cctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAAc                                 <  1:656381/70‑1 (MQ=255)
                                        ctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGt                               <  1:2346840/71‑1 (MQ=255)
                                        ctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGt                               <  1:2802575/71‑1 (MQ=255)
                                          gcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGa                             >  1:2849347/1‑71 (MQ=255)
                                           cTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGAc                            <  1:3657697/71‑1 (MQ=255)
                                           cTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGAc                            <  1:2873814/71‑1 (MQ=255)
                                           cTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGAc                            <  1:2569047/71‑1 (MQ=255)
                                              aCCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGcc                         >  1:3642335/1‑71 (MQ=255)
                                               ccTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGcc                         <  1:2003626/70‑1 (MQ=255)
                                               ccTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCg                        <  1:1585018/71‑1 (MQ=255)
                                               ccTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCg                        <  1:3475404/71‑1 (MQ=255)
                                               ccTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCg                        <  1:1798176/71‑1 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTAGCGTCAACGACTACCTg                                        >  1:2649134/1‑54 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGcc                         >  1:3709759/1‑69 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCg                        >  1:572529/1‑70 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCg                        >  1:1722504/1‑70 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                       <  1:1809848/71‑1 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                       <  1:303433/71‑1 (MQ=255)
                                                cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                       <  1:2561230/71‑1 (MQ=255)
                                                  gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                      >  1:2649424/1‑70 (MQ=255)
                                                  gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                      >  1:1089072/1‑70 (MQ=255)
                                                  gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                      >  1:752969/1‑70 (MQ=255)
                                                  gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                      >  1:3035325/1‑70 (MQ=255)
                                                   aaCGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaaaa                    <  1:3119978/71‑1 (MQ=255)
                                                   aaCGCACTAACCGGTAAAGGCGTCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaaa                     <  1:2933344/69‑1 (MQ=255)
                                                    aCGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaac                   <  1:1386094/71‑1 (MQ=255)
                                                        aCTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAAcc               <  1:3200255/71‑1 (MQ=255)
                                                            aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaac                   <  1:1433589/63‑1 (MQ=255)
                                                            aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaac                   <  1:500332/63‑1 (MQ=255)
                                                            aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTcc           >  1:93062/1‑71 (MQ=255)
                                                            aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTcc           >  1:2912696/1‑71 (MQ=255)
                                                             ccGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCg          <  1:2263508/71‑1 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaacaa                 >  1:3489152/1‑61 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       <  1:1736961/71‑1 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       <  1:2531559/71‑1 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       >  1:1625837/1‑71 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       <  1:154378/71‑1 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       >  1:3318555/1‑71 (MQ=255)
                                                                gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACATGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       >  1:577675/1‑71 (MQ=255)
                                                                   aaGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       <  1:2073850/68‑1 (MQ=255)
                                                                   aaGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg       <  1:2424753/68‑1 (MQ=255)
                                                                   aaGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGttt    <  1:1512766/71‑1 (MQ=255)
                                                                     ggCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                       <  1:2824296/50‑1 (MQ=255)
                                                                     ggCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGa  <  1:872024/71‑1 (MQ=255)
                                                                     |                                                                      
ATGCGTACCGGTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGA  >  minE/93322‑93461

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: