Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,323 G→A 5.9% F40F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,3230GA 5.9% 267.3 / 4.4 84F40F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base A (3/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGTTTTCAGTTCAT  >  minE/1934253‑1934393
                                                                      |                                                                      
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)
       cAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCg                                                                  <  1:277073/70‑1 (MQ=255)
       cAGTTCCTGGTTATCCTTCGCATCAAACTCTTCGTGCAACTGGTCAGAGGTATGGTCGCGGCGAAAATca                                                                  <  1:2108073/70‑2 (MQ=255)
        aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt                                                                >  1:1236295/1‑71 (MQ=255)
        aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt                                                                >  1:2157153/1‑71 (MQ=255)
          ttcttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCg                                                              <  1:3722214/71‑1 (MQ=255)
            cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                                             >  1:1519415/1‑70 (MQ=255)
            cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                                             >  1:1925090/1‑70 (MQ=255)
               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa                                                          <  1:2734720/70‑1 (MQ=255)
               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGaa                                                          <  1:2298148/70‑1 (MQ=255)
               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGa                                                           <  1:161816/69‑1 (MQ=255)
               cGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGa                                                           <  1:2677430/69‑1 (MQ=255)
                gTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAgg                                                        <  1:2219731/71‑1 (MQ=255)
                 ttCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGc                                                       <  1:699005/71‑1 (MQ=255)
                   ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa                                                     >  1:3139375/1‑71 (MQ=255)
                   ctctTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCaa                                                     >  1:3811096/1‑71 (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:3192608/71‑1 (MQ=255)
                                        gTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCg                                <  1:3372527/71‑1 (MQ=255)
                                          gCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGc                               >  1:3658725/1‑70 (MQ=255)
                                           cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                                             >  1:3890924/1‑39 (MQ=255)
                                           cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                                             >  1:3449737/1‑39 (MQ=255)
                                           cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                                             >  1:3060033/1‑39 (MQ=255)
                                           cAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                                                             >  1:943347/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)
                                               tGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCtt                         >  1:468530/1‑71 (MQ=255)
                                                ggTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTtc                        >  1:1695775/1‑71 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAg                                     >  1:1636447/1‑54 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAg                                     >  1:2959797/1‑54 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAg                                     >  1:371338/1‑54 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCac                    >  1:3331198/1‑71 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCac                    >  1:3586517/1‑71 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCac                    >  1:618605/1‑71 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCac                    <  1:3822475/71‑1 (MQ=255)
                                                    agagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCa                     >  1:2733781/1‑70 (MQ=255)
                                                     gagGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCac                    >  1:3076082/1‑70 (MQ=255)
                                                       ggTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGt                                   <  1:1087381/53‑1 (MQ=255)
                                                         tATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGc                           <  1:1244526/59‑1 (MQ=255)
                                                         tATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGc                           <  1:459906/59‑1 (MQ=255)
                                                          aTGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGAcgcg              >  1:3908423/1‑71 (MQ=255)
                                                             aTCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAg                            <  1:1427527/54‑1 (MQ=255)
                                                             aTCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGttt           >  1:3563643/1‑71 (MQ=255)
                                                             aTCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGtt            <  1:487502/70‑1 (MQ=255)
                                                             aTCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCGCGAGCAgg                                    >  1:2019470/1‑46 (MQ=255)
                                                              tCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGtttt          >  1:26038/1‑71 (MQ=255)
                                                               cgcgACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGtttt          >  1:70609/1‑70 (MQ=255)
                                                                gcgACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGtttt          >  1:764264/1‑69 (MQ=255)
                                                                gcgACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGtttt          >  1:229987/1‑69 (MQ=255)
                                                                gcgACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGTTTTc         <  1:647699/70‑1 (MQ=255)
                                                                  gACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCt                          >  1:2543128/1‑51 (MQ=255)
                                                                  gACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGTTTTCAGt      >  1:2031295/1‑71 (MQ=255)
                                                                    cGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGtt            >  1:3750730/1‑63 (MQ=255)
                                                                     ggAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAgg                                    >  1:3452578/1‑38 (MQ=255)
                                                                     ggAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAg                                     <  1:2169710/37‑1 (MQ=255)
                                                                     ggAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAg                                     <  1:1889372/37‑1 (MQ=255)
                                                                      gAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCa                             <  1:1358920/44‑1 (MQ=255)
                                                                      gAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGTTTTCAGTTCAt  <  1:2750403/71‑1 (MQ=255)
                                                                      gAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGTTTTCAGTTCAt  <  1:470870/71‑1 (MQ=255)
                                                                      |                                                                      
GCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGCGCAGGTTCGCCAGCTTCTCACGACGCGTTTTCAGTTCAT  >  minE/1934253‑1934393

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: