Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,687 G→A 16.2% F252F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6870GA16.2% 189.1 / 26.6 80F252F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (38/29);  new base A (7/6);  total (45/35)
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.86e-01

TGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAG  >  minE/1933617‑1933755
                                                                      |                                                                    
tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                      <  1:3220104/71‑1 (MQ=255)
tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                      <  1:2740201/71‑1 (MQ=255)
tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                      <  1:249991/71‑1 (MQ=255)
  aaCTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                      >  1:3250418/1‑69 (MQ=255)
   aCTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCaaaa                                                                   <  1:3645128/71‑2 (MQ=255)
    ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                  >  1:3504748/1‑71 (MQ=255)
    ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                  <  1:843694/71‑1 (MQ=255)
    ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                  >  1:3670150/1‑71 (MQ=255)
      ctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                                                                >  1:1724942/1‑71 (MQ=255)
       tGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGccac                                                               >  1:3473895/1‑71 (MQ=255)
           ttATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCgcctac                                                           <  1:922462/71‑1 (MQ=255)
           ttATGACGTACGGAAATACCTTCGTTAGGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCaccaac                                                           <  1:1063126/71‑1 (MQ=255)
              tGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                        <  1:1300805/71‑1 (MQ=255)
               gACGTACGGATATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                                    >  1:1279412/1‑58 (MQ=255)
               gACGTACGGATATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                                    >  1:1274816/1‑58 (MQ=255)
                aCGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccacc                                                       <  1:3086059/70‑1 (MQ=255)
                 cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                        >  1:3630647/1‑68 (MQ=255)
                 cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                        >  1:1376230/1‑68 (MQ=255)
                  gTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGcca                                                                >  1:189450/1‑59 (MQ=255)
                    aCGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACg                                                  >  1:3649797/1‑71 (MQ=255)
                     cagAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAg                                                     >  1:629623/3‑67 (MQ=255)
                     cagAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACgt                                                 >  1:3141714/3‑70 (MQ=255)
                     cGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGc                                                 >  1:2166419/1‑71 (MQ=255)
                      ggAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCt                                                >  1:1701368/1‑71 (MQ=255)
                      ggAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCt                                                >  1:2088152/1‑71 (MQ=255)
                      agAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACgtt                                                <  1:3260411/70‑3 (MQ=255)
                      agAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACgtt                                                <  1:1161142/70‑3 (MQ=255)
                       gAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCtt                                               >  1:2154592/1‑71 (MQ=255)
                        aaaTACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGAc                                                   <  1:3478955/66‑1 (MQ=255)
                              cTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTAc                                        <  1:2315807/71‑1 (MQ=255)
                               ttCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCaaaac                                                                  >  1:2638097/1‑42 (MQ=25)
                               ttCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATa                                         >  1:3514976/1‑69 (MQ=255)
                                tcattACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAg                                      >  1:3712641/4‑71 (MQ=255)
                                 cGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGt                                     >  1:2433687/1‑71 (MQ=255)
                                  attACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTAcaacca                                                      <  1:3163864/52‑1 (MQ=255)
                                   ttACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCt                                                <  1:3123190/58‑1 (MQ=255)
                                   ttACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTc                                   <  1:1908235/71‑1 (MQ=255)
                                    tACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGa                                                    >  1:1414325/1‑53 (MQ=255)
                                     aCGGAAGTTACGGTTGATTTCGAATACCCGTTCaaaa                                                                   >  1:2408405/1‑36 (MQ=25)
                                     aCGGAAGTTACGGTTGATTTCGAATACCCGTTCaaaa                                                                   >  1:984208/1‑36 (MQ=25)
                                     aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTcc                                  <  1:495593/70‑1 (MQ=255)
                                       ggAAGTTACGGTTGATTTCGAATACACTCTCGAAGc                                                                  <  1:2014571/36‑1 (MQ=255)
                                       ggAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                  <  1:3189086/36‑1 (MQ=255)
                                       ggAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                  >  1:855054/1‑36 (MQ=255)
                                       ggAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgc                               >  1:623995/1‑71 (MQ=255)
                                        gAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcg                              <  1:792838/71‑1 (MQ=255)
                                         aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                        <  1:2790031/44‑1 (MQ=255)
                                         aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCTGACGCTTGAGGTACAGTTCCGgcgc                             <  1:2946860/71‑1 (MQ=255)
                                         aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                             >  1:2645146/1‑71 (MQ=255)
                                         aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                             >  1:3912495/1‑71 (MQ=255)
                                          aGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                             >  1:1843283/1‑70 (MQ=255)
                                           gTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                        <  1:1967341/42‑1 (MQ=255)
                                           gTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACAAGACGCTTGAGGTACAGTTCCGGCGCGa                           >  1:1364136/1‑71 (MQ=255)
                                              aCGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATAc                        >  1:3718751/1‑71 (MQ=255)
                                               cGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATAc                        <  1:1133437/70‑1 (MQ=255)
                                               cGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                       <  1:616068/71‑1 (MQ=255)
                                                ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                       >  1:1168678/1‑70 (MQ=255)
                                                 gTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCa                     <  1:1078146/71‑1 (MQ=255)
                                                  ttGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                      <  1:3474315/69‑1 (MQ=255)
                                                  ttGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAg                    >  1:3582482/1‑71 (MQ=255)
                                                  ttGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAg                    >  1:3080684/1‑71 (MQ=255)
                                                    gATTTCGAATACACGCTCGAAGCCACCAACCACCAGa                                                    >  1:2640767/1‑37 (MQ=255)
                                                    gATTTCGAATACACGCTAGAAGCCACCAACCACCAGa                                                    >  1:1152665/1‑37 (MQ=255)
                                                     aTTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                 >  1:3416192/1‑71 (MQ=255)
                                                     aTTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                 <  1:2983683/71‑1 (MQ=255)
                                                      tttCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATa                 >  1:911576/1‑70 (MQ=255)
                                                       ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTcc                                  >  1:1181936/1‑52 (MQ=255)
                                                       ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa               >  1:1290634/1‑71 (MQ=255)
                                                       ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa               >  1:1286058/1‑71 (MQ=255)
                                                       ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa               >  1:3801662/1‑71 (MQ=255)
                                                         cGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACAt              <  1:3335801/70‑1 (MQ=255)
                                                         cGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACAt              <  1:638973/70‑1 (MQ=255)
                                                         cGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATg             >  1:1181427/1‑71 (MQ=255)
                                                          gAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGt            >  1:587896/1‑71 (MQ=255)
                                                          gAATACACGCTCGAAGCCACCAACCACCAGACGCTTAGGTACAGTTCCGgcgc                             <  1:339026/53‑1 (MQ=255)
                                                           aaTACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGt            <  1:1217053/70‑1 (MQ=255)
                                                             tACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATcta         <  1:2333770/71‑3 (MQ=255)
                                                             tACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATct          >  1:1983765/1‑69 (MQ=255)
                                                             tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAg                                      <  1:3768193/42‑1 (MQ=255)
                                                             tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTc           >  1:1214196/1‑69 (MQ=255)
                                                             tACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCga         <  1:86483/71‑1 (MQ=255)
                                                              acacGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                 <  1:2422356/62‑1 (MQ=255)
                                                               cacGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCgaga       <  1:3669173/71‑1 (MQ=255)
                                                                  gCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATcc    <  1:1093881/71‑1 (MQ=255)
                                                                  gCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATc     >  1:893998/1‑70 (MQ=255)
                                                                    tCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCg                                 <  1:2364790/40‑1 (MQ=255)
                                                                    tCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAg  >  1:571229/1‑71 (MQ=255)
                                                                     cGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCgaga       >  1:126340/1‑65 (MQ=255)
                                                                      |                                                                    
TGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAG  >  minE/1933617‑1933755

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: