Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 288,729 A→T 5.6% E458V (GAA→GTA)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE288,7290AT 5.6% 362.5 / 8.6 122E458V (GAA→GTA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (59/56);  new base T (3/4);  total (62/60)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 123 total aligned reads displayed.

CCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTG  >  minE/288662‑288799
                                                                   |                                                                      
ccGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGt                                                                     >  1:1809984/1‑71 (MQ=255)
 cGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGAAAGt                                                                     <  1:1799857/70‑1 (MQ=255)
     ctctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTTCCAGGTAGTggt                                                                  <  1:263352/69‑2 (MQ=255)
     ctctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTTCCAGGTAGTggt                                                                  <  1:1656018/69‑2 (MQ=255)
       ctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTTCCAGGTAGTggtcg                                                                >  1:715264/1‑66 (MQ=38)
       ctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTTCCAGGTAGTggtcg                                                                >  1:2268091/1‑66 (MQ=38)
       ctGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTTCCAGGTAGTggtcg                                                                >  1:270446/1‑66 (MQ=38)
          cACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACg                                                             <  1:435618/69‑1 (MQ=255)
            cTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGAtct                                                         <  1:248960/71‑1 (MQ=255)
            cTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGAtct                                                         <  1:1980611/71‑1 (MQ=255)
             tGCTTGAAGTGCTGGATCCAGAGCAGAATGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGAtct                                                         >  1:162841/1‑70 (MQ=255)
                    aGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACg                                                 <  1:957546/71‑1 (MQ=255)
                    aGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACg                                                 <  1:1011376/71‑1 (MQ=255)
                    aGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACg                                                 <  1:1489812/71‑1 (MQ=255)
                     gTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGt                                                <  1:309700/71‑1 (MQ=255)
                     gTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGt                                                <  1:2090850/71‑1 (MQ=255)
                      tGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGt                                                >  1:1570036/1‑70 (MQ=255)
                        cTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGAt                                             >  1:3728884/1‑71 (MQ=255)
                        cTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGAt                                             >  1:3605419/1‑71 (MQ=255)
                        cTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGAt                                             >  1:1897685/1‑71 (MQ=255)
                           gATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                          >  1:3097959/1‑71 (MQ=255)
                           gATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                          >  1:3434322/1‑71 (MQ=255)
                           gATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                          >  1:3148970/1‑71 (MQ=255)
                           gATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                          >  1:3496089/1‑71 (MQ=255)
                           gATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                          >  1:1714986/1‑71 (MQ=255)
                           gATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGt                                           <  1:2970246/70‑1 (MQ=255)
                            aTCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGt                                           <  1:236178/69‑1 (MQ=255)
                             tCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg                                        <  1:2600358/71‑1 (MQ=255)
                               cAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                       <  1:2559638/70‑1 (MQ=255)
                               cAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTc                                      >  1:2131936/1‑71 (MQ=255)
                                agagCAGAACGTAGCGTTCAGCGACCACTTCCAGGTAGTggt                                                                  <  1:1331399/42‑2 (MQ=25)
                                 gagCAGAACGTAGCGTTCAGCGACCACTACCTGGTAGTGGATTACGATCTCAGCGACGTGATGt                                           <  1:2941984/64‑1 (MQ=255)
                                 gagCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGt                                           <  1:3530521/64‑1 (MQ=255)
                                 gagCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGt                                           <  1:314766/64‑1 (MQ=255)
                                 gagCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                       >  1:2389418/1‑68 (MQ=255)
                                 gagCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGt                                       >  1:1470581/1‑68 (MQ=255)
                                   gCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGttt                                         >  1:3307799/1‑64 (MQ=255)
                                   gCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGa                                  >  1:3489867/1‑71 (MQ=255)
                                    cAGAACTTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTg                                               >  1:7383/1‑57 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTg                                               >  1:2533156/1‑57 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTg                                               >  1:2529059/1‑57 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTg                                               >  1:3339815/1‑57 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTg                                               >  1:225053/1‑57 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                 <  1:837256/71‑1 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                 >  1:3687133/1‑71 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                 <  1:1684620/71‑1 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                 <  1:198546/71‑1 (MQ=255)
                                    cAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGAc                                 <  1:1985673/71‑1 (MQ=255)
                                     aGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACg                                <  1:2455930/71‑1 (MQ=255)
                                      gAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                               >  1:2481253/1‑71 (MQ=255)
                                      gAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                               <  1:3639557/71‑1 (MQ=255)
                                      gAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                               >  1:3573963/1‑71 (MQ=255)
                                      gAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                               >  1:1689513/1‑71 (MQ=255)
                                      gAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGt                               <  1:3255065/71‑1 (MQ=255)
                                        aCGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCg                             <  1:3568150/71‑1 (MQ=255)
                                             gCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                        <  1:488267/71‑1 (MQ=255)
                                              cGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                        >  1:1087323/1‑70 (MQ=255)
                                              cGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTc                        >  1:2371012/1‑70 (MQ=255)
                                               gTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTcc                       <  1:2157796/70‑1 (MQ=255)
                                                ttCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACg                                <  1:2739716/60‑1 (MQ=255)
                                                ttCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACg                                <  1:462747/60‑1 (MQ=255)
                                                ttCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCAt                     <  1:2308193/71‑1 (MQ=255)
                                                 tCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTcc                       >  1:865482/1‑68 (MQ=255)
                                                 tCAGCGACCACTACCTGGAAGTGGATTACGATCTAAGCGACGTGATGTTTGTCGCGACGTCGAACTcc                       >  1:633998/1‑68 (MQ=255)
                                                  cAGCTACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                            <  1:2161602/46‑1 (MQ=255)
                                                  cAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                            <  1:2703085/46‑1 (MQ=255)
                                                  cAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                            <  1:1529337/46‑1 (MQ=255)
                                                  cAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATg                    >  1:2321913/1‑70 (MQ=255)
                                                  cAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATg                    <  1:3428254/70‑1 (MQ=255)
                                                   aGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGa                                  >  1:234160/1‑55 (MQ=255)
                                                   aGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGa                                  >  1:1383863/1‑55 (MQ=255)
                                                    gCGACCACTACCTGTAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAAc                          <  1:145716/62‑1 (MQ=255)
                                                    gCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg                                        <  1:2388214/48‑1 (MQ=255)
                                                    gCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg                                        <  1:2914004/48‑1 (MQ=255)
                                                    gCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAAc                          <  1:2783438/62‑1 (MQ=255)
                                                    gCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAAc                          <  1:754069/62‑1 (MQ=255)
                                                    gCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGaa                  >  1:2264817/1‑70 (MQ=255)
                                                     cGACCACTACCTGGAAGTGGATTACGATCTCAGCGa                                                   >  1:517040/1‑36 (MQ=255)
                                                      gACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAt               >  1:754249/1‑71 (MQ=255)
                                                      gACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAt               >  1:3648463/1‑71 (MQ=255)
                                                      gACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAt               <  1:2490446/71‑1 (MQ=255)
                                                      gACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAt               >  1:1285453/1‑71 (MQ=255)
                                                       aCCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGtt                                          <  1:1970709/43‑1 (MQ=255)
                                                       aCCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACAtt              >  1:862534/1‑71 (MQ=255)
                                                        ccACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTc             <  1:125309/71‑1 (MQ=255)
                                                        ccACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTc             >  1:225168/1‑71 (MQ=255)
                                                         cACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTg                                        >  1:2873051/1‑43 (MQ=255)
                                                         cACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCg                             <  1:2503184/54‑1 (MQ=255)
                                                         cACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCg                             <  1:3254447/54‑1 (MQ=255)
                                                         cACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTcc            <  1:2726428/71‑1 (MQ=255)
                                                          aCTACCTGGAAGTGGATTACGATCTCAGCGACGTGATg                                            <  1:1563380/38‑1 (MQ=255)
                                                            tACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCgg          >  1:116113/1‑70 (MQ=255)
                                                            tACCTGGAAGTGGATTACGATCTCAGCGAAGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCgg          >  1:421123/1‑70 (MQ=255)
                                                              ccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcg                                     >  1:526434/1‑41 (MQ=255)
                                                              ccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcg                                     >  1:1374110/1‑41 (MQ=255)
                                                              ccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTcg                                     >  1:2303938/1‑41 (MQ=255)
                                                              ccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCa        <  1:2407534/70‑1 (MQ=255)
                                                              ccTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCa        <  1:2510541/70‑1 (MQ=255)
                                                                 ggAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgc    >  1:427033/1‑71 (MQ=255)
                                                                  gAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCgct   <  1:1312125/71‑1 (MQ=255)
                                                                   |                                                                      
CCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTG  >  minE/288662‑288799

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: