Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,705 C→T 18.4% intergenic (+1473/‑306) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7050CT18.4% 271.6 / 49.3 114intergenic (+1473/‑306)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/53);  new base T (14/7);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 120 total aligned reads displayed.

TAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGT  >  minE/2572635‑2572775
                                                                      |                                                                      
taAACAGGTGGCTAATCCCAAGGGGATCTATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCac                                                                        >  1:1816853/1‑71 (MQ=255)
    cAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctac                                                                    <  1:2688241/71‑1 (MQ=255)
             aaTCCCAAGGTGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                                                           >  1:1612121/1‑71 (MQ=255)
             aaTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctac                                                                    <  1:3029301/62‑1 (MQ=255)
               tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCtt                                                                 >  1:2499087/1‑63 (MQ=255)
               tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                          >  1:1021027/1‑70 (MQ=255)
               tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                          >  1:1269454/1‑70 (MQ=255)
               tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                         <  1:2772588/71‑1 (MQ=255)
               tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                         <  1:2158085/71‑1 (MQ=255)
               tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAAATCAACCTACCTTCCAGCCCa                                                         <  1:3674630/71‑1 (MQ=255)
                 ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGGCCAAt                                                       >  1:3284214/1‑71 (MQ=255)
                 ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                        >  1:1063436/1‑70 (MQ=255)
                 ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                       >  1:2753554/1‑71 (MQ=255)
                  cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                         <  1:933206/68‑1 (MQ=255)
                   aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctacct                                                                  <  1:1733727/58‑1 (MQ=255)
                   aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                     >  1:2666170/1‑71 (MQ=255)
                     ggggATCGATTCCAGAAATCGGGTGCAACGGATACTGATAAAAATTCCacc                                                                       <  1:3147181/51‑1 (MQ=255)
                      gggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAg                                                  >  1:1632056/1‑71 (MQ=255)
                          tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                              <  1:3901341/71‑1 (MQ=255)
                          tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                              <  1:734034/71‑1 (MQ=255)
                           cGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                               >  1:3383275/1‑69 (MQ=255)
                            gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                              >  1:716211/1‑53 (MQ=255)
                            gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAg                                            >  1:1264899/1‑71 (MQ=255)
                             aTTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                           <  1:3702400/71‑1 (MQ=255)
                              ttCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                          <  1:874905/55‑1 (MQ=255)
                               tCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                          <  1:605181/70‑1 (MQ=255)
                                ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                         <  1:3787628/70‑1 (MQ=255)
                                 cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                       <  1:390349/71‑1 (MQ=255)
                                  aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                       <  1:1872515/54‑1 (MQ=255)
                                  aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                      <  1:914487/71‑1 (MQ=255)
                                   gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                     >  1:3445080/1‑71 (MQ=255)
                                    aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                              >  1:2142621/1‑45 (MQ=255)
                                     aaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                         <  1:1237055/65‑1 (MQ=255)
                                        cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGccc                                                          >  1:3972475/1‑45 (MQ=255)
                                        cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                       >  1:1617017/1‑64 (MQ=255)
                                         gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                       >  1:702095/1‑47 (MQ=255)
                                         gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                               >  1:3634672/1‑71 (MQ=255)
                                           gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                     >  1:2946995/1‑47 (MQ=255)
                                           gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                             <  1:782865/71‑1 (MQ=255)
                                            ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                               <  1:1275520/52‑1 (MQ=255)
                                            ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                         <  1:3534905/58‑1 (MQ=255)
                                              cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                           <  1:3681761/70‑1 (MQ=255)
                                              cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                            >  1:1241624/1‑67 (MQ=255)
                                              cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                          >  1:3689940/1‑71 (MQ=255)
                                              cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                          >  1:3488594/1‑71 (MQ=255)
                                               aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                          <  1:3263023/70‑2 (MQ=255)
                                               aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                          <  1:3629112/70‑2 (MQ=255)
                                                aCGGATACTGATACATATTCCATCCACCTCTCAACCCAATCgg                                                    >  1:2531630/1‑43 (MQ=255)
                                                aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                       <  1:3335460/40‑1 (MQ=255)
                                                 cGGATACTGATAAAAATTCCACCTCCCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                       >  1:2784264/1‑71 (MQ=255)
                                                 cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                       >  1:1169944/1‑71 (MQ=255)
                                                 cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                       <  1:1144066/71‑1 (MQ=255)
                                                  ggATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                      >  1:1446459/1‑71 (MQ=255)
                                                  ggATACTGATAAAAATTCCACCTACCTTCCAGccc                                                          >  1:1016167/1‑35 (MQ=255)
                                                  ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                     >  1:2394455/1‑40 (MQ=255)
                                                   gATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCg                                                     <  1:545833/39‑1 (MQ=25)
                                                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtg                                          <  1:3245108/49‑1 (MQ=255)
                                                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                       >  1:2921150/1‑52 (MQ=255)
                                                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                     >  1:532667/1‑70 (MQ=255)
                                                    aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                     >  1:1447469/1‑70 (MQ=255)
                                                     tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                       >  1:86929/1‑51 (MQ=255)
                                                     tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                    >  1:1461899/1‑70 (MQ=255)
                                                      aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGa                                                   <  1:1864114/38‑1 (MQ=12)
                                                       cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                    >  1:2670562/1‑52 (MQ=255)
                                                       cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                 <  1:3380212/71‑1 (MQ=255)
                                                       cTGATAAAAATCCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                    >  1:2403125/1‑52 (MQ=255)
                                                        tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                       <  1:2543605/48‑1 (MQ=255)
                                                        tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                    <  1:3353523/67‑1 (MQ=255)
                                                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgtgat                                        >  1:3220834/1‑44 (MQ=255)
                                                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                    >  1:1525899/1‑50 (MQ=255)
                                                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                    >  1:346663/1‑50 (MQ=255)
                                                         gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                    >  1:2645672/1‑50 (MQ=255)
                                                         gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                        <  1:3902469/46‑3 (MQ=255)
                                                         gATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                     >  1:2062969/1‑49 (MQ=255)
                                                         gATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                    <  1:504292/66‑1 (MQ=255)
                                                          aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                        >  1:4064961/1‑43 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGATTGTTATACCGCCCCTGCTGCGGTTAATGtt              <  1:2329000/71‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                       <  1:2963637/46‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                         <  1:156782/60‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                         <  1:3080931/60‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                     <  1:2890094/64‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt              <  1:3621980/71‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt               <  1:3625784/70‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGt               <  1:2755390/70‑1 (MQ=255)
                                                          aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTAGTGTTATACCGCCCCTGCTGCGGTTAATGtt              <  1:3137623/70‑1 (MQ=255)
                                                           tAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                             >  1:2141774/1‑55 (MQ=255)
                                                           tAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTa             <  1:3945637/71‑1 (MQ=255)
                                                            aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAg                                            >  1:3628962/1‑39 (MQ=255)
                                                            aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                     >  1:763830/1‑62 (MQ=255)
                                                             aaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                       >  1:1564126/1‑43 (MQ=255)
                                                             aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                         <  1:610657/57‑1 (MQ=255)
                                                               aaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa         >  1:2694211/1‑71 (MQ=255)
                                                                 ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata       <  1:1348121/71‑1 (MQ=255)
                                                                 ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata       <  1:351529/71‑1 (MQ=255)
                                                                   ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa      >  1:429961/1‑70 (MQ=255)
                                                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata    <  1:612605/71‑1 (MQ=255)
                                                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat     <  1:686262/70‑1 (MQ=255)
                                                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat     <  1:691148/70‑1 (MQ=255)
                                                                    cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat     <  1:1454696/70‑1 (MQ=255)
                                                                      cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt  <  1:723006/71‑1 (MQ=255)
                                                                      |                                                                      
TAAACAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGT  >  minE/2572635‑2572775

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: