Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,749,136 G→C 8.3% G403G (GGG→GGC secY → preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,749,1360GC 8.3% 440.8 / 22.1 144G403G (GGG→GGCsecYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (77/55);  new base C (4/8);  total (81/63)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 151 total aligned reads displayed.

TATTACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTC  >  minE/2749067‑2749204
                                                                     |                                                                    
tatTACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTggg                                                                      <  1:4826057/70‑1 (MQ=255)
tatTACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGa                                                                     >  1:44200/1‑71 (MQ=255)
tatTACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGa                                                                     >  1:2101766/1‑71 (MQ=255)
  ttACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGAcc                                                                   >  1:1059333/1‑71 (MQ=255)
  ttACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGAcc                                                                   >  1:1819393/1‑71 (MQ=255)
     ccTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTc                                                                 <  1:1199163/70‑1 (MQ=255)
       tttATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTggg                                                                      <  1:2652480/63‑1 (MQ=255)
       tttATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTggg                                                                      <  1:701733/63‑1 (MQ=255)
       tttATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTggg                                                                      <  1:1283854/63‑1 (MQ=255)
         tATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACt                                                              <  1:4937699/69‑1 (MQ=255)
         tATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACt                                                              <  1:318387/69‑1 (MQ=255)
         tATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGGTCTACTTCGGTGGGACCTCACt                                                              <  1:3764691/69‑1 (MQ=255)
                cTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCgt                                                     >  1:3277515/1‑71 (MQ=255)
                      ccGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCAc                                                               >  1:2685803/1‑55 (MQ=255)
                      ccGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGttgttg                                                >  1:2680075/1‑70 (MQ=255)
                      ccGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGttgttg                                                <  1:3959419/70‑1 (MQ=255)
                       cGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGttgttg                                                <  1:2493194/69‑1 (MQ=255)
                              cATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:2856254/69‑1 (MQ=255)
                              cATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:3242683/69‑1 (MQ=255)
                              cATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:2001491/69‑1 (MQ=255)
                              cATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:4131207/69‑1 (MQ=255)
                              cATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:4420140/69‑1 (MQ=255)
                              cATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:1333971/69‑1 (MQ=255)
                                 gCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:4588896/66‑1 (MQ=255)
                                   gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                  <  1:623377/71‑1 (MQ=255)
                                   gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                  <  1:3372880/71‑1 (MQ=255)
                                   gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                  <  1:709671/71‑1 (MQ=255)
                                   gTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGAc                                  <  1:840211/71‑1 (MQ=255)
                                      aTGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATc                                                       >  1:4639277/1‑47 (MQ=255)
                                      aTGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                >  1:925416/1‑70 (MQ=255)
                                      aTGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                >  1:3521659/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                   >  1:4153319/1‑65 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                   >  1:4875228/1‑65 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGa                                   >  1:990965/1‑65 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                              >  1:289450/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                              >  1:1573622/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                              >  1:5157511/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                              >  1:4707745/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                              >  1:3750444/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTa                              >  1:2485480/1‑70 (MQ=255)
                                        gCAATGAAAGTACCGTTCTACTTAGGTGGGACCTCACTGCt                                                           >  1:1170159/1‑41 (MQ=255)
                                         cAATGAAAGTACCGTTCTACTTCGGTGGGACCTCa                                                                >  1:2065539/1‑35 (MQ=255)
                                         cAATGAAAGTACCGTTCTACTTCGGTGGGACCTCa                                                                >  1:1239595/1‑35 (MQ=255)
                                         cAATGAAAGTACCGTTCTACTTCGGTGGGACCTCa                                                                >  1:3070460/1‑35 (MQ=255)
                                         cAATGAAAGTACCGTTCTACTTCGGTGGGACCTCa                                                                >  1:15770/1‑35 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGc                          >  1:1920914/1‑70 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                         >  1:3650084/1‑71 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                         >  1:608828/1‑71 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                         >  1:4630681/1‑71 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                         >  1:3991493/1‑71 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                         >  1:2069657/1‑71 (MQ=255)
                                            tGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCt                         >  1:1740008/1‑71 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:5085184/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:1434510/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:68063/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:1243909/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:2824878/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:1042127/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCtt                                                          >  1:422229/1‑36 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTTATTATGGACTTTATGGCTc                        <  1:565287/70‑1 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATTGACTTTATGGCTc                        <  1:853024/70‑1 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTc                        <  1:4756991/70‑1 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTc                        <  1:3901876/70‑1 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTc                        <  1:2028481/70‑1 (MQ=255)
                                              aaaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCa                       <  1:5173154/71‑1 (MQ=255)
                                               aaGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCaa                      <  1:491119/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)
                                                  tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:3883139/49‑1 (MQ=255)
                                                  tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGAt                                         <  1:2579221/49‑1 (MQ=255)
                                                  tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                <  1:3174251/58‑1 (MQ=255)
                                                  tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                <  1:4515746/58‑1 (MQ=255)
                                                  tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                <  1:2672345/58‑1 (MQ=255)
                                                  tACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                <  1:4667503/58‑1 (MQ=255)
                                                   aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTAt                                      >  1:1341593/1‑51 (MQ=255)
                                                   aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATgg                           >  1:2359240/1‑62 (MQ=255)
                                                   aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATg                            >  1:3558737/1‑61 (MQ=255)
                                                   aCCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATg                            >  1:2863725/1‑61 (MQ=255)
                                                    ccGTTCTACTTCGGTGGGTCCTCACTGCTTATCGTTGTTGTCGTGATTATGGACtt                                <  1:3157542/56‑1 (MQ=255)
                                                    ccGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCa                 <  1:2954252/71‑1 (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:5156547/2‑69 (MQ=255)
                                                      cttCTACTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTGCa                 >  1:4182873/2‑69 (MQ=255)
                                                       ttCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAAc              <  1:904888/71‑1 (MQ=255)
                                                         cTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAActc            >  1:5167/1‑71 (MQ=255)
                                                           aCTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAActct           >  1:3354370/1‑70 (MQ=255)
                                                           aCTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAActct           >  1:275706/1‑70 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTTATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                   <  1:4761051/62‑1 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                   <  1:1390712/62‑1 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                   <  1:2450033/62‑1 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                   <  1:3645042/62‑1 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                   <  1:1605746/62‑1 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAAGTg                   <  1:1696585/62‑1 (MQ=255)
                                                           aCTTCGGTGGCACGTCGCTCCTGATCATCGTGGTGGTGACGATGGACTTCATGGCACAACTg                   <  1:2228080/62‑1 (MQ=255)
                                                                   gggACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTc  >  1:3145537/1‑71 (MQ=255)
                                                                   gggACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTc  >  1:2124711/1‑71 (MQ=255)
                                                                     |                                                                    
TATTACCTTTATCTGCCTGATCCCGGAGTTCATGCGTGATGCAATGAAAGTACCGTTCTACTTCGGTGGGACCTCACTGCTTATCGTTGTTGTCGTGATTATGGACTTTATGGCTCAAGTGCAAACTCTGATGATGTC  >  minE/2749067‑2749204

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