Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,699 C→A 16.6% V248V (GTG→GTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6990CA16.6% 185.5 / 16.4 66V248V (GTG→GTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base A (7/4);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGG  >  minE/1933631‑1933768
                                                                    |                                                                     
tGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                                     <  1:1300805/71‑1 (MQ=255)
  aCGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccacc                                                                    <  1:3086059/70‑1 (MQ=255)
   cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                                     >  1:1376230/1‑68 (MQ=255)
   cGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAaccac                                                                     >  1:3630647/1‑68 (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: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)
                 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)
                       aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTcc                                               <  1:495593/70‑1 (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:3080684/1‑71 (MQ=255)
                                    ttGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAg                                 >  1:3582482/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:3801662/1‑71 (MQ=255)
                                         ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa                            >  1:1286058/1‑71 (MQ=255)
                                         ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa                            >  1:1290634/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)
                                                          aGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGcgcg           >  1:2184303/1‑71 (MQ=255)
                                                           gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt          <  1:2949407/71‑1 (MQ=255)
                                                            ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATa                                      >  1:3723678/1‑42 (MQ=255)
                                                            ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt          >  1:3802675/1‑70 (MQ=255)
                                                             caccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAg                                 <  1:3878224/46‑1 (MQ=255)
                                                             caccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAg                                 <  1:1007358/46‑1 (MQ=255)
                                                              accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATAc                                     <  1:489179/41‑1 (MQ=255)
                                                               ccaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGAt                   <  1:2319111/58‑1 (MQ=255)
                                                                caaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtgg     >  1:1888811/1‑71 (MQ=255)
                                                                caaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtgg     <  1:2314149/71‑1 (MQ=255)
                                                                 aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCgaga                    >  1:3477225/1‑55 (MQ=255)
                                                                 aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAAGTACATGTCGAGATCCAGCGCGTTAtggt    >  1:2978053/1‑71 (MQ=255)
                                                                  accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtg      >  1:135577/1‑68 (MQ=255)
                                                                   ccaccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtgg  <  1:2061418/71‑1 (MQ=255)
                                                                   ccaccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtgg  <  1:1134267/71‑1 (MQ=255)
                                                                   ccaccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtgg  <  1:2897533/71‑1 (MQ=255)
                                                                    |                                                                     
TGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGG  >  minE/1933631‑1933768

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: