Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,284 G→A 6.3% F53F (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,2840GA 6.3% 245.7 / 4.6 79F53F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base A (3/2);  total (40/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.43e-01

GGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGC  >  minE/1934214‑1934354
                                                                      |                                                                      
gggTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                                                                        <  1:3623730/71‑1 (MQ=255)
  gTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                                      <  1:1904250/71‑1 (MQ=255)
   tcatcaTGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                                     >  1:2079312/1‑71 (MQ=255)
    catcatGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                                     >  1:2449958/1‑70 (MQ=255)
       catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCt                                                                  <  1:3895473/70‑1 (MQ=255)
       catGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTg                                                                 >  1:1004369/1‑71 (MQ=255)
        atGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGc                                                                >  1:3580720/1‑71 (MQ=255)
             gCCAGCAACGGCGACTTCTATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa                                                           <  1:3632756/71‑1 (MQ=255)
             gCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                                      >  1:430752/1‑60 (MQ=255)
              ccAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                                     <  1:590094/60‑1 (MQ=255)
               cAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAAt                                                         >  1:797984/1‑71 (MQ=255)
                aGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                                    >  1:989762/1‑59 (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)
                  cAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                                      >  1:1360918/1‑55 (MQ=255)
                    aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCa                                                           <  1:3577816/64‑1 (MQ=255)
                    aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTg                                                       <  1:2617371/68‑1 (MQ=255)
                    aCGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTc                                                    <  1:2290633/71‑1 (MQ=255)
                          aCTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAgg                                               <  1:1001321/70‑1 (MQ=255)
                           cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                                              >  1:1893917/1‑54 (MQ=255)
                           cTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGt                                                              >  1:3077641/1‑54 (MQ=255)
                              cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAt                                                                     <  1:1798553/44‑1 (MQ=255)
                              cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTg                                                             >  1:1663257/1‑52 (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:2930972/1‑71 (MQ=255)
                              cGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                                          >  1:874366/1‑71 (MQ=255)
                               gATGTTCAGCGCTTCCAGTTCTTCGGTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGa                                          <  1:1010110/70‑1 (MQ=255)
                                aTGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCaga                                                 >  1:1201739/1‑62 (MQ=255)
                                 tGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGaa                                                                      >  1:3409677/1‑40 (MQ=255)
                                   ttCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTc                                                                   <  1:2001492/41‑1 (MQ=255)
                                    tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCg                                                                        <  1:810349/35‑1 (MQ=255)
                                    tCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATcgc                                      >  1:177321/1‑69 (MQ=255)
                                     cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                                    >  1:2138374/1‑38 (MQ=255)
                                     cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                                   >  1:3419140/1‑71 (MQ=255)
                                     cAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGAc                                   >  1:3699312/1‑71 (MQ=255)
                                     cAGAGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAAtt                                                                    >  1:1711027/1‑38 (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:1236295/1‑71 (MQ=255)
                                               aGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGtt                         >  1:2157153/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:1925090/1‑70 (MQ=255)
                                                   cttcGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCgg                      >  1:1519415/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:2488411/71‑1 (MQ=255)
                                                                      gAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcgc  <  1:3377989/71‑1 (MQ=255)
                                                                      gAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTcg   <  1:3494579/70‑1 (MQ=255)
                                                                      |                                                                      
GGGTCATCATGCGGCCAGCAACGGCGACTTCGATGTTCAGCGCTTCCAGTTCTTCGTTCTCTTTGCCGTCGAATTCTGCGTGCAATTGGTCAGAGGTATGATCGCGACGGAAATCGTTCGGGAAGGCAATCCCCTGCTCGC  >  minE/1934214‑1934354

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: