Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 93,594 G→A 39.9% A192A (GCG→GCA 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,5940GA39.9% 65.2 / 89.8 130A192A (GCG→GCAsecApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base G (33/45);  new base A (17/35);  total (50/80)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 132 total aligned reads displayed.

AGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCT  >  minE/93531‑93662
                                                               |                                                                    
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:1859370/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:214950/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:2172728/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:2203809/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:2775026/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:2712713/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:1569820/69‑1 (MQ=255)
aGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAg                                                                 <  1:2208213/69‑1 (MQ=255)
       aTCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCTCGACAACATGGCGTTCAGCCCTg                                                            <  1:3417384/67‑1 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:513830/1‑69 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:721967/1‑69 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:2198900/1‑69 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:1438957/1‑69 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:591538/1‑69 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:2600116/1‑69 (MQ=255)
          aCTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTgaagaa                                                       >  1:2958027/1‑69 (MQ=255)
            ttACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACg                                                     >  1:534901/1‑69 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                   <  1:1429139/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                   <  1:598004/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                   <  1:1404814/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:3013392/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:2175322/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:3157239/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:3171165/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:2043344/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:3350462/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:971524/67‑1 (MQ=255)
                ggTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:1291285/67‑1 (MQ=255)
                 gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                   <  1:1893465/66‑1 (MQ=255)
                 gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                   <  1:1162826/66‑1 (MQ=255)
                 gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACgtg                                                   <  1:1043112/66‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          <  1:255712/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTTCAACGTa                                          <  1:1605897/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1826658/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:49065/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:3225193/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1239876/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2935515/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2825793/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2640282/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:567088/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:683776/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:832255/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1482419/69‑1 (MQ=255)
                       acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:880210/69‑1 (MQ=255)
                                aCGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:526566/60‑1 (MQ=255)
                                  ggCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTg                                                            <  1:120100/40‑1 (MQ=255)
                                       tGACTAGCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:439655/1‑69 (MQ=255)
                                       tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATgcgc                           >  1:2148118/1‑68 (MQ=255)
                                       tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:3310844/1‑69 (MQ=255)
                                       tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:1200442/1‑69 (MQ=255)
                                       tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:3418405/1‑69 (MQ=255)
                                       tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:296164/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACATCGTAAACt                                      >  1:3200227/1‑56 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:112571/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:947887/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:3098195/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1476632/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2145714/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2865794/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:280082/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2295761/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1309555/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2699484/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2394200/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1624388/1‑69 (MQ=255)
                                        gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2437881/1‑69 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2547755/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2240424/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:1432293/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:922983/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2394425/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2025417/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:3100502/1‑48 (MQ=255)
                                            aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:3093506/1‑48 (MQ=255)
                                             ccTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                   <  1:1408823/38‑1 (MQ=255)
                                                gcgcGAGAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1809960/69‑1 (MQ=255)
                                                gcgcGACATCATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2168117/69‑1 (MQ=255)
                                                gcgcGACAAGATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2925674/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:740159/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2460020/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:422515/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1098980/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1776150/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1996787/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2359720/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2368702/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2368703/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:3137878/69‑1 (MQ=255)
                                                gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:433236/69‑1 (MQ=255)
                                                    gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:1915813/69‑1 (MQ=255)
                                                    gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:757058/69‑1 (MQ=255)
                                                    gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:1614071/69‑1 (MQ=255)
                                                    gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:1080815/69‑1 (MQ=255)
                                                             gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:2245357/1‑69 (MQ=255)
                                                             gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:1811316/1‑69 (MQ=255)
                                                             gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:420196/1‑69 (MQ=255)
                                                             gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:1411652/1‑69 (MQ=255)
                                                               gTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCt  >  1:368903/1‑69 (MQ=255)
                                                               gTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCt  >  1:962602/1‑69 (MQ=255)
                                                               |                                                                    
AGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCT  >  minE/93531‑93662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: