Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 704,318 T→G 17.8% S109S (TCA→TCC) ‡ ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE704,3180TG17.8% 151.5 / 36.8 84S109S (TCA→TCC) ‡ymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/37);  new base G (8/7);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCG  >  minE/704250‑704388
                                                                    |                                                                      
tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAAt                                                                        >  1:4654760/1‑69 (MQ=255)
  aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAtt                                                                    >  1:4414167/1‑71 (MQ=255)
    gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg                                                                  >  1:2473338/1‑71 (MQ=255)
       tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAttt                                                                   <  1:4087923/67‑1 (MQ=255)
        aCCAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCtt                                                              <  1:5372011/71‑1 (MQ=255)
         ccAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCtt                                                              >  1:985316/1‑70 (MQ=255)
         ccAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCtt                                                              <  1:2763889/70‑1 (MQ=255)
          cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAttt                                                                   <  1:4834842/64‑1 (MQ=255)
          cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:2303699/71‑1 (MQ=255)
          cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:4194489/71‑1 (MQ=255)
           aaaCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGt                                                           <  1:5401449/71‑1 (MQ=255)
            accgcgACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTtata                                                          >  1:4471478/3‑71 (MQ=255)
            accGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTtata                                                          >  1:3816527/3‑71 (MQ=255)
              cGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGATTGATTTGTCTTTGTATg                                                        <  1:2826470/71‑1 (MQ=255)
                  aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg                                                    <  1:298368/71‑1 (MQ=255)
                   cTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa                                                   <  1:4845245/71‑1 (MQ=255)
                    ttGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                             >  1:1376976/1‑60 (MQ=255)
                    ttGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                             >  1:2282672/1‑60 (MQ=255)
                    ttGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATc                                                                <  1:4057529/57‑1 (MQ=255)
                     tGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCt                                                               >  1:3091942/1‑57 (MQ=255)
                     tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                 >  1:3675203/1‑71 (MQ=255)
                     tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                 >  1:1750856/1‑71 (MQ=255)
                     tGCGGCAGCCAGTAACTCACTTCCCACAAACGCAGATTGAGATCGAAGGCTTTATCTTtata                                                          >  1:3758542/1‑62 (MQ=255)
                      gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg                                                >  1:2796996/1‑71 (MQ=255)
                       cGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcg                                                <  1:2704464/70‑1 (MQ=255)
                       cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                               <  1:4096640/71‑1 (MQ=255)
                       cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                               <  1:3651759/71‑1 (MQ=255)
                          tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:1210079/55‑1 (MQ=255)
                          tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:2474168/55‑1 (MQ=255)
                          tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:4445834/55‑1 (MQ=255)
                          tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:1648098/55‑1 (MQ=255)
                               aGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCttta                                                            >  1:1135682/1‑49 (MQ=255)
                               aGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCttta                                                            >  1:1501238/1‑49 (MQ=255)
                               aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                       <  1:2746868/71‑1 (MQ=255)
                               aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                       <  1:2124752/71‑1 (MQ=255)
                                gTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCttta                                                                  <  1:5629868/43‑2 (MQ=255)
                                gTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCttta                                                                  <  1:4888498/43‑2 (MQ=255)
                                gTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCttta                                                                  <  1:3257370/43‑2 (MQ=255)
                                gTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                      <  1:3557014/71‑1 (MQ=255)
                                gTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCAGAAAt                                      <  1:2269246/71‑1 (MQ=255)
                                 tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                             <  1:4772255/47‑1 (MQ=255)
                                 tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                             <  1:5454464/47‑1 (MQ=255)
                                 tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                             <  1:1188925/47‑1 (MQ=255)
                                 tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                        >  1:1528892/1‑68 (MQ=255)
                                 tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                        >  1:130526/1‑68 (MQ=255)
                                  aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACt                                                     >  1:3500166/1‑54 (MQ=255)
                                  aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACt                                                     >  1:1236879/1‑54 (MQ=255)
                                  aaCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt                                    >  1:1163262/1‑71 (MQ=255)
                                   aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAttt                                                                   >  1:2115823/1‑39 (MQ=255)
                                   aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGt                                                           >  1:5538718/1‑47 (MQ=255)
                                   aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt                                   >  1:1807894/1‑71 (MQ=255)
                                   aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt                                   >  1:3100519/1‑71 (MQ=255)
                                   aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt                                    >  1:4132175/1‑70 (MQ=255)
                                     cttcttCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTcc                                 >  1:4450791/1‑71 (MQ=255)
                                      ttcttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                      >  1:3540141/1‑65 (MQ=255)
                                      ttcttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCAGCTGAAATCTTcc                                 >  1:2706530/1‑70 (MQ=255)
                                       tcttcCCACAGGCGAAGTTTAAAATCGAGTGAttt                                                                   <  1:2149124/35‑1 (MQ=255)
                                        cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                              >  1:3531913/1‑71 (MQ=255)
                                        cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                              >  1:4357608/1‑71 (MQ=255)
                                         ttcCCACAGGCGAATTTTAAAATCGAATGATTTGTCTTTg                                                            <  1:1648101/40‑1 (MQ=255)
                                              acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTa                                                          >  1:822542/1‑37 (MQ=255)
                                              acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt                                                         <  1:2463480/38‑1 (MQ=255)
                                              acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                        <  1:4293144/55‑1 (MQ=255)
                                              acaGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCta                        >  1:428365/1‑71 (MQ=255)
                                               caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGGTGAAATCTTCCCACTGGc                          <  1:4277656/68‑1 (MQ=255)
                                               caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                      >  1:875528/1‑56 (MQ=255)
                                               caGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCtat                       <  1:1734292/71‑1 (MQ=255)
                                                   cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt                                           >  1:2227201/1‑47 (MQ=255)
                                                   cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt                                           >  1:1814948/1‑47 (MQ=255)
                                                    gAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa                                                   >  1:3783593/1‑38 (MQ=255)
                                                     aaGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAAt                                      <  1:5651640/50‑1 (MQ=255)
                                                      aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                              >  1:4380159/1‑57 (MQ=255)
                                                      aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                <  1:5130217/71‑1 (MQ=255)
                                                      aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                <  1:3769852/71‑1 (MQ=255)
                                                       gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                          <  1:4611536/44‑1 (MQ=255)
                                                       gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                          <  1:4574332/44‑1 (MQ=255)
                                                        tttAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGc                 <  1:3889830/68‑1 (MQ=255)
                                                          tAAAATCGAATGATTTGTCTTTGTATGACTGATcgcc                                              >  1:1557522/1‑37 (MQ=255)
                                                          tAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCa                               >  1:2202457/1‑52 (MQ=255)
                                                          tAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCgtgt             <  1:2930932/70‑1 (MQ=255)
                                                            aaaTCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTcc                                 >  1:3072299/1‑48 (MQ=255)
                                                            aaaTCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTgcg          <  1:4934951/71‑1 (MQ=255)
                                                               tCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGcac       <  1:3229673/71‑1 (MQ=255)
                                                                 gAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGcaca      <  1:4553275/70‑1 (MQ=255)
                                                                   aTGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCt                                    >  1:2240798/1‑38 (MQ=255)
                                                                    tGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg  >  1:5656822/1‑71 (MQ=255)
                                                                    tGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg  >  1:3307530/1‑71 (MQ=255)
                                                                    |                                                                      
TCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCG  >  minE/704250‑704388

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