Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 93,402 A→C 5.2% V128V (GTA→GTC) ‡ 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,4020AC 5.2% 327.7 / 11.1 115V128V (GTA→GTC) ‡secApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base A (51/58);  new base C (2/4);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 116 total aligned reads displayed.

GTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGG  >  minE/93332‑93470
                                                                      |                                                                    
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:1397037/71‑1 (MQ=255)
   aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGtt                                                                   <  1:1695193/71‑1 (MQ=255)
   aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                                    >  1:2511688/1‑70 (MQ=255)
   aaGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                                    >  1:2987721/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)
       aaaaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCg                                                               >  1:1454036/1‑71 (MQ=255)
         aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCggt                                                             <  1:172975/71‑3 (MQ=21)
         aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAACTACGTGCACGCCTTTACCggt                                                             <  1:2775245/71‑3 (MQ=21)
         aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGt                                                              >  1:3172550/1‑70 (MQ=255)
         aaaCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGAGCACGTAGTTACCGt                                                              >  1:2473023/1‑70 (MQ=255)
           aCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGtaac                                                           <  1:1243325/71‑4 (MQ=255)
           aCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTc                                                             <  1:2044694/69‑1 (MQ=255)
           aCCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTc                                                             <  1:1297624/69‑1 (MQ=255)
            cccTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAAc                                                          <  1:1650267/71‑1 (MQ=255)
             ccTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAAc                                                          <  1:1678077/70‑1 (MQ=255)
                gACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTAc                                                                 >  1:2972092/1‑60 (MQ=255)
                gACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTAc                                                                 >  1:1495129/1‑60 (MQ=255)
                 aCCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGAc                                                       <  1:331779/69‑1 (MQ=255)
                 aCCGCAAAGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGt                                                                    <  1:2476413/56‑1 (MQ=255)
                   cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                                  >  1:2195694/1‑56 (MQ=255)
                   cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                                  >  1:2808460/1‑56 (MQ=255)
                   cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                                  >  1:2256060/1‑56 (MQ=255)
                   cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                                  >  1:380900/1‑56 (MQ=255)
                   cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTa                                                                  >  1:618937/1‑56 (MQ=255)
                   cGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTAcc                                                   >  1:775852/1‑71 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                         >  1:658442/1‑63 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                         >  1:750950/1‑63 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                         >  1:1776948/1‑63 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACg                                                         >  1:3231135/1‑63 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTg                                                 <  1:3360839/71‑1 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTg                                                 <  1:750505/71‑1 (MQ=255)
                     cAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTg                                                 <  1:2371858/71‑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:92963/71‑1 (MQ=255)
                         gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                             <  1:2260362/71‑1 (MQ=255)
                          ctgcctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAAc                                                          >  1:1626588/1‑57 (MQ=255)
                             cctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAAc                                          <  1:656381/70‑1 (MQ=255)
                             cctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAAc                                          <  1:491967/70‑1 (MQ=255)
                              ctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACTTAGTTACCGTCAACGACTACCTGGCGCAACGt                                        <  1:2222797/71‑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:2569047/71‑1 (MQ=255)
                                 cTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGAc                                     <  1:2873814/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:2561230/71‑1 (MQ=255)
                                      cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                                <  1:3640165/71‑1 (MQ=255)
                                      cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                                <  1:303433/71‑1 (MQ=255)
                                      cTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGa                                <  1:1809848/71‑1 (MQ=255)
                                        gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                               >  1:752969/1‑70 (MQ=255)
                                        gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                               >  1:1089072/1‑70 (MQ=255)
                                        gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                               >  1:3035325/1‑70 (MQ=255)
                                        gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                               >  1:2649424/1‑70 (MQ=255)
                                         aaCGCACTAACCGGTAAAGGCGTCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaaa                              <  1:2933344/69‑1 (MQ=255)
                                           cGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaac                            >  1:1541816/1‑70 (MQ=255)
                                              aCTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAAcc                        <  1:3200255/71‑1 (MQ=255)
                                                 aaCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTc                     >  1:39067/1‑71 (MQ=255)
                                                  aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaac                            <  1:1433589/63‑1 (MQ=255)
                                                  aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaac                            <  1:500332/63‑1 (MQ=255)
                                                  aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTcc                    >  1:2912696/1‑71 (MQ=255)
                                                   ccGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCg                   <  1:2263508/71‑1 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaacaa                          >  1:3489152/1‑61 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAaacaa                          >  1:290578/1‑61 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                >  1:3318555/1‑71 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                <  1:1736961/71‑1 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                <  1:154378/71‑1 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                >  1:1625837/1‑71 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                <  1:2531559/71‑1 (MQ=255)
                                                      gTAAAGGCGTGCACGTAGTTACCGTCAACGACTACATGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                >  1:577675/1‑71 (MQ=255)
                                                         aaGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                <  1:2424753/68‑1 (MQ=255)
                                                         aaGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTg                <  1:2073850/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)
                                                              gTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAAtt        >  1:3304905/1‑71 (MQ=255)
                                                                 cACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCt     <  1:2915237/71‑1 (MQ=255)
                                                                   cGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTg   <  1:2342479/71‑1 (MQ=255)
                                                                    gTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTc                     <  1:1298759/52‑1 (MQ=255)
                                                                    gTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTgg  >  1:2797227/1‑71 (MQ=255)
                                                                    gTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTgg  >  1:1750245/1‑71 (MQ=255)
                                                                    gTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTgg  >  1:349657/1‑71 (MQ=255)
                                                                     tAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaaaa                             >  1:3092527/1‑43 (MQ=255)
                                                                     tAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaaaa                             >  1:122264/1‑43 (MQ=255)
                                                                      |                                                                    
GTGAAGGAAAAACCCTGACCGCAACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGG  >  minE/93332‑93470

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: