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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7070TC18.4% 257.4 / 52.7 114intergenic (+1475/‑304)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/53);  new base C (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 = 5.12e-01

Only 100 of 122 total aligned reads displayed.

CAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTA  >  minE/2572639‑2572777
                                                                    |                                                                      
cAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctac                                                                      <  1:2688241/71‑1 (MQ=255)
         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)
             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)
                aGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                        >  1:2702818/1‑69 (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)
                        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)
                            ccAGAAATCGGGGTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCCGAGCCTGAgtgt                                           <  1:3865189/70‑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)
                                aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                >  1:2142621/1‑45 (MQ=255)
                                 aaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                           <  1:1237055/65‑1 (MQ=255)
                                   tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         >  1:1477177/1‑65 (MQ=255)
                                    cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGccc                                                            >  1:3972475/1‑45 (MQ=255)
                                    cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                         >  1:1617017/1‑64 (MQ=255)
                                     gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                 >  1:3634672/1‑71 (MQ=255)
                                     ggattCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGccc                                 >  1:3595513/4‑71 (MQ=255)
                                       gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:2946995/1‑47 (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:3488594/1‑71 (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: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)
                                                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)
                                                 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)
                                                   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)
                                                     gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgtgat                                          >  1:3220834/1‑44 (MQ=255)
                                                     gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                      >  1:2645672/1‑50 (MQ=255)
                                                     gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                      >  1:346663/1‑50 (MQ=255)
                                                     gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATAc                                      >  1:1525899/1‑50 (MQ=255)
                                                     gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                          <  1:3902469/46‑3 (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)
                                                      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)
                                                       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)
                                                             ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata         <  1:1348121/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: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)
                                                                   ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                   <  1:1095040/55‑1 (MQ=255)
                                                                   ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt    >  1:1536064/1‑70 (MQ=255)
                                                                    tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGATAATCTTACCGATAATAGTta  >  1:1392362/1‑71 (MQ=255)
                                                                    |                                                                      
CAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTA  >  minE/2572639‑2572777

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: