Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,344:1 +A 24.1% intergenic (+112/‑1667) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3441.A24.1% 302.4 / 165.9 141intergenic (+112/‑1667)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (45/62);  new base A (20/14);  total (65/76)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 155 total aligned reads displayed.

CTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTC  >  minE/2571274‑2571413
                                                                        |                                                                      
   attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     >  1:3597404/1‑70 (MQ=255)
   attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:1810774/70‑1 (MQ=255)
     tattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                                  >  1:1101002/1‑71 (MQ=255)
      attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                                  >  1:2657581/1‑70 (MQ=255)
       tttttGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                                  >  1:3723745/1‑69 (MQ=255)
             cATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                           <  1:4267138/70‑1 (MQ=255)
               tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                                                        >  1:2352113/1‑71 (MQ=255)
                atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:3364244/57‑1 (MQ=255)
                atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:3553326/57‑1 (MQ=255)
                atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                                                   >  1:4371958/1‑59 (MQ=255)
                atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                                       >  1:4021442/1‑71 (MQ=255)
                 ttACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCtctc                                                      <  1:2285486/71‑1 (MQ=255)
                    cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTTTGTTAATAAGTTc                                                                         <  1:687733/51‑1 (MQ=255)
                    cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                                         <  1:4607018/51‑1 (MQ=255)
                    cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTgca                                                                   >  1:5408570/1‑53 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:284437/52‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:5198668/52‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:1216093/52‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATg                                                                 <  1:5365132/56‑1 (MQ=255)
                       cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGa                                                >  1:2016698/1‑71 (MQ=255)
                        aTTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAt                                               <  1:327776/71‑1 (MQ=255)
                         ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                                              >  1:598891/1‑71 (MQ=255)
                           aGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAttt                                             <  1:2451652/70‑1 (MQ=255)
                              ttCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                                                >  1:4025417/1‑48 (MQ=255)
                                tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGAtttctt                                          <  1:4252406/69‑1 (MQ=255)
                                tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCTgg                                                                    <  1:4008848/42‑1 (MQ=16)
                                cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                                         <  1:4416152/39‑1 (MQ=255)
                                cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                                                       >  1:4491163/1‑40 (MQ=255)
                                cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                                                       >  1:3069491/1‑40 (MQ=255)
                                 gatgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                      <  1:1022779/71‑1 (MQ=255)
                                  atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTc                                                                         <  1:3842416/38‑1 (MQ=12)
                                  atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATc                                       >  1:3425292/1‑71 (MQ=255)
                                  atgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCAtt                                                   <  1:1042987/57‑1 (MQ=255)
                                    gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     >  1:4875433/1‑37 (MQ=255)
                                    gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                      <  1:3671395/68‑1 (MQ=255)
                                    gatTTTACTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                                   >  1:734041/1‑71 (MQ=255)
                                    gatATTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     >  1:5530527/1‑37 (MQ=255)
                                     tttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTc                                   >  1:5093691/2‑71 (MQ=255)
                                     atTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                                                   >  1:3172798/1‑38 (MQ=255)
                                     atTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCtctc                                                      <  1:2802838/51‑1 (MQ=255)
                                     atTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                                   <  1:2984344/70‑1 (MQ=255)
                                      tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                           <  1:720328/45‑1 (MQ=255)
                                      tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt                                 <  1:119029/71‑1 (MQ=255)
                                      ttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                                  >  1:1490904/1‑71 (MQ=255)
                                       tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCtt                                 >  1:2110719/1‑71 (MQ=255)
                                        tttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                <  1:3487164/70‑1 (MQ=255)
                                        tttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                <  1:4212067/70‑1 (MQ=255)
                                        tttCTTTAATATTCATTATGTTAATAAGTTCTCGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt                               >  1:5527992/1‑71 (MQ=255)
                                         ttccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCtt                                 <  1:1754179/67‑1 (MQ=255)
                                         ttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                                         >  1:2233766/1‑60 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:529762/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:4726710/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:5635651/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                                <  1:2759612/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                                <  1:1522883/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                                <  1:4309757/68‑1 (MQ=255)
                                          tCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                              <  1:5695919/70‑1 (MQ=255)
                                            cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             <  1:116806/71‑1 (MQ=255)
                                            cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             >  1:4842888/1‑71 (MQ=255)
                                             tttAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTc                             <  1:475903/69‑1 (MQ=255)
                                               tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGcc                         >  1:837254/1‑71 (MQ=255)
                                               tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGGTCTGcc                         >  1:5014420/1‑71 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                                         >  1:2221760/1‑56 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                                         >  1:4814237/1‑56 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCt                            >  1:3759709/1‑69 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTgt                          >  1:5697330/1‑70 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:2674080/1‑71 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:2356605/1‑71 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:3672754/1‑71 (MQ=255)
                                                aaTATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                            <  1:2164482/67‑1 (MQ=255)
                                                aaTATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTgat                              >  1:2205858/1‑64 (MQ=255)
                                                 atatTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCg                        <  1:1877765/70‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATc                                       >  1:5095556/1‑52 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                            <  1:2512567/63‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                            <  1:1716614/63‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                            <  1:1611576/63‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAgt                    <  1:5305193/71‑2 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                     <  1:3328060/70‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                     <  1:1864496/70‑1 (MQ=255)
                                                      cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                  >  1:1549294/1‑71 (MQ=255)
                                                       agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             >  1:5207379/4‑60 (MQ=255)
                                                       agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             >  1:1731674/4‑60 (MQ=255)
                                                       aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                      >  1:482338/1‑50 (MQ=255)
                                                       aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATg                 <  1:4985725/71‑1 (MQ=255)
                                                       aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTATATCCGTCTTGTTCTGCCGCAGGAt                  <  1:5026670/70‑1 (MQ=255)
                                                        ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                                         <  1:398453/46‑1 (MQ=255)
                                                        ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                              <  1:127074/57‑1 (MQ=255)
                                                        ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                              <  1:4167997/57‑1 (MQ=255)
                                                           tGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc             <  1:4784700/71‑1 (MQ=255)
                                                           tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:4295897/53‑1 (MQ=255)
                                                           tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:4290170/53‑1 (MQ=255)
                                                             ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt            <  1:1536982/70‑1 (MQ=255)
                                                              taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                                         >  1:414453/1‑40 (MQ=255)
                                                              taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTg                           >  1:4514969/1‑54 (MQ=255)
                                                              taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTg                           >  1:4344942/1‑54 (MQ=255)
                                                              taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta           >  1:4223154/1‑70 (MQ=255)
                                                                ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCg                                     <  1:4732441/42‑1 (MQ=255)
                                                                ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt                               <  1:4440026/48‑1 (MQ=255)
                                                                ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata         >  1:2060905/1‑70 (MQ=255)
                                                                    gTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTc                             <  1:1266936/46‑1 (MQ=255)
                                                                        |                                                                      
CTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTC  >  minE/2571274‑2571413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: