Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,934,020 G→A 18.5% T141T (ACC→ACT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,934,0200GA18.5% 171.8 / 29.8 81T141T (ACC→ACTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (28/38);  new base A (8/7);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCC  >  minE/1933953‑1934088
                                                                   |                                                                    
ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                                                   >  1:5676856/1‑70 (MQ=255)
 cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                                                  <  1:4136137/71‑3 (MQ=255)
 cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                                                  <  1:1974970/71‑3 (MQ=255)
 cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                                                   <  1:4003635/70‑2 (MQ=255)
  tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTg                                                                   >  1:285824/1‑69 (MQ=255)
    aTCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                               <  1:4825226/71‑1 (MQ=255)
    aTCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGGGCAGt                                                               >  1:1274071/1‑71 (MQ=255)
    aTCCTGCAAGCCGTGGAATTTATACGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGt                                                               >  1:2917538/1‑71 (MQ=255)
     tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                               >  1:4916423/1‑70 (MQ=255)
     tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                               >  1:3631352/1‑70 (MQ=255)
     tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                                                               >  1:2325026/1‑70 (MQ=255)
     tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTg                                                              <  1:4174745/71‑1 (MQ=255)
     tCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                                                              <  1:3500613/71‑1 (MQ=255)
     tCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTg                                                              <  1:1761591/71‑1 (MQ=255)
     tCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGAGCAGTg                                                              <  1:735382/71‑1 (MQ=255)
              ccGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     <  1:4668860/71‑1 (MQ=255)
               catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                                    >  1:2573148/3‑70 (MQ=255)
               catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                                    >  1:1855730/3‑70 (MQ=255)
               catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                                                    >  1:4588118/3‑70 (MQ=255)
               cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa                                                      <  1:4822578/69‑1 (MQ=255)
                atGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTg                                                              <  1:1258198/59‑1 (MQ=255)
                 tGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAGCTc                                                  <  1:3823357/71‑1 (MQ=255)
                 tGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCg                                                                      <  1:2863152/51‑1 (MQ=255)
                  ggAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGa                                                        >  1:2926615/1‑64 (MQ=255)
                   gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGc                                                <  1:191712/71‑1 (MQ=255)
                      tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                                                             <  1:3767021/55‑1 (MQ=255)
                      tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCgg                                             >  1:2880493/1‑71 (MQ=255)
                       ttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                                       <  1:2890008/60‑1 (MQ=255)
                            cGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg                                       <  1:5295937/71‑1 (MQ=255)
                            cGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg                                       <  1:3564349/71‑1 (MQ=255)
                              gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTc                                                  <  1:4178431/58‑1 (MQ=255)
                              gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                                     >  1:1499819/1‑71 (MQ=255)
                              gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                                     >  1:5415451/1‑71 (MQ=255)
                                aGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtttt                                         >  1:4613305/1‑65 (MQ=255)
                                 gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                        <  1:861621/49‑1 (MQ=255)
                                 gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                        <  1:1941728/49‑1 (MQ=255)
                                 gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                        <  1:4663411/49‑1 (MQ=255)
                                 gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg                                  <  1:215619/71‑1 (MQ=255)
                                 gCGGACGCAGTGCTTTGGGCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                        <  1:3394994/49‑1 (MQ=255)
                                  cGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGa                                 <  1:5582937/71‑1 (MQ=255)
                                    gACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAAc                               <  1:4641451/71‑1 (MQ=255)
                                     aCGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGaa                                >  1:2462441/1‑69 (MQ=255)
                                        cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     <  1:3062072/45‑1 (MQ=255)
                                        cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     <  1:3005811/45‑1 (MQ=255)
                                        cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     <  1:4849127/45‑1 (MQ=255)
                                        cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     <  1:1403694/45‑1 (MQ=255)
                                          gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                                       <  1:3277571/41‑1 (MQ=255)
                                          gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                                    <  1:5414204/44‑1 (MQ=255)
                                          gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                                    <  1:5187435/44‑1 (MQ=255)
                                          gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGc                            <  1:5055446/68‑1 (MQ=255)
                                          gTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGc                            <  1:1988341/68‑1 (MQ=255)
                                           tGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtt                                           <  1:2294863/52‑1 (MQ=255)
                                            gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                       >  1:4105097/1‑71 (MQ=255)
                                            gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                       >  1:47924/1‑71 (MQ=255)
                                            gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                       >  1:1074442/1‑71 (MQ=255)
                                            gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                       >  1:854649/1‑71 (MQ=255)
                                            gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAcc                       >  1:4893703/1‑71 (MQ=255)
                                            gCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTAc                        >  1:1799409/1‑70 (MQ=255)
                                             cTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     >  1:5202966/1‑40 (MQ=255)
                                             cTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                                                     >  1:1796875/1‑40 (MQ=255)
                                             cTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCt                      >  1:5323365/1‑71 (MQ=255)
                                               ttGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                                                        <  1:2630518/35‑1 (MQ=255)
                                                tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACa                              >  1:591886/1‑60 (MQ=255)
                                                tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCttt                    <  1:2978676/70‑1 (MQ=255)
                                                tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCtt                     <  1:4657059/69‑1 (MQ=255)
                                                tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTc                   >  1:2457260/1‑71 (MQ=255)
                                                tGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTc                   >  1:5371789/1‑71 (MQ=255)
                                                 ggTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTcg                  >  1:1529964/1‑71 (MQ=255)
                                                     agcagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGcc              <  1:5575288/71‑1 (MQ=255)
                                                     agcagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGcc              >  1:5550570/1‑71 (MQ=255)
                                                       cagTCGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg                                  <  1:848158/49‑1 (MQ=255)
                                                       cagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGgttttggttttg                                  <  1:969624/49‑1 (MQ=255)
                                                       cagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCtt                          <  1:4008047/57‑1 (MQ=255)
                                                        agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg           >  1:3497594/1‑71 (MQ=255)
                                                        agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg           >  1:4204476/1‑71 (MQ=255)
                                                        agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg           >  1:5254303/1‑71 (MQ=255)
                                                        agACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg           <  1:4847325/71‑1 (MQ=255)
                                                         gACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttt                                          <  1:3367423/39‑1 (MQ=255)
                                                         gACGCAACTCGGTGCAGTGGATAGACAGTTAGCCggttt                                          <  1:4013623/39‑1 (MQ=255)
                                                            gCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGa         >  1:1397935/1‑69 (MQ=255)
                                                                 tCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGcc  >  1:4870298/1‑71 (MQ=255)
                                                                  cGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAAc                               >  1:2572662/1‑41 (MQ=255)
                                                                   |                                                                    
CCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCC  >  minE/1933953‑1934088

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: