Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,681 A→C 23.6% R254R (CGT→CGG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6810AC23.6% 112.2 / 57.6 85R254R (CGT→CGGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (37/28);  new base C (10/10);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGA  >  minE/1933614‑1933750
                                                                   |                                                                     
ttgTGAACTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCgt                                                                    <  1:2639753/69‑2 (MQ=255)
  gTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTc                                                                  >  1:1442235/1‑71 (MQ=255)
   tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                 <  1:249991/71‑1 (MQ=255)
   tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                 <  1:2740201/71‑1 (MQ=255)
   tGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                 <  1:3220104/71‑1 (MQ=255)
     aaCTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCg                                                                 >  1:3250418/1‑69 (MQ=255)
      aCTCAGGATTATGGCGAACAGAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCaaaa                                                              <  1:3645128/71‑2 (MQ=255)
       ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                             >  1:3670150/1‑71 (MQ=255)
       ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                             <  1:843694/71‑1 (MQ=255)
       ctctGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                             >  1:3504748/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:1274816/1‑58 (MQ=255)
                  gACGTACGGATATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGaa                                                               >  1:1279412/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:2088152/1‑71 (MQ=255)
                         ggAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCt                                           >  1:1701368/1‑71 (MQ=255)
                         agAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACgtt                                           <  1:1161142/70‑3 (MQ=255)
                         agAAATACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACgtt                                           <  1:3260411/70‑3 (MQ=255)
                          gAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCtt                                          >  1:2154592/1‑71 (MQ=255)
                           aaaTACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGAc                                              <  1:3478955/66‑1 (MQ=255)
                              tACCTTCATTACGGAAGTTACGGTTGATTTCGAATACCCgt                                                                    <  1:3495609/41‑2 (MQ=38)
                                 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)
                                     attACGGAAGTTACGGTTGATTTCGAATACCCGTTc                                                                  <  1:962862/35‑1 (MQ=37)
                                     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)
                                                                   |                                                                     
TGGTGAACTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGA  >  minE/1933614‑1933750

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: