Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 704,311 A→G 16.6% F112L (TTT→CTT) ‡ ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE704,3110AG16.6% 196.7 / 30.0 90F112L (TTT→CTT) ‡ymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base G (8/7);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCG  >  minE/704242‑704380
                                                                     |                                                                     
cAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaa                                                                       <  1:4148528/70‑1 (MQ=255)
  gCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATc                                                                    <  1:1958778/71‑1 (MQ=255)
       atCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGaa                                                                 >  1:2059709/1‑69 (MQ=255)
       atCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGaa                                                                 >  1:4044390/1‑69 (MQ=255)
        tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAAt                                                                >  1:4654760/1‑69 (MQ=255)
          aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAAt                                                                     >  1:2205433/1‑62 (MQ=255)
          aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAtt                                                            >  1:4414167/1‑71 (MQ=255)
          aCGTTTACCAAACGCAACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAAt                                                                     >  1:3409610/1‑62 (MQ=255)
            gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg                                                          >  1:2473338/1‑71 (MQ=255)
              ttACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaaa                                                                      >  1:311831/1‑57 (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)
                  cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaaa                                                                      >  1:3623456/1‑53 (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)
                            ttGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaa                                                                       <  1:4581591/42‑1 (MQ=255)
                            ttGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaa                                                                       <  1:603560/42‑1 (MQ=255)
                            ttGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                     >  1:2282672/1‑60 (MQ=255)
                            ttGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                     >  1:1376976/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)
                              gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAAt                                                                     >  1:570752/1‑42 (MQ=255)
                              gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAAt                                                                     >  1:4923259/1‑42 (MQ=255)
                              gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAAt                                                                     >  1:5088302/1‑42 (MQ=255)
                              gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGaa                                                                 >  1:3399473/1‑46 (MQ=255)
                              gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGaa                                                                 >  1:1240711/1‑46 (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:2124752/71‑1 (MQ=255)
                                       aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                               <  1:2746868/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:5454464/47‑1 (MQ=255)
                                         tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                     <  1:1188925/47‑1 (MQ=255)
                                         tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                     <  1:4772255/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:4357608/1‑71 (MQ=255)
                                                cttcCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCAc                      >  1:3531913/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:1814948/1‑47 (MQ=255)
                                                           cGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCt                                   >  1:2227201/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)
                                                                     |                                                                     
CAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCG  >  minE/704242‑704380

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: