Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,708 A→C 17.5% intergenic (+1476/‑303) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7080AC17.5% 210.8 / 58.9 114intergenic (+1476/‑303)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/54);  new base C (14/6);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 2.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTA  >  minE/2572639‑2572777
                                                                     |                                                                     
cAGGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctac                                                                      <  1:2688241/71‑1 (MQ=255)
         aaTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCacctac                                                                      <  1:3029301/62‑1 (MQ=255)
           tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCtt                                                                   >  1:2499087/1‑63 (MQ=255)
           tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                            >  1:1021027/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:4001960/71‑1 (MQ=255)
                      tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                <  1:734034/71‑1 (MQ=255)
                      tCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                <  1:3901341/71‑1 (MQ=255)
                       cGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                 >  1:3383275/1‑69 (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)
                            ccAGAAATCGGGGTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCCGAGCCTGAgtgt                                           <  1:3865189/70‑1 (MQ=255)
                             cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         <  1:567787/71‑1 (MQ=255)
                             cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         <  1:390349/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)
                                 aaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                           <  1:1237055/65‑1 (MQ=255)
                                   tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                         >  1:1624669/1‑49 (MQ=255)
                                   tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         >  1:1477177/1‑65 (MQ=255)
                                    cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAta                                         >  1:1617017/1‑64 (MQ=255)
                                     gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                         >  1:702095/1‑47 (MQ=255)
                                     gggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                 >  1:3634672/1‑71 (MQ=255)
                                     ggattCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGccc                                 >  1:3595513/4‑71 (MQ=255)
                                       gTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:2946995/1‑47 (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:3488594/1‑71 (MQ=255)
                                          cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                            >  1:3689940/1‑71 (MQ=255)
                                          cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                            >  1:2362725/1‑71 (MQ=255)
                                           aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                            <  1:3629112/70‑2 (MQ=255)
                                           aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagca                            <  1:3263023/70‑2 (MQ=255)
                                            aCGGATACTGATACATATTCCATCCACCTCTCAACCCAATCgg                                                      >  1:2531630/1‑43 (MQ=255)
                                            aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                         <  1:3335460/40‑1 (MQ=255)
                                             cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                         <  1:1144066/71‑1 (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)
                                                aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         >  1:2921150/1‑52 (MQ=255)
                                                aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                       >  1:1447469/1‑70 (MQ=255)
                                                 tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         >  1:86929/1‑51 (MQ=255)
                                                 tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                      >  1:1461899/1‑70 (MQ=255)
                                                   cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAc                                      >  1:2670562/1‑52 (MQ=255)
                                                   cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAAt                   <  1:3380212/71‑1 (MQ=255)
                                                    tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                <  1:328303/41‑1 (MQ=255)
                                                    tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         <  1:2543605/48‑1 (MQ=255)
                                                     gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgtgat                                          >  1:3220834/1‑44 (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)
                                                      aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                           <  1:3080931/60‑1 (MQ=255)
                                                      aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                           <  1:156782/60‑1 (MQ=255)
                                                      aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                       <  1:3131141/64‑1 (MQ=255)
                                                      aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                       <  1:2890094/64‑1 (MQ=255)
                                                      aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                <  1:1950267/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)
                                                             ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata         <  1:351529/71‑1 (MQ=255)
                                                             ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata         <  1:1348121/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: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: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)
                                                                     |                                                                     
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: