Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,122,104 C→A 28.5% intergenic (+33/+81) yraP → / ← yraQ hypothetical protein/predicted permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,122,1040CA28.5% 132.3 / 48.6 91intergenic (+33/+81)yraP/yraQhypothetical protein/predicted permease
Reads supporting (aligned to +/- strand):  ref base C (40/25);  new base A (9/17);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 3.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAT  >  minE/2122039‑2122164
                                                                 |                                                            
cGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCcgcg                                                           >  1:968474/1‑69 (MQ=255)
cGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCcgcg                                                           >  1:2382637/1‑69 (MQ=255)
cGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCcgcg                                                           >  1:2449500/1‑69 (MQ=255)
cGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCcgcg                                                           >  1:1630521/1‑69 (MQ=255)
cGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCcgcg                                                           >  1:142835/1‑69 (MQ=255)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:725578/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:1910455/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:2670112/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:2618662/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:750544/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:2850775/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:829008/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:3047725/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:365835/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:507930/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:2441750/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:1526564/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:2393759/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:825095/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:2298668/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:593634/69‑1 (MQ=38)
  ggTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCt                                                        <  1:61284/69‑1 (MQ=38)
            gCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                               <  1:1631147/69‑1 (MQ=255)
            gCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                               <  1:2892143/69‑1 (MQ=255)
            gCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                               <  1:672239/69‑1 (MQ=255)
            gCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                               <  1:1898642/69‑1 (MQ=255)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:2485855/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:2840175/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:2807964/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:2811926/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:2816663/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:2823748/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:523639/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:1503075/1‑69 (MQ=21)
                 tacgtttaTTAAATAGCAGCCCTTGTAATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTACAAc                                         >  1:1424804/1‑69 (MQ=21)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:2741980/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:2737845/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:855760/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:2386239/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:1068275/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:1083570/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:1371224/68‑1 (MQ=255)
                              tAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:1601892/68‑1 (MQ=255)
                                          gTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCg                              <  1:251922/56‑1 (MQ=255)
                                              tGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCg              >  1:21561/1‑68 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCTTAGGTCAGATAAGGCGttt           >  1:1590723/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:35188/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:170833/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2878169/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1659332/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2905808/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:23644/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2973060/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2976917/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:3106398/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:3285370/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:3350021/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2907771/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:674770/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1250255/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:773875/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:837740/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:946746/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1030179/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1755711/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2287663/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2124976/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2099067/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2796905/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2767876/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1907027/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1940030/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:1753468/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGttt           >  1:2668276/1‑69 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGtt            >  1:1157982/1‑68 (MQ=255)
                                                ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGtt            >  1:250291/1‑68 (MQ=255)
                                                     tGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgc       >  1:390504/1‑68 (MQ=255)
                                                     tGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgc       >  1:125058/1‑68 (MQ=255)
                                                     tGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgc       >  1:2533517/1‑68 (MQ=255)
                                                     tGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgc       >  1:1890447/1‑68 (MQ=255)
                                                      gCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgccg     <  1:1780945/69‑1 (MQ=255)
                                                      gCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgccg     <  1:1531513/69‑1 (MQ=255)
                                                      gCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgccg     <  1:1261367/69‑1 (MQ=255)
                                                      gCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgccg     <  1:2434041/69‑1 (MQ=255)
                                                      gCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgccg     <  1:2749247/69‑1 (MQ=255)
                                                       cGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcgccg     <  1:2052635/68‑1 (MQ=255)
                                                         aCGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAt  >  1:2041814/1‑69 (MQ=255)
                                                         aCGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAt  >  1:1251851/1‑69 (MQ=255)
                                                         aCGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAt  >  1:291073/1‑69 (MQ=255)
                                                         aCGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAt  >  1:247162/1‑69 (MQ=255)
                                                           gCTTGCCGCGTCTTATCAGGCCTACAACGATACaa                                  <  1:2348367/35‑1 (MQ=255)
                                                           gCTTGCCGCGTCTTATCAGGCCTACAACGATACaa                                  <  1:1781399/35‑1 (MQ=255)
                                                           gCTTGCCGCGTCTTATCAGGCCTACAACGATACaa                                  <  1:1099307/35‑1 (MQ=255)
                                                           gCTTGCCGCGTCTTATCAGGCCTACAACGATACaa                                  <  1:3300417/35‑1 (MQ=255)
                                                           gCTTGCCGCGTCTTATCAGGCCTACAACGATACaa                                  <  1:1565443/35‑1 (MQ=255)
                                                           gCTTGCCGCGTCTTATCAGGCCTACAACGATACaa                                  <  1:3352640/35‑1 (MQ=255)
                                                                 |                                                            
CGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAT  >  minE/2122039‑2122164

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: