Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 704,324 G→A 14.3% D107D (GAC→GAT ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE704,3240GA14.3% 204.9 / 21.3 84D107D (GAC→GATymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base A (8/4);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTAT  >  minE/704254‑704391
                                                                      |                                                                   
gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg                                                                     >  1:2473338/1‑71 (MQ=255)
    aCCAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCtt                                                                 <  1:5372011/71‑1 (MQ=255)
     ccAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCtt                                                                 >  1:985316/1‑70 (MQ=255)
     ccAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCtt                                                                 <  1:2763889/70‑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:1750856/1‑71 (MQ=255)
                 tGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATc                                                    >  1:3675203/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:3651759/71‑1 (MQ=255)
                   cGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCTTTATACGTTTGATcgc                                                  <  1:4096640/71‑1 (MQ=255)
                      tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                               <  1:2474168/55‑1 (MQ=255)
                      tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                               <  1:1210079/55‑1 (MQ=255)
                      tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                               <  1:4445834/55‑1 (MQ=255)
                      tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                                                               <  1:1648098/55‑1 (MQ=255)
                           aGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCttta                                                               >  1:1501238/1‑49 (MQ=255)
                           aGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGAAGGCTTTATCttta                                                               >  1:1135682/1‑49 (MQ=255)
                           aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                          <  1:2746868/71‑1 (MQ=255)
                           aGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaaa                                          <  1:2124752/71‑1 (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:1188925/47‑1 (MQ=255)
                             tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCttt                                                                <  1:5454464/47‑1 (MQ=255)
                             tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                           >  1:130526/1‑68 (MQ=255)
                             tAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGaa                                           >  1:1528892/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)
                               aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGt                                                              >  1:5538718/1‑47 (MQ=255)
                               aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt                                      >  1:3100519/1‑71 (MQ=255)
                               aCCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCtt                                      >  1:1807894/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)
                                    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: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:3769852/71‑1 (MQ=255)
                                                  aGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCg                   <  1:5130217/71‑1 (MQ=255)
                                                   gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                             <  1:4574332/44‑1 (MQ=255)
                                                   gTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTg                                             <  1:4611536/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)
                                                                 gATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCg     >  1:4534918/1‑70 (MQ=255)
                                                                  aTTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGta   >  1:3034772/1‑71 (MQ=255)
                                                                  aTTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGta   >  1:877800/1‑71 (MQ=255)
                                                                   tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat  <  1:1405437/71‑1 (MQ=255)
                                                                   tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat  <  1:4540870/71‑1 (MQ=255)
                                                                   tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat  <  1:3776153/71‑1 (MQ=255)
                                                                   tttGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGtat  <  1:1142498/71‑1 (MQ=255)
                                                                      |                                                                   
GTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGCTGAAATCTTCCCACTGGCTATATTTGCGTGTGCGCACATCCGTAT  >  minE/704254‑704391

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: