Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,768 A→G 19.7% intergenic (+1536/‑243) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7680AG19.7% 247.1 / 52.8 122intergenic (+1536/‑243)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/52);  new base G (14/10);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 130 total aligned reads displayed.

TTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTC  >  minE/2572700‑2572837
                                                                    |                                                                     
ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGat                                                                      >  1:3171790/1‑70 (MQ=255)
  ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                                                                    >  1:4211728/1‑70 (MQ=255)
  ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                                                                    <  1:5065269/70‑1 (MQ=255)
      ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                                                                >  1:2631853/1‑70 (MQ=255)
        acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                              >  1:2439551/1‑70 (MQ=255)
         cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa                                                                       <  1:3796863/60‑1 (MQ=255)
         cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                              <  1:4763221/69‑1 (MQ=255)
          ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                           <  1:5513053/71‑1 (MQ=255)
             ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                         <  1:5045414/70‑1 (MQ=255)
              cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                         >  1:4094487/1‑69 (MQ=255)
              cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                         >  1:2144679/1‑69 (MQ=255)
              cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                                                        <  1:2124118/70‑1 (MQ=255)
                 ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                    <  1:1764035/71‑1 (MQ=255)
                 cccAATCGGGTCCGGGGTGATATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcc                                                    >  1:64164/1‑71 (MQ=255)
                  cccaTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGccc                                                   >  1:1898110/4‑71 (MQ=255)
                  ccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGccc                                                   >  1:2761371/1‑71 (MQ=255)
                   cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                  <  1:28268/71‑1 (MQ=255)
                    caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                                 <  1:2543541/70‑1 (MQ=255)
                     aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGa                                                >  1:3363095/1‑71 (MQ=255)
                      tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                               >  1:2551937/1‑71 (MQ=255)
                        ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                              >  1:5074099/1‑70 (MQ=255)
                           tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                          <  1:2827269/70‑1 (MQ=255)
                           tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                          <  1:5675327/70‑1 (MQ=255)
                            ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                          <  1:5209863/54‑1 (MQ=255)
                            ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatca                                         >  1:1458927/1‑71 (MQ=255)
                             cTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                        <  1:1215386/71‑1 (MQ=255)
                              tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                          <  1:1550200/52‑1 (MQ=255)
                              tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                                                        >  1:3606836/1‑54 (MQ=255)
                              tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                        >  1:4255053/1‑70 (MQ=255)
                              tGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                        >  1:4848809/1‑70 (MQ=255)
                               gggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcatca                                      <  1:3320751/69‑1 (MQ=255)
                               gAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAg                                                                 <  1:5317765/44‑1 (MQ=255)
                                ggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACTcctcatcat                                     <  1:3500579/70‑1 (MQ=255)
                                aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAg                                                                 <  1:375271/43‑1 (MQ=255)
                                aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:4024081/70‑1 (MQ=255)
                                aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                     <  1:9140/71‑1 (MQ=255)
                                 gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                              >  1:4311450/1‑61 (MQ=255)
                                 gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                              >  1:3075723/1‑61 (MQ=255)
                                 gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                              >  1:1850727/1‑61 (MQ=255)
                                  taatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                        <  1:1756528/63‑1 (MQ=255)
                                   gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                        <  1:3750460/65‑1 (MQ=255)
                                   gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                        <  1:3995348/65‑1 (MQ=255)
                                   aatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCAtgct                                  >  1:3625482/3‑71 (MQ=255)
                                    ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                  <  1:459614/54‑1 (MQ=255)
                                    ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                  <  1:5120786/54‑1 (MQ=255)
                                    ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                     <  1:1729571/67‑1 (MQ=255)
                                    ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                  <  1:1963142/70‑1 (MQ=255)
                                    ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                  <  1:704595/70‑1 (MQ=255)
                                    atataCCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcccc                                                  >  1:1408611/2‑53 (MQ=255)
                                    atataCCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACTCCTCATCAtgctg                                 <  1:1260426/70‑2 (MQ=255)
                                    atataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                          >  1:2377207/2‑62 (MQ=255)
                                      ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtata                                                            >  1:3483156/1‑42 (MQ=255)
                                      ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtata                                                            >  1:1599657/1‑42 (MQ=255)
                                      ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtat                                                             >  1:4843679/1‑41 (MQ=255)
                                      ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATagtact                               >  1:311533/1‑71 (MQ=255)
                                       taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                    >  1:3381715/1‑49 (MQ=255)
                                       taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                    >  1:3598182/1‑49 (MQ=255)
                                       taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCATGctgctg                              >  1:1084867/1‑71 (MQ=255)
                                         ccGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTActgct                            <  1:3679169/71‑1 (MQ=255)
                                            ccccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg                         >  1:1675088/1‑71 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                                >  1:4879658/1‑46 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                                >  1:707945/1‑46 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                                >  1:4672802/1‑46 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                  <  1:4925133/60‑1 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactactact                            <  1:2410494/66‑1 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCgg                        >  1:1158939/1‑70 (MQ=255)
                                              ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                       <  1:527883/71‑1 (MQ=255)
                                               ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:3793395/55‑1 (MQ=255)
                                               ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:878015/55‑1 (MQ=255)
                                               ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:2079703/55‑1 (MQ=255)
                                               ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:1131489/55‑1 (MQ=255)
                                               ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCa                      <  1:3059692/71‑1 (MQ=255)
                                                tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:303116/54‑1 (MQ=255)
                                                tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:4509804/54‑1 (MQ=255)
                                                tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:4503456/54‑1 (MQ=255)
                                                tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:105533/54‑1 (MQ=255)
                                                tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                      <  1:3481780/54‑1 (MQ=255)
                                                tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                     >  1:3132992/1‑71 (MQ=255)
                                                 gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                        <  1:1960161/51‑1 (MQ=255)
                                                 gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                        <  1:946158/51‑1 (MQ=255)
                                                 gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                       >  1:3048195/1‑68 (MQ=255)
                                                  ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGt                    <  1:4689544/70‑1 (MQ=255)
                                                    gcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                       >  1:5064551/1‑65 (MQ=255)
                                                     cGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtg                 <  1:4764864/70‑1 (MQ=255)
                                                       gTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgatg              <  1:3982140/71‑1 (MQ=255)
                                                        ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat               <  1:2148182/69‑1 (MQ=255)
                                                        ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat               <  1:3604419/69‑1 (MQ=255)
                                                        ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat               <  1:5478436/69‑1 (MQ=255)
                                                         tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa            >  1:610744/1‑71 (MQ=255)
                                                         tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa            >  1:587680/1‑71 (MQ=255)
                                                          aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactact                               >  1:5012677/1‑51 (MQ=255)
                                                          aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa            >  1:1630218/1‑70 (MQ=255)
                                                            tGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtga                >  1:1852637/1‑64 (MQ=255)
                                                               tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                    >  1:143625/1‑41 (MQ=255)
                                                               tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                    >  1:2378101/1‑41 (MQ=255)
                                                               tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                    >  1:379287/1‑41 (MQ=255)
                                                                 ccGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                  <  1:4970106/41‑1 (MQ=255)
                                                                 ccGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTa         >  1:1186802/1‑66 (MQ=255)
                                                                 ccGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCt    >  1:678455/1‑71 (MQ=255)
                                                                   gATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCttt  <  1:2099569/71‑2 (MQ=255)
                                                                    |                                                                     
TTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTC  >  minE/2572700‑2572837

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: