Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,690 G→A 16.4% G251G (GGC→GGT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,6900GA16.4% 190.1 / 24.5 79G251G (GGC→GGTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base A (7/6);  total (44/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 = 1.00e+00

CTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGT  >  minE/1933621‑1933760
                                                                     |                                                                      
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)
            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: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)
                                   ggAAGTTACGGTTGATTTCGAATACACTCTCGAAGc                                                                       <  1:2014571/36‑1 (MQ=255)
                                   ggAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                       >  1:855054/1‑36 (MQ=255)
                                   ggAAGTTACGGTTGATTTCGAATACACGCTCGAAGc                                                                       <  1:3189086/36‑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:3912495/1‑71 (MQ=255)
                                     aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                                  >  1:2645146/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:1290634/1‑71 (MQ=255)
                                                   ttCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa                    >  1:1286058/1‑71 (MQ=255)
                                                     cGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACAt                   <  1:638973/70‑1 (MQ=255)
                                                     cGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACAt                   <  1:3335801/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)
                                                                     |                                                                      
CTCTGGGTTATGACGTACGGAAATACCTTCGTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGT  >  minE/1933621‑1933760

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: