Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,711 T→C 13.3% intergenic (+1479/‑300) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7110TC13.3% 299.6 / 33.8 113intergenic (+1479/‑300)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/55);  new base C (10/5);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 130 total aligned reads displayed.

AATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAG  >  minE/2572648‑2572780
                                                               |                                                                     
aaTCCCAAGGTGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                                                                >  1:1612121/1‑71 (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)
    ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                             >  1:1063436/1‑70 (MQ=255)
     cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                              <  1:933206/68‑1 (MQ=255)
      aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctacct                                                                       <  1:1733727/58‑1 (MQ=255)
      aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                          >  1:2666170/1‑71 (MQ=255)
       aGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                           >  1:2702818/1‑69 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                   <  1:3901341/71‑1 (MQ=255)
             tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                   <  1:734034/71‑1 (MQ=255)
               gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                   >  1:716211/1‑53 (MQ=255)
               gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAg                                                 >  1:1264899/1‑71 (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)
                   ccAGAAATCGGGGTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCCGAGCCTGAgtgt                                              <  1:3865189/70‑1 (MQ=255)
                    cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                            <  1:390349/71‑1 (MQ=255)
                    cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                            <  1:567787/71‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                            <  1:1872515/54‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                           <  1:914487/71‑1 (MQ=255)
                     aGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtat                                           <  1:131062/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)
                           cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                            >  1:1617017/1‑64 (MQ=255)
                            gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                            >  1:702095/1‑47 (MQ=255)
                            ggattCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGccc                                    >  1:3595513/4‑71 (MQ=255)
                              gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCct                                  <  1:782865/71‑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:2362725/1‑71 (MQ=255)
                                 cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                               >  1:3689940/1‑71 (MQ=255)
                                  aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                               <  1:3263023/70‑2 (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)
                                     ggATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                           >  1:1446459/1‑71 (MQ=255)
                                     ggATACTGATAAAAATTCCACCTACCTTCCAGccc                                                               >  1:1016167/1‑35 (MQ=255)
                                       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)
                                        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)
                                          cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                      <  1:3380212/71‑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)
                                            gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                         >  1:2645672/1‑50 (MQ=255)
                                            gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                         >  1:346663/1‑50 (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:3080931/60‑1 (MQ=255)
                                             aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                              <  1:156782/60‑1 (MQ=255)
                                             aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                          <  1:3131141/64‑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)
                                                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)
                                                      ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa           >  1:429961/1‑70 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata         <  1:1199011/71‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat          <  1:1454696/70‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat          <  1:686262/70‑1 (MQ=255)
                                                       cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat          <  1:691148/70‑1 (MQ=255)
                                                         cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt       <  1:723006/71‑1 (MQ=255)
                                                          ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                      <  1:1095040/55‑1 (MQ=255)
                                                          ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt       >  1:1536064/1‑70 (MQ=255)
                                                          ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt       >  1:1540990/1‑70 (MQ=255)
                                                           tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGATAATCTTACCGATAATAGTta     >  1:1392362/1‑71 (MQ=255)
                                                            acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     >  1:2032962/1‑70 (MQ=255)
                                                            acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     >  1:2413933/1‑70 (MQ=255)
                                                             cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     <  1:1552513/69‑1 (MQ=255)
                                                             cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta     <  1:1166089/69‑1 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                       >  1:66270/1‑50 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                       >  1:2355988/1‑50 (MQ=255)
                                                              ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg  <  1:2932678/71‑1 (MQ=255)
                                                               ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                      >  1:3446/1‑50 (MQ=255)
                                                               ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                      >  1:8035/1‑50 (MQ=255)
                                                               |                                                                     
AATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAG  >  minE/2572648‑2572780

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: