Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,299 T→C 8.3% Q48Q (CAA→CAG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,2990TC 8.3% 232.6 / 10.1 84Q48Q (CAA→CAGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (35/42);  new base C (4/3);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCA  >  minE/1934229‑1934366
                                                                      |                                                                   
cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAt                                                                     >  1:797984/1‑71 (MQ=255)
 aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt                                                                    >  1:56651/1‑71 (MQ=255)
 aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAtt                                                                    <  1:2212748/71‑1 (MQ=255)
  gCAACGGCGACTTCGATGTTCAGTGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                                                   <  1:3438780/71‑1 (MQ=255)
  gCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                                                   <  1:994345/71‑1 (MQ=255)
     aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                                                   <  1:2617371/68‑1 (MQ=255)
     aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTc                                                                <  1:2290633/71‑1 (MQ=255)
           aCTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAgg                                                           <  1:1001321/70‑1 (MQ=255)
               cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGt                                                                 >  1:3374137/1‑60 (MQ=255)
               cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTAt                                                        >  1:2489960/1‑69 (MQ=255)
               cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATg                                                       <  1:19321/70‑1 (MQ=255)
               cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                                                      >  1:874366/1‑71 (MQ=255)
               cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                                                      >  1:2930972/1‑71 (MQ=255)
                gATGTTCAGCGCTTCCAGTTCTTCGGTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                                                      <  1:1010110/70‑1 (MQ=255)
                 aTGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCaga                                                             >  1:1201739/1‑62 (MQ=255)
                     tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATcgc                                                  >  1:177321/1‑69 (MQ=255)
                      cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                                               >  1:3419140/1‑71 (MQ=255)
                      cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                                               >  1:3699312/1‑71 (MQ=255)
                       aGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACg                                              >  1:2465454/1‑71 (MQ=255)
                        gcgcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACgg                                             <  1:249032/71‑1 (MQ=255)
                         cgcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGa                                            <  1:640171/71‑1 (MQ=255)
                          gcTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGaa                                           <  1:95399/71‑1 (MQ=255)
                          gattCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGaaa                                           <  1:3352812/69‑1 (MQ=255)
                           cTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTTGTCAGAGGTATGATCGCGACGGaa                                           <  1:1306342/70‑1 (MQ=255)
                            ttCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGaaaa                                          >  1:3717359/1‑70 (MQ=255)
                            ttCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAAt                                         >  1:1094387/1‑71 (MQ=255)
                             tCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATc                                        >  1:1628646/1‑71 (MQ=255)
                             tCCAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAAt                                         >  1:1441068/1‑70 (MQ=255)
                              ccAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAAt                                         <  1:914869/69‑1 (MQ=255)
                               cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGa                                                >  1:135864/1‑61 (MQ=255)
                               cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCg                                       <  1:277073/70‑1 (MQ=255)
                               cAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAATca                                       <  1:2108073/70‑2 (MQ=255)
                               cAGTTATTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                                                   >  1:1832902/1‑42 (MQ=255)
                                aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt                                     >  1:2157153/1‑71 (MQ=255)
                                aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt                                     >  1:1236295/1‑71 (MQ=255)
                                  ttcttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCg                                   <  1:3722214/71‑1 (MQ=255)
                                    cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa                                                               <  1:3893408/41‑1 (MQ=255)
                                    cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCa                                                               <  1:514162/41‑1 (MQ=255)
                                    cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                  >  1:1519415/1‑70 (MQ=255)
                                    cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                  >  1:1925090/1‑70 (MQ=255)
                                     ttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCagag                                                            <  1:1117770/43‑1 (MQ=255)
                                       cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa                               <  1:2298148/70‑1 (MQ=255)
                                       cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa                               <  1:2734720/70‑1 (MQ=255)
                                       cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGa                                <  1:161816/69‑1 (MQ=255)
                                       cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGa                                <  1:2677430/69‑1 (MQ=255)
                                        gTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATcgc                                                  >  1:1207895/1‑50 (MQ=255)
                                        gTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAgg                             <  1:2219731/71‑1 (MQ=255)
                                         ttCTCTTTGCCGTCGAATTCTGCGTGCATTTGGTCAGAGGTATGa                                                      <  1:2630468/45‑1 (MQ=255)
                                         ttCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                                                      <  1:3478289/45‑1 (MQ=255)
                                         ttCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGc                            <  1:699005/71‑1 (MQ=255)
                                           ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa                          >  1:3811096/1‑71 (MQ=255)
                                           ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa                          >  1:3139375/1‑71 (MQ=255)
                                           atctTTGACGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGAt                                                     >  1:2815884/2‑44 (MQ=255)
                                              tttGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa                          <  1:2659247/68‑1 (MQ=255)
                                              tttGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAAt                         <  1:1168111/69‑1 (MQ=255)
                                               ttGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCa                           >  1:1427207/1‑66 (MQ=255)
                                                tGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATcc                       <  1:2087645/69‑1 (MQ=255)
                                                tGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATcc                       <  1:1442494/69‑1 (MQ=255)
                                                      cGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcg               >  1:1531808/1‑71 (MQ=255)
                                                       gAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa                          >  1:2011953/1‑59 (MQ=255)
                                                       gAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcgc              <  1:3377989/71‑1 (MQ=255)
                                                       gAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcgc              <  1:2488411/71‑1 (MQ=255)
                                                       gAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcg               <  1:3494579/70‑1 (MQ=255)
                                                        aaTTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGa                                            <  1:364758/40‑1 (MQ=255)
                                                         aTTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcgcg             >  1:2798313/1‑70 (MQ=255)
                                                         aTTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCgcgc            <  1:2632238/71‑1 (MQ=255)
                                                         aTTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCgcgc            <  1:1695486/71‑1 (MQ=255)
                                                          ttCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa                               >  1:2455871/1‑51 (MQ=255)
                                                          ttCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcgcg             >  1:2041155/1‑69 (MQ=255)
                                                          ttCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCa           <  1:3076383/71‑1 (MQ=255)
                                                           tCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAATCAttgg                                   <  1:6512/46‑3 (MQ=16)
                                                           tCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAATCAtt                                     >  1:2034962/1‑44 (MQ=16)
                                                            cTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt                                     >  1:390649/1‑43 (MQ=255)
                                                            cTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAgg         >  1:2717760/1‑71 (MQ=255)
                                                               cGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGt        <  1:2388826/69‑1 (MQ=255)
                                                               cGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTc      <  1:3751951/71‑1 (MQ=255)
                                                               cGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTc      <  1:3018823/71‑1 (MQ=255)
                                                               cGTGCAATGGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTc      >  1:38206/1‑71 (MQ=255)
                                                                gTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCg     <  1:3372527/71‑1 (MQ=255)
                                                                gTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCg     <  1:3192608/71‑1 (MQ=255)
                                                                  gCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGc    >  1:3658725/1‑70 (MQ=255)
                                                                   cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                  >  1:3890924/1‑39 (MQ=255)
                                                                   cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                  >  1:3060033/1‑39 (MQ=255)
                                                                   cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                  >  1:943347/1‑39 (MQ=255)
                                                                   cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                  >  1:3449737/1‑39 (MQ=255)
                                                                   cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCa                           >  1:1120957/1‑46 (MQ=255)
                                                                   cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCa  <  1:3539421/71‑1 (MQ=255)
                                                                      ttGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATc                        <  1:488857/46‑1 (MQ=255)
                                                                      |                                                                   
CAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCA  >  minE/1934229‑1934366

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: