Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,338 T→G 30.5% intergenic (+106/‑1673) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3380TG30.5% 133.2 / 129.9 141intergenic (+106/‑1673)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/59);  new base G (21/22);  total (60/81)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

Only 100 of 156 total aligned reads displayed.

CTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATA  >  minE/2571274‑2571406
                                                                 |                                                                      
cTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                                  >  1:901165/1‑71 (MQ=255)
  aaTTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              >  1:2212674/1‑71 (MQ=255)
   attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTtaat                                                                        <  1:1100764/62‑1 (MQ=255)
   attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              >  1:3597404/1‑70 (MQ=255)
   attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              <  1:4295718/70‑1 (MQ=255)
    ttattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                                                    <  1:4140373/65‑1 (MQ=255)
    ttattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                                                    <  1:1623633/65‑1 (MQ=255)
     tattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                           >  1:1101002/1‑71 (MQ=255)
      attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                           >  1:2657581/1‑70 (MQ=255)
      attTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                           >  1:2151432/1‑70 (MQ=255)
       tttttGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                           >  1:3723745/1‑69 (MQ=255)
            gCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                    >  1:764500/1‑71 (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)
                atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                                            >  1:4377557/1‑59 (MQ=255)
                atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                                >  1:4021442/1‑71 (MQ=255)
                 ttACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTtaataa                                                                      >  1:4959213/1‑50 (MQ=255)
                    cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTTTGTTAATAAGTTc                                                                  <  1:687733/51‑1 (MQ=255)
                    cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                                  <  1:4607018/51‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                                                    <  1:5675227/48‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              <  1:284437/52‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              <  1:692641/52‑1 (MQ=255)
                     ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATg                                                          <  1:5365132/56‑1 (MQ=255)
                       cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                                         <  1:1926925/55‑1 (MQ=255)
                        aTTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAt                                        <  1:327776/71‑1 (MQ=255)
                         ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                           >  1:1240828/1‑51 (MQ=255)
                          tAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGTAtt                                       <  1:225473/70‑1 (MQ=255)
                           aGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAttt                                      <  1:2451652/70‑1 (MQ=255)
                              ttCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                                         >  1:4025417/1‑48 (MQ=255)
                                tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGAtttctt                                   <  1:4252406/69‑1 (MQ=255)
                                tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCt                                                                <  1:4702250/40‑1 (MQ=12)
                                tgatgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCTgg                                                             <  1:4008848/42‑1 (MQ=16)
                                cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                                                >  1:3069491/1‑40 (MQ=255)
                                cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                                                >  1:4491163/1‑40 (MQ=255)
                                 gatgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                               <  1:1022779/71‑1 (MQ=255)
                                  atgatTTTTCCTTTAATATCCAGCATGTTaagaag                                                                     >  1:4082362/1‑35 (MQ=9)
                                  atgatTTTTCCTTTAATATCCAGCATGTTaagaag                                                                     >  1:1241167/1‑35 (MQ=9)
                                  atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTc                                                                  <  1:3842416/38‑1 (MQ=12)
                                  atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATc                                >  1:3425292/1‑71 (MQ=255)
                                  atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCTgg                                                             >  1:1356806/1‑41 (MQ=16)
                                  atgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCAtt                                            <  1:1042987/57‑1 (MQ=255)
                                    gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              >  1:4875433/1‑37 (MQ=255)
                                    gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                               <  1:3671395/68‑1 (MQ=255)
                                    gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                            >  1:3264253/1‑71 (MQ=255)
                                    gatTTTACTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                            >  1:734041/1‑71 (MQ=255)
                                    gatATTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                              >  1:5530527/1‑37 (MQ=255)
                                      tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                    <  1:1745162/45‑1 (MQ=255)
                                      tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt                          <  1:119029/71‑1 (MQ=255)
                                       tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                           >  1:1847902/1‑70 (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:2353435/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                         <  1:2989061/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                         <  1:5635651/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                         <  1:2033084/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                         <  1:2759612/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                         <  1:4254032/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                         <  1:1522883/68‑1 (MQ=255)
                                          tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                         <  1:4309757/68‑1 (MQ=255)
                                            cTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt                       <  1:5031128/69‑1 (MQ=255)
                                            cTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                     >  1:4191709/1‑71 (MQ=255)
                                            cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                      >  1:4842888/1‑71 (MQ=255)
                                            cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                      >  1:1524999/1‑71 (MQ=255)
                                             tttAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                         <  1:5480766/67‑1 (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:3672754/1‑71 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                   >  1:4743825/1‑71 (MQ=255)
                                               tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                   >  1:2674080/1‑71 (MQ=255)
                                                aaTATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                     <  1:2164482/67‑1 (MQ=255)
                                                 atatTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCg                 <  1:1877765/70‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATc                                >  1:5095556/1‑52 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                     <  1:1611576/63‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                     <  1:2512567/63‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                     <  1:1716614/63‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAgt             <  1:5305193/71‑2 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg              <  1:3328060/70‑1 (MQ=255)
                                                    ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg              <  1:1175257/70‑1 (MQ=255)
                                                      cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt           >  1:1549294/1‑71 (MQ=255)
                                                      cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt           >  1:2803109/1‑71 (MQ=255)
                                                       agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                      >  1:5207379/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:4167997/57‑1 (MQ=255)
                                                         tATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt           >  1:5231582/1‑68 (MQ=255)
                                                           tGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc      <  1:4784700/71‑1 (MQ=255)
                                                           tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                         <  1:4290170/53‑1 (MQ=255)
                                                             ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt     <  1:596909/70‑1 (MQ=255)
                                                              taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTg                    >  1:4514969/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:4326404/70‑1 (MQ=255)
                                                                ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata  >  1:2060905/1‑70 (MQ=255)
                                                                 |                                                                      
CTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATA  >  minE/2571274‑2571406

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: