Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,941 G→A 20.4% intergenic (+1709/‑70) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,9410GA20.4% 194.0 / 42.9 88intergenic (+1709/‑70)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/29);  new base A (10/8);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GATAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTT  >  minE/2572872‑2573009
                                                                     |                                                                    
gATAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGgaga                                                                     <  1:1281686/71‑1 (MQ=255)
 aTAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGAt                                                                    >  1:1447287/1‑71 (MQ=255)
  tAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGAtt                                                                   >  1:492922/1‑71 (MQ=255)
   aaaTCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTc                                                                  >  1:771026/1‑71 (MQ=255)
    aaTCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTca                                                                 <  1:3636097/71‑2 (MQ=255)
     aTCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGt                                                                <  1:821986/71‑1 (MQ=255)
      tCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGt                                                                >  1:177536/1‑70 (MQ=255)
        aTTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGt                                                                >  1:1110705/1‑68 (MQ=255)
         tttGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCg                                                                 >  1:3902057/1‑66 (MQ=255)
         tttGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGt                                                                >  1:246649/1‑67 (MQ=255)
         tttGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGt                                                                <  1:1264292/67‑1 (MQ=255)
         tttGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGa                                                            <  1:2472896/71‑1 (MQ=255)
           tggttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGat                                                           <  1:1002078/70‑1 (MQ=255)
             gttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGatat                                                         <  1:57418/70‑1 (MQ=255)
              ttggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCa                                                       <  1:63747/71‑1 (MQ=255)
                ggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGTCTTAAGGAGATTCGTAAGATATCAGc                                                     >  1:1727338/1‑71 (MQ=255)
                ggCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAg                                                      <  1:739253/70‑1 (MQ=255)
                  cgcACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGGTATCAGCCa                                                   >  1:1727925/1‑71 (MQ=255)
                   gcACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCAc                                                  <  1:1078111/71‑1 (MQ=255)
                       tGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGata                                                          <  1:1371497/59‑1 (MQ=255)
                       tGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGata                                                          <  1:3737198/59‑1 (MQ=255)
                       tGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCa                                                   <  1:2788121/66‑1 (MQ=255)
                        ggAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACta                                                >  1:3208371/1‑68 (MQ=255)
                         gAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACt                                                 >  1:737501/1‑66 (MQ=255)
                         gAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATAc                                             <  1:971236/70‑1 (MQ=255)
                            taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGat                                         <  1:1326741/71‑1 (MQ=255)
                            taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGat                                         <  1:2784413/71‑1 (MQ=255)
                            taatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGat                                         <  1:3307480/71‑1 (MQ=255)
                             aatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCa                                                   >  1:920237/1‑60 (MQ=255)
                             aatCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGat                                         >  1:3126183/1‑70 (MQ=255)
                              atCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCa                                                   >  1:4004731/1‑59 (MQ=255)
                              atCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACta                                                >  1:1679580/1‑62 (MQ=255)
                              atCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGata                                        >  1:2179051/1‑70 (MQ=255)
                              atCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGata                                        >  1:1247628/1‑70 (MQ=255)
                              atCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGata                                        >  1:3046975/1‑70 (MQ=255)
                                 tGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGata                                                          >  1:808254/1‑49 (MQ=255)
                                 tGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAg                                                      >  1:1833989/1‑53 (MQ=255)
                                 tGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGa                                                            <  1:1513734/47‑1 (MQ=255)
                                   aGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGAtata                                      <  1:569023/67‑1 (MQ=255)
                                     ttttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTaa                                                              <  1:4042406/41‑1 (MQ=255)
                                     ttttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATataaataa                                 >  1:2049871/1‑70 (MQ=255)
                                     ttttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATataaataa                                 >  1:448166/1‑70 (MQ=255)
                                     ttttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                <  1:2975705/71‑1 (MQ=255)
                                     ttttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                <  1:2676961/71‑1 (MQ=255)
                                       ttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                >  1:1643217/1‑69 (MQ=255)
                                       ttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                >  1:1389589/1‑69 (MQ=255)
                                       ttATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                >  1:2470154/1‑69 (MQ=255)
                                       ttATCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAaataat                                >  1:1803001/1‑69 (MQ=255)
                                        tATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAaataat                                >  1:2386328/1‑68 (MQ=255)
                                         aTCACTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGa                            >  1:1344677/1‑71 (MQ=255)
                                           cACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACt                          >  1:51271/1‑71 (MQ=255)
                                            aCTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGata                                                          >  1:696353/1‑38 (MQ=255)
                                            aCTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTc                         <  1:1428669/71‑1 (MQ=255)
                                            aCTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTc                         <  1:740301/71‑1 (MQ=255)
                                             cTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCAc                                                  <  1:3033308/45‑1 (MQ=255)
                                             cTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCAc                                                  <  1:1062372/45‑1 (MQ=255)
                                             cTACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCa                        <  1:3306696/71‑1 (MQ=255)
                                              tACAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGcc                                                    <  1:3222670/42‑1 (MQ=255)
                                               aCAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACt                      >  1:2973097/1‑71 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCa                                                   >  1:451991/1‑41 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCa                                                   >  1:3635932/1‑41 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACtat                                               >  1:2177977/1‑45 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGata                                        >  1:1386190/1‑52 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGata                                        >  1:1619485/1‑52 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCAc                       <  1:2831702/69‑1 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCAc                       <  1:1242675/69‑1 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACtt                     >  1:2552011/1‑71 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACtt                     >  1:538782/1‑71 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACtt                     >  1:918127/1‑71 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCAcct                     >  1:476675/1‑69 (MQ=255)
                                                cAGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCAcct                     >  1:1077330/1‑69 (MQ=255)
                                                 aGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTg                    >  1:1587684/1‑71 (MQ=255)
                                                 aGCTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACt                                                 >  1:1348933/1‑42 (MQ=255)
                                                   cTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGc                   >  1:3968793/1‑70 (MQ=255)
                                                   cTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGc                   >  1:2919862/1‑70 (MQ=255)
                                                   cTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGc                   >  1:2924507/1‑70 (MQ=255)
                                                   cTTATAGAGGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAAtaataa                              <  1:2763269/59‑1 (MQ=255)
                                                     tataGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCa                        >  1:2424803/1‑63 (MQ=255)
                                                      ataGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCaa                 <  1:3243392/69‑1 (MQ=255)
                                                        agagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAAc               >  1:3257545/1‑69 (MQ=255)
                                                        agagGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGa                                          >  1:1204386/1‑42 (MQ=255)
                                                        agagGCTTAAGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCAc                       >  1:3289423/1‑61 (MQ=255)
                                                         gagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCa             <  1:4101527/70‑1 (MQ=255)
                                                         gagGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCa             <  1:649174/70‑1 (MQ=255)
                                                          agGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGAtata                                      <  1:2797572/44‑1 (MQ=255)
                                                            gCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGa           <  1:3395544/69‑1 (MQ=255)
                                                              ttAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACgg        >  1:1554100/1‑70 (MQ=255)
                                                                aaGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGt       <  1:2355925/69‑1 (MQ=255)
                                                                aaGGAAATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACCTGTAAACCAGACGGt       <  1:161199/69‑1 (MQ=255)
                                                                   gagaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAAt    <  1:1037036/69‑1 (MQ=255)
                                                                    agaTTCGTAAGATATCAGCCACTATACCGATATAAAtaataa                              <  1:3456616/42‑1 (MQ=255)
                                                                    agaTTCGTAAGATATCAGCCACTATACCGATATAAAtaataa                              <  1:2463906/42‑1 (MQ=255)
                                                                    agaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAAttt  >  1:996181/1‑70 (MQ=255)
                                                                     gaTTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCa             >  1:4003622/1‑58 (MQ=255)
                                                                     |                                                                    
GATAAATCATTTGGTTGGCGCACTGGAATAATCTGAGTTTTATCACTACAGCTTATAGAGGCTTAAGGAGATTCGTAAGATATCAGCCACTATACCGATATAAATAATAAGACTCACTTGCAAACCAGACGGTAATTT  >  minE/2572872‑2573009

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: