Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,118,446 G→A 5.9% P164S (CCA→TCA)  ydiU ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,118,4460GA 5.9% 205.3 / 2.5 67P164S (CCA→TCA) ydiUconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base A (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCC  >  minE/1118376‑1118515
                                                                      |                                                                     
cggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTgg                                                                       <  1:1890978/71‑1 (MQ=255)
cggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTgg                                                                       <  1:4574089/71‑1 (MQ=255)
 ggTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGc                                                                      >  1:4319783/1‑71 (MQ=255)
       aaaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgcc                                                                >  1:2341774/1‑71 (MQ=255)
        aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                               >  1:934468/1‑71 (MQ=255)
        aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                               >  1:2702688/1‑71 (MQ=255)
        aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTccgccg                                                               >  1:849264/1‑71 (MQ=255)
        aaaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                    >  1:274434/1‑66 (MQ=255)
         aaaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGt                                                              <  1:4438858/71‑1 (MQ=255)
          aaTGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGtt                                                             >  1:4217380/1‑71 (MQ=255)
            tGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTc                                                           >  1:2794939/1‑71 (MQ=255)
              ttCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                         <  1:2121411/71‑1 (MQ=255)
                cGAAATTACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                        <  1:2178618/70‑1 (MQ=255)
                cGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCg                                                        <  1:1764701/70‑1 (MQ=255)
                cGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                       >  1:1689918/1‑71 (MQ=255)
                cGAAATGACCAAAGCGCAGATGACTTGGTGCAACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                       >  1:3379215/1‑71 (MQ=255)
                  aaaTGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATa                                                     >  1:2361756/1‑71 (MQ=255)
                       aCCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTg                                                <  1:3155759/71‑1 (MQ=255)
                          aaaGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTccc                                                         <  1:4083480/59‑1 (MQ=255)
                            aGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                    >  1:4126212/1‑46 (MQ=255)
                              cgcAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCTTTTCCCGATACACTGGCGAATc                                         <  1:2357817/71‑1 (MQ=255)
                                cAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCg                                        <  1:2794847/70‑1 (MQ=255)
                                 aGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTcc                                                          >  1:1650615/1‑51 (MQ=255)
                                   aTGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTc                                                                    >  1:2647574/1‑39 (MQ=255)
                                    tGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATaca                                                   >  1:2616093/1‑55 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTcc                                                          >  1:3202578/1‑47 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGa                                                       >  1:4550111/1‑50 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                   >  1:3649497/1‑70 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                   >  1:2817458/1‑70 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                   >  1:1691475/1‑70 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGt                                   >  1:562364/1‑70 (MQ=255)
                                     gACTTGGTGCCACACGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGCTGGt                                   <  1:714957/70‑1 (MQ=255)
                                      aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                  <  1:4536341/70‑1 (MQ=255)
                                      aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                  <  1:1694466/70‑1 (MQ=255)
                                      aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                 <  1:4274794/71‑1 (MQ=255)
                                      aCTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGa                                 <  1:280709/71‑1 (MQ=255)
                                       cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTg                                  <  1:4316446/69‑1 (MQ=255)
                                       cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGAc                                <  1:139018/71‑1 (MQ=255)
                                         tGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                              >  1:3008774/1‑71 (MQ=255)
                                          ggTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                              >  1:348860/1‑70 (MQ=255)
                                           gTGCCACACGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGCTGGTGACGAt                             <  1:3725507/70‑1 (MQ=255)
                                              ccACACGCATCAGCATCTCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACtt                         <  1:3090433/71‑1 (MQ=255)
                                                acacGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGc                                       >  1:2334425/1‑55 (MQ=255)
                                                acacGCATCAGCATCGCGCCTGACTCCACCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTa                        >  1:753958/1‑70 (MQ=255)
                                                  acGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcg                     >  1:3109144/1‑71 (MQ=255)
                                                   cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGc                                              >  1:2168606/1‑45 (MQ=255)
                                                   cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATAc                           >  1:1435511/1‑64 (MQ=255)
                                                   cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                    <  1:1617045/71‑1 (MQ=255)
                                                   cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                    <  1:1846979/71‑1 (MQ=255)
                                                   cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                    >  1:4545402/1‑71 (MQ=255)
                                                   cGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgc                    <  1:2171750/71‑1 (MQ=255)
                                                    gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                   >  1:709789/1‑71 (MQ=255)
                                                    gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                   >  1:1716141/1‑71 (MQ=255)
                                                    gCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                   <  1:4379989/71‑1 (MQ=255)
                                                     cATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcgc                  <  1:4565444/71‑1 (MQ=255)
                                                     cATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAcgcg                   <  1:3078624/70‑1 (MQ=255)
                                                        cAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGc                                       <  1:3693050/47‑1 (MQ=255)
                                                        cAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAAc                      >  1:3690947/1‑64 (MQ=255)
                                                         aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGTTGGTGACGATACTTAACGCGCggg               >  1:4446537/1‑70 (MQ=255)
                                                         aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTaa                       <  1:306003/62‑1 (MQ=255)
                                                         aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgt              >  1:2724991/1‑71 (MQ=255)
                                                             tCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGgtcgt           >  1:4296942/1‑70 (MQ=255)
                                                               gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTTACGATACTTAACGCGCGGGTCGTCgg        <  1:1426267/71‑1 (MQ=255)
                                                               gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg        <  1:3197870/71‑1 (MQ=255)
                                                               gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCgg        <  1:933094/71‑1 (MQ=255)
                                                               gcgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTCGTGACGATACTTAACGCGCGGGTCGTCgg        <  1:4412169/71‑1 (MQ=255)
                                                                cgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa       <  1:2947960/71‑1 (MQ=255)
                                                                cgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa       <  1:4305943/71‑1 (MQ=255)
                                                                cgcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGa       <  1:1705868/71‑1 (MQ=255)
                                                                     ggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGa                              >  1:3225689/1‑43 (MQ=255)
                                                                     ggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGcc  >  1:3325683/1‑71 (MQ=255)
                                                                      gCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAAt     >  1:153602/1‑67 (MQ=255)
                                                                      |                                                                     
CGGTAGTAAAAATGTTCGAAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCC  >  minE/1118376‑1118515

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: