Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,773 A→G 19.6% intergenic (+1541/‑238) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7730AG19.6% 245.2 / 52.1 122intergenic (+1541/‑238)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 = 4.64e-01

Only 100 of 131 total aligned reads displayed.

CCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTG  >  minE/2572705‑2572839
                                                                    |                                                                  
cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                                                                  <  1:2116581/71‑1 (MQ=255)
   acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                                >  1:2439551/1‑70 (MQ=255)
    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                                <  1:5310379/69‑1 (MQ=255)
    cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                                <  1:4763221/69‑1 (MQ=255)
     ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                             <  1:5513053/71‑1 (MQ=255)
     ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                             >  1:1943673/1‑71 (MQ=255)
       tCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                             >  1:1695969/1‑69 (MQ=255)
        ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                            <  1:5113808/69‑1 (MQ=255)
         cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                           >  1:4094487/1‑69 (MQ=255)
         cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                           >  1:1811317/1‑69 (MQ=255)
         cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                                                          <  1:2124118/70‑1 (MQ=255)
            ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      >  1:2273757/1‑71 (MQ=255)
             ccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGccc                                                     >  1:2761371/1‑71 (MQ=255)
              cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                    <  1:28268/71‑1 (MQ=255)
               caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                                   <  1:3753246/70‑1 (MQ=255)
               caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                                   <  1:2543541/70‑1 (MQ=255)
                aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGa                                                  >  1:3363095/1‑71 (MQ=255)
                aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCACTGa                                                  >  1:3547067/1‑71 (MQ=255)
                 tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                                 >  1:4835131/1‑71 (MQ=255)
                 tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                                 >  1:2551937/1‑71 (MQ=255)
                   ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                                >  1:5074099/1‑70 (MQ=255)
                   ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc                                               <  1:3677966/71‑1 (MQ=255)
                      tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                            <  1:5675327/70‑1 (MQ=255)
                      tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                            <  1:2827269/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)
                          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)
                           aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        >  1:4364238/1‑70 (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)
                            gtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcat                                          <  1:4888930/67‑1 (MQ=255)
                             taatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                          <  1:1756528/63‑1 (MQ=255)
                              gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:3750460/65‑1 (MQ=255)
                              aatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCAtgct                                    >  1:3625482/3‑71 (MQ=255)
                               ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                    <  1:459614/54‑1 (MQ=255)
                               ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                    <  1:4608315/54‑1 (MQ=255)
                               ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                       <  1:1729571/67‑1 (MQ=255)
                               ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                    <  1:704595/70‑1 (MQ=255)
                               ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                    <  1:1963142/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:1599657/1‑42 (MQ=255)
                                 ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATagtact                                 >  1:311533/1‑71 (MQ=255)
                                  taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      >  1:3598182/1‑49 (MQ=255)
                                  taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      >  1:3381715/1‑49 (MQ=255)
                                  taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCATGctgctg                                >  1:1084867/1‑71 (MQ=255)
                                   aCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:3216827/60‑1 (MQ=255)
                                       ccccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg                           >  1:1675088/1‑71 (MQ=255)
                                         ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                                  >  1:4879658/1‑46 (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:1131489/55‑1 (MQ=255)
                                          ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:878015/55‑1 (MQ=255)
                                          ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCa                        <  1:3059692/71‑1 (MQ=255)
                                           tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:105533/54‑1 (MQ=255)
                                           tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:4509804/54‑1 (MQ=255)
                                           tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:3481780/54‑1 (MQ=255)
                                           tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:303116/54‑1 (MQ=255)
                                           tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                       >  1:3132992/1‑71 (MQ=255)
                                            gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:1960161/51‑1 (MQ=255)
                                            gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                         >  1:3048195/1‑68 (MQ=255)
                                             ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGt                      <  1:4689544/70‑1 (MQ=255)
                                             ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTc                     >  1:67307/1‑71 (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)
                                                    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)
                                                          tACCGATAATAGTTATAGTACAGCCCTTATTCATCATCAta                                      >  1:5017773/1‑41 (MQ=255)
                                                          tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                      >  1:2378101/1‑41 (MQ=255)
                                                          tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                      >  1:143625/1‑41 (MQ=255)
                                                          tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                      >  1:379287/1‑41 (MQ=255)
                                                          tACCGATAATAGTTATAGTACAGCCCTGATTAATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTg        >  1:3958912/1‑71 (MQ=255)
                                                           aCCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtga                  >  1:86042/1‑60 (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)
                                                                taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta  >  1:548902/1‑68 (MQ=255)
                                                                taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta  >  1:5567465/1‑68 (MQ=255)
                                                                taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta  <  1:3806569/71‑4 (MQ=255)
                                                                  ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg                           >  1:2202934/1‑44 (MQ=255)
                                                                  ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                       <  1:2247157/48‑1 (MQ=255)
                                                                  ataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg                       <  1:4014853/48‑1 (MQ=255)
                                                                    |                                                                  
CCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTG  >  minE/2572705‑2572839

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: