Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 491,553 A→G 41.1% S74P (TCG→CCG)  ybgS ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE491,5530AG41.1% 24.7 / 52.2 73S74P (TCG→CCG) ybgSconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/31);  new base G (12/18);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TTTGCACATGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACATTTT  >  minE/491485‑491621
                                                                    |                                                                    
tttGCACATGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGa                                                                      >  1:2054328/1‑69 (MQ=255)
        tGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATgcagca                                                              >  1:3118543/1‑69 (MQ=255)
        tGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATgcagca                                                              >  1:2145434/1‑69 (MQ=255)
        tGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATgcagca                                                              >  1:2113499/1‑69 (MQ=255)
        tGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATgcagc                                                               >  1:3108860/1‑68 (MQ=255)
        tGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATgcagc                                                               >  1:1658462/1‑68 (MQ=255)
         ggTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCAt                                                             <  1:2603478/69‑1 (MQ=255)
         ggTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCAt                                                             <  1:734273/69‑1 (MQ=255)
         ggTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCAt                                                             <  1:2109600/69‑1 (MQ=255)
         ggTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCAt                                                             <  1:762395/69‑1 (MQ=255)
          gtgtTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCAtt                                                            >  1:1536689/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGAt                                                                    >  1:2734975/1‑48 (MQ=39)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTgcc                                                       >  1:3260901/1‑61 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:2286076/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:642777/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:446788/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:1942571/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:1686868/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:1913039/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:508155/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:3238984/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:2589075/1‑69 (MQ=255)
                       ctTCATCTTTGGTCATGCCGTCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAAc                                               >  1:1233973/1‑69 (MQ=255)
                          cATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCgg                                            >  1:2531795/1‑69 (MQ=255)
                          cATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCgg                                            >  1:2833525/1‑69 (MQ=255)
                            tCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGt                                           <  1:2905716/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGt                                           <  1:2765522/67‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:2728220/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:2686655/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:1355129/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:71393/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:36467/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:895279/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:2439743/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:2497568/68‑1 (MQ=255)
                             cTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTa                                          <  1:1765737/68‑1 (MQ=255)
                                  gTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTAtt                                        <  1:2646603/65‑1 (MQ=255)
                                  gTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTAtt                                        <  1:1365132/65‑1 (MQ=255)
                                  gTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTAc                                     <  1:967966/68‑1 (MQ=255)
                                      ttCCGGCATTGTTCATTGAAGAACCATCCGAATGCAGCAGTGTGCCGCCAGACCCGGt                                           <  1:261442/58‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:3142893/62‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:3144094/62‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:3280036/62‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:1011490/62‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:2564687/62‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:531736/62‑1 (MQ=255)
                                        ccGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGTCGCCAGAACCGGTATTTAc                                     <  1:1523151/62‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCATAACCGGTATTTACCCCGt                                <  1:1680010/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:2988470/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:3315829/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:3301363/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:935190/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:931244/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:1231692/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:1595711/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:737853/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:3090178/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:1775606/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:1840546/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:547563/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:2804408/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:28037/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:524764/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:2708141/64‑1 (MQ=255)
                                           tCATTGTTCATTGAAGAACCATCCGGATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGt                                <  1:2703405/64‑1 (MQ=255)
                                                gTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGAc                      >  1:2573297/1‑69 (MQ=255)
                                                  tCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTAc                                     <  1:3330725/52‑1 (MQ=255)
                                                  tCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTAc                                     <  1:995023/52‑1 (MQ=255)
                                                  tCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTAc                                     <  1:2512587/52‑1 (MQ=255)
                                                  tCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTAc                                     <  1:1270871/52‑1 (MQ=255)
                                                  tCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTAc                                     <  1:3310366/52‑1 (MQ=255)
                                                       gaagaaCCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTAtt                >  1:2805533/1‑68 (MQ=255)
                                                       gaagaaCCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTAtt                >  1:2410083/1‑68 (MQ=255)
                                                                  cGAATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAt     >  1:763677/1‑68 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:703413/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1031397/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1116987/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1204100/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1205367/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1230036/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1666977/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:1682905/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:2133135/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:661910/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:2422985/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:54688/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:2476519/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:2595640/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:2983911/69‑1 (MQ=255)
                                                                    aaTGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACAtttt  <  1:3194690/69‑1 (MQ=255)
                                                                    |                                                                    
TTTGCACATGGTGTTCTTGTGCTCTTCATCTTTGGTCATTCCGTCATTGTTCATTGAAGAACCATCCGAATGCAGCATTGTGCCGCCAGAACCGGTATTTACCCCGTTATTGTCGACGTTATTTGGCGCGACATTTT  >  minE/491485‑491621

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