Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,986 G→A 40.9% L153L (CTG→TTG) ‡ lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9860GA40.9% 55.3 / 99.1 93L153L (CTG→TTG) ‡lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (21/34);  new base A (16/22);  total (37/56)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTG  >  minE/1933918‑1934056
                                                                    |                                                                      
gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                                                      <  1:3141770/71‑1 (MQ=255)
gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                                                      <  1:4376784/71‑1 (MQ=255)
gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                                                      <  1:1631048/71‑1 (MQ=255)
gatgagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCg                                                                      <  1:468999/71‑1 (MQ=255)
   gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                                                                    >  1:3111569/1‑70 (MQ=255)
   gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGa                                                                    >  1:4377217/1‑70 (MQ=255)
   gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                                                                   >  1:2724307/1‑71 (MQ=255)
   gagatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                                                                   >  1:2646964/1‑71 (MQ=255)
    agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAc                                                                   <  1:570914/70‑1 (MQ=255)
      atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGc                                                                 >  1:3726896/1‑70 (MQ=255)
        cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                              >  1:5245815/1‑71 (MQ=255)
        cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                              >  1:462623/1‑71 (MQ=255)
        cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                              >  1:2749457/1‑71 (MQ=255)
        cGAGATAACGCTGACGAT‑‑CGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                                                            >  1:421839/1‑71 (MQ=255)
          agaTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                                                            <  1:2135753/71‑1 (MQ=255)
            aTAACGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACgt                                                                 <  1:5205120/64‑2 (MQ=255)
              aaCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtt                                                          >  1:3160493/1‑69 (MQ=255)
              aaCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtt                                                          >  1:4850361/1‑69 (MQ=255)
               aCGTTGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCtttag                                                       <  1:3315533/71‑3 (MQ=255)
                 gttGACGATAACGGACTTCCTGATCCTGCAGACCATGGGATTTATCTGGTAAAGGACGTAGTGCTTTAGTc                                                     <  1:880468/69‑1 (MQ=255)
                 gCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCtt                                                          <  1:3626070/66‑1 (MQ=255)
                  ttGACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCttta                                                        <  1:4037612/66‑2 (MQ=255)
                  cTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCtt                                                          <  1:4503008/65‑1 (MQ=255)
                    gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                   <  1:408624/70‑1 (MQ=255)
                    gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                   >  1:460533/1‑70 (MQ=255)
                    gACGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                   >  1:1287451/1‑70 (MQ=255)
                     aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCgg                                                                     <  1:3369224/51‑1 (MQ=255)
                     aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCgg                                                                     <  1:2177350/51‑1 (MQ=255)
                     aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca                                                 <  1:871873/71‑1 (MQ=255)
                     aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca                                                 <  1:5361150/71‑1 (MQ=255)
                     aCGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcagc                                                  <  1:4004532/70‑1 (MQ=255)
                        aTAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACg                                                                  <  1:3778771/51‑1 (MQ=255)
                        aTAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                                                            <  1:4950430/57‑1 (MQ=255)
                           gcgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGc                                            <  1:3250009/70‑1 (MQ=255)
                             gcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAAc                                         >  1:3602980/1‑71 (MQ=255)
                              gacTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCt                                        <  1:4835933/69‑1 (MQ=255)
                              gacTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCt                                        <  1:3181762/69‑1 (MQ=255)
                              cgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGt                                                      >  1:2058026/1‑57 (MQ=255)
                              cgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                        >  1:1744131/1‑71 (MQ=255)
                               gcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTc                                       <  1:3085453/71‑1 (MQ=255)
                               acTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                   <  1:3137036/58‑1 (MQ=255)
                               acTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                   <  1:161760/58‑1 (MQ=255)
                                cTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCACCt                                        >  1:1213546/1‑69 (MQ=255)
                                 ttCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                            <  1:3471179/64‑1 (MQ=255)
                                 ttCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                            <  1:4853049/64‑1 (MQ=255)
                                  tCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCt                                                           <  1:806994/48‑1 (MQ=255)
                                   ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCttta                                                        >  1:5470677/1‑49 (MQ=255)
                                   ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTcag                                     <  1:1144142/69‑3 (MQ=255)
                                   ccTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                   >  1:5676856/1‑70 (MQ=255)
                                   ccTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGAc                                              >  1:94270/1‑60 (MQ=255)
                                    cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTca                                      >  1:4077343/1‑66 (MQ=255)
                                    cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                  <  1:1974970/71‑3 (MQ=255)
                                    cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                  <  1:4136137/71‑3 (MQ=255)
                                    cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGta                                   <  1:4003635/70‑2 (MQ=255)
                                     tGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCa                                           >  1:2312386/1‑61 (MQ=255)
                                     tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTg                                   >  1:285824/1‑69 (MQ=255)
                                       aTCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                               <  1:4825226/71‑1 (MQ=255)
                                       aTCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTgg                                                       <  1:2292799/47‑1 (MQ=255)
                                       aTCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGGGCAGt                               >  1:1274071/1‑71 (MQ=255)
                                       aTCCTGCAAGCCGTGGAATTTATACGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGt                               >  1:2917538/1‑71 (MQ=255)
                                        tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                               >  1:3631352/1‑70 (MQ=255)
                                        tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                               >  1:2325026/1‑70 (MQ=255)
                                        tCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGt                               >  1:4916423/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)
                                          cTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCttt                                                         <  1:1921866/42‑1 (MQ=255)
                                                 ccGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                     <  1:4668860/71‑1 (MQ=255)
                                                  catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                    >  1:1855730/3‑70 (MQ=255)
                                                  catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                    >  1:4588118/3‑70 (MQ=255)
                                                  catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                    >  1:2573148/3‑70 (MQ=255)
                                                  cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACa                      <  1:4822578/69‑1 (MQ=255)
                                                   atGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTg                              <  1:1258198/59‑1 (MQ=255)
                                                    tGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                            <  1:3170555/45‑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)
                                                       aaTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTc                                       <  1:386645/47‑1 (MQ=255)
                                                         tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                        <  1:2585115/44‑1 (MQ=255)
                                                         tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                        <  1:3768983/44‑1 (MQ=255)
                                                         tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACt                                        <  1:5386285/44‑1 (MQ=255)
                                                         tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                             <  1:3767021/55‑1 (MQ=255)
                                                         tttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCgg             >  1:2880493/1‑71 (MQ=255)
                                                          ttATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                       <  1:2890008/60‑1 (MQ=255)
                                                           tATCCCGCAGCGGACGCAGTGCTTTGGTCAGCAGACg                                             >  1:4659781/1‑37 (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:5415451/1‑71 (MQ=255)
                                                                 gCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt     >  1:1499819/1‑71 (MQ=255)
                                                                   aGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggtttt         >  1:4613305/1‑65 (MQ=255)
                                                                    gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                        <  1:1941728/49‑1 (MQ=255)
                                                                    gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                        <  1:861621/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)
                                                                    |                                                                      
GATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTG  >  minE/1933918‑1934056

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: