Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,713 C→T 16.4% intergenic (+1481/‑298) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7130CT16.4% 329.7 / 47.5 128intergenic (+1481/‑298)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (48/59);  new base T (14/7);  total (62/66)
Fisher's exact test for biased strand distribution p-value = 9.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 131 total aligned reads displayed.

AATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGT  >  minE/2572648‑2572781
                                                                 |                                                                    
aaTCCCAAGGTGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                                                                 >  1:1612121/1‑71 (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)
    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)
      aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                           >  1:2666170/1‑71 (MQ=255)
         gggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAg                                                        >  1:1632056/1‑71 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                     >  1:694563/1‑70 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                    <  1:734034/71‑1 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                    <  1:4001960/71‑1 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                    <  1:3901341/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:779791/70‑1 (MQ=255)
                   ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                               <  1:3787628/70‑1 (MQ=255)
                    cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                             <  1:390349/71‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                            <  1:131062/71‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                            <  1:914487/71‑1 (MQ=255)
                      gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                           >  1:3445080/1‑71 (MQ=255)
                          tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                             >  1:1624669/1‑49 (MQ=255)
                          tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                             >  1:1477177/1‑65 (MQ=255)
                           cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGccc                                                                >  1:3972475/1‑45 (MQ=255)
                            gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                     >  1:3634672/1‑71 (MQ=255)
                            ggattCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGccc                                     >  1:3595513/4‑71 (MQ=255)
                              gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                                   <  1:782865/71‑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:3629112/70‑2 (MQ=255)
                                   aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                             <  1:3335460/40‑1 (MQ=255)
                                   aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGtt                                              <  1:2544727/55‑1 (MQ=255)
                                    cGGATACTGATAAAAATTCCACCTCCCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                             >  1:2784264/1‑71 (MQ=255)
                                    cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                             >  1:1169944/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)
                                       aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGATCCTGAGTGTTATACCGCCCCTGCTGCGGtt                          <  1:2890548/71‑1 (MQ=255)
                                       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)
                                        tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                          >  1:1461899/1‑70 (MQ=255)
                                         aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGa                                                         <  1:1864114/38‑1 (MQ=12)
                                          cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                          >  1:2670562/1‑52 (MQ=255)
                                          cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                        <  1:364774/70‑1 (MQ=255)
                                          cTGATAAAAATCCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                          >  1:2403125/1‑52 (MQ=255)
                                           tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                    <  1:328303/41‑1 (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:2645672/1‑50 (MQ=255)
                                            gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                          >  1:346663/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)
                                             aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                    <  1:1950267/71‑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)
                                                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)
                                                    ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGat              >  1:3399694/1‑70 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata          <  1:1199011/71‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata          <  1:612605/71‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat           <  1:691148/70‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat           <  1:1454696/70‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat           <  1:686262/70‑1 (MQ=255)
                                                         cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt        <  1:723006/71‑1 (MQ=255)
                                                          ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                       <  1:1095040/55‑1 (MQ=255)
                                                          ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt        >  1:1540990/1‑70 (MQ=255)
                                                          ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt        >  1:1536064/1‑70 (MQ=255)
                                                           tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGATAATCTTACCGATAATAGTta      >  1:1392362/1‑71 (MQ=255)
                                                            acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta      >  1:2413933/1‑70 (MQ=255)
                                                            acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta      >  1:2032962/1‑70 (MQ=255)
                                                             cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta      <  1:1552513/69‑1 (MQ=255)
                                                             cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta      <  1:1166089/69‑1 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                        >  1:2355988/1‑50 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg   <  1:2932678/71‑1 (MQ=255)
                                                               ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                       >  1:8035/1‑50 (MQ=255)
                                                                 ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                               >  1:2983549/1‑40 (MQ=255)
                                                                 ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                    >  1:1564481/1‑51 (MQ=255)
                                                                 ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt  <  1:650788/69‑1 (MQ=255)
                                                                 |                                                                    
AATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGT  >  minE/2572648‑2572781

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: