Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,993 C→T 42.5% L150L (CTG→CTA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,9930CT42.5% 39.7 / 96.0 87L150L (CTG→CTAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (17/33);  new base T (16/21);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

ATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAAC  >  minE/1933924‑1934059
                                                                     |                                                                  
atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGc                                                                    >  1:3726896/1‑70 (MQ=255)
  cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                                 >  1:462623/1‑71 (MQ=255)
  cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                                 >  1:2749457/1‑71 (MQ=255)
  cGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGt                                                                 >  1:5245815/1‑71 (MQ=255)
  cGAGATAACGCTGACGAT‑‑CGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                                                               >  1:421839/1‑71 (MQ=255)
    agaTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGc                                                               <  1:2135753/71‑1 (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)
               aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca                                                    <  1:5361150/71‑1 (MQ=255)
               aCGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTcagca                                                    <  1:871873/71‑1 (MQ=255)
               aCGATAACGGACTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcagc                                                     <  1:4004532/70‑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:3181762/69‑1 (MQ=255)
                        gacTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCt                                           <  1:4835933/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:161760/58‑1 (MQ=255)
                         acTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTcag                                                      <  1:3137036/58‑1 (MQ=255)
                          cTTCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCACCt                                           >  1:1213546/1‑69 (MQ=255)
                           ttCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                               <  1:4853049/64‑1 (MQ=255)
                           ttCCTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGgcgc                                               <  1:3471179/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:4136137/71‑3 (MQ=255)
                              cTGATCCTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGtac                                     <  1:1974970/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:4916423/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)
                                    cTGCAGACCATGGAATTTATCTGGTAAAGGACGTAGTGCttt                                                            <  1:1921866/42‑1 (MQ=255)
                                           ccGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                        <  1:4668860/71‑1 (MQ=255)
                                            catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                       >  1:2573148/3‑70 (MQ=255)
                                            catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                       >  1:4588118/3‑70 (MQ=255)
                                            catGGAATTTATCTGGTAAAGGACGTAGTGCTTTAGTCAGCAGGCGCAGCTCAGTACAGTGAATGGAAAgc                       >  1:1855730/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:3564349/71‑1 (MQ=255)
                                                         cGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg          <  1:5295937/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:4663411/49‑1 (MQ=255)
                                                              gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGa                           <  1:861621/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:4849127/45‑1 (MQ=255)
                                                                     cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                        <  1:1403694/45‑1 (MQ=255)
                                                                     cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                        <  1:3005811/45‑1 (MQ=255)
                                                                     cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAg                        <  1:3062072/45‑1 (MQ=255)
                                                                     |                                                                  
ATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAAC  >  minE/1933924‑1934059

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: