Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 93,618 G→A 43.5% Q200Q (CAG→CAA 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,6180GA43.5% 25.2 / 66.0 88Q200Q (CAG→CAAsecApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base A (16/22);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCT  >  minE/93554‑93662
                                                                |                                            
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          <  1:255712/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTTCAACGTa                                          <  1:1605897/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1673298/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2386515/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2640282/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2825793/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2120677/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:2935515/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1826658/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:3225193/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:3345140/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:43884/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1482419/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:832255/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:880210/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1218897/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1239876/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:1393753/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:683776/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:49065/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:567088/69‑1 (MQ=255)
acaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:525530/69‑1 (MQ=255)
         aCGGCTTTGACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTa                                          <  1:526566/60‑1 (MQ=255)
                tGACTAGCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:439655/1‑69 (MQ=255)
                tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATgcgc                           >  1:2148118/1‑68 (MQ=255)
                tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:3418405/1‑69 (MQ=255)
                tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:3310844/1‑69 (MQ=255)
                tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:296164/1‑69 (MQ=255)
                tGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                          >  1:1200442/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACATCGTAAACt                                      >  1:3200227/1‑56 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:996995/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2394200/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:947887/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2699484/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:280082/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2865794/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:3098195/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2437881/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2295761/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1476632/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1570643/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1624388/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:1309555/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2237874/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:2145714/1‑69 (MQ=255)
                 gACTACCTGCGCGACAACATGGCATTCAGCCCTGAAGAACGTGTACAACGTAAACTGCACTATGCGCtg                         >  1:112571/1‑69 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2025417/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:1432293/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:1450388/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:692623/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:3100502/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:1944935/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:3093506/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2193532/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2394425/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2240424/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:2547755/1‑48 (MQ=255)
                     aCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTa                                          >  1:922983/1‑48 (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:2368702/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2368703/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2359720/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1426968/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:433236/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:422515/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1098980/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:740159/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1776150/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2460020/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:3137878/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2604446/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1996787/69‑1 (MQ=255)
                         gcgcGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:2850128/69‑1 (MQ=255)
                             gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:757058/69‑1 (MQ=255)
                             gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:1080815/69‑1 (MQ=255)
                             gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:1614071/69‑1 (MQ=255)
                             gACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg             <  1:1915813/69‑1 (MQ=255)
                                   aTGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTcc       <  1:244342/69‑1 (MQ=255)
                                      gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:1411652/1‑69 (MQ=255)
                                      gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:420196/1‑69 (MQ=255)
                                      gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:1811316/1‑69 (MQ=255)
                                      gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:2156995/1‑69 (MQ=255)
                                      gCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc    >  1:2245357/1‑69 (MQ=255)
                                        gTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCt  >  1:368903/1‑69 (MQ=255)
                                        gTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCt  >  1:962602/1‑69 (MQ=255)
                                           cAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGa                 <  1:1003142/51‑1 (MQ=255)
                                                                |                                            
ACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCT  >  minE/93554‑93662

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: