Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,345 T→G 23.5% intergenic (+113/‑1666) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3450TG23.5% 205.4 / 96.6 145intergenic (+113/‑1666)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/65);  new base G (20/14);  total (66/79)
Fisher's exact test for biased strand distribution p-value = 8.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 150 total aligned reads displayed.

AATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTC  >  minE/2571276‑2571413
                                                                       |                                                                     
aaTTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     >  1:2212674/1‑71 (MQ=255)
 attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:4483016/70‑1 (MQ=255)
 attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:4295718/70‑1 (MQ=255)
 attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:1810774/70‑1 (MQ=255)
   tattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGAt                                                                  >  1:1101002/1‑71 (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:2722097/1‑71 (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)
              atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                                       >  1:5669681/1‑71 (MQ=255)
               ttACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCtctc                                                      <  1:2285486/71‑1 (MQ=255)
                  cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTgca                                                                   >  1:5408570/1‑53 (MQ=255)
                   ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                                     <  1:692641/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)
                     cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                                                <  1:1926925/55‑1 (MQ=255)
                       ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                                              >  1:4836689/1‑71 (MQ=255)
                        tAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGTAtt                                              <  1:225473/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:4491163/1‑40 (MQ=255)
                               gatgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                      <  1:1022779/71‑1 (MQ=255)
                                atgatTTTTCCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATc                                       >  1:3425292/1‑71 (MQ=255)
                                atgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCAtt                                                   <  1:1042987/57‑1 (MQ=255)
                                  gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                      <  1:3671395/68‑1 (MQ=255)
                                  gatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                                   >  1:3264253/1‑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)
                                    tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                                           <  1:1745162/45‑1 (MQ=255)
                                    tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt                                 <  1:119029/71‑1 (MQ=255)
                                    ttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                                  >  1:1490904/1‑71 (MQ=255)
                                     tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                                  >  1:1847902/1‑70 (MQ=255)
                                      tttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                <  1:4212067/70‑1 (MQ=255)
                                      tttCTTTAATATTCATTATGTTAATAAGTTCTCGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt                               >  1:5527992/1‑71 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:2353435/68‑1 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:4726710/68‑1 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:327799/68‑1 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:2989061/68‑1 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:529762/68‑1 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:5635651/68‑1 (MQ=255)
                                        tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:2033084/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)
                                          cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             <  1:116806/71‑1 (MQ=255)
                                          cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             >  1:4842888/1‑71 (MQ=255)
                                           tttAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTc                             <  1:475903/69‑1 (MQ=255)
                                             tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGGTCTGcc                         >  1:5014420/1‑71 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                                         >  1:4814237/1‑56 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                                         >  1:2221760/1‑56 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTgt                          >  1:5697330/1‑70 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:5528958/1‑71 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:2674080/1‑71 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:2356605/1‑71 (MQ=255)
                                             tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc                          >  1:4743825/1‑71 (MQ=255)
                                              aaTATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                            <  1:2164482/67‑1 (MQ=255)
                                              aaTATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTgat                              >  1:3672967/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:1611576/63‑1 (MQ=255)
                                                  ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAgt                    <  1:5305193/71‑2 (MQ=255)
                                                  ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                     <  1:3328060/70‑1 (MQ=255)
                                                    cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                  >  1:2803109/1‑71 (MQ=255)
                                                     agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             >  1:1731674/4‑60 (MQ=255)
                                                     agcaTGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                             >  1:5207379/4‑60 (MQ=255)
                                                     aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                  <  1:3930766/70‑1 (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)
                                                         tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:380857/53‑1 (MQ=255)
                                                         tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                                <  1:4290170/53‑1 (MQ=255)
                                                           ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                                      <  1:4938778/44‑1 (MQ=255)
                                                           ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta           <  1:1616543/71‑1 (MQ=255)
                                                           ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt            <  1:1536982/70‑1 (MQ=255)
                                                           ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt            <  1:596909/70‑1 (MQ=255)
                                                           ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt            <  1:3411049/70‑1 (MQ=255)
                                                            taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                                         >  1:414453/1‑40 (MQ=255)
                                                            taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCt                                  >  1:479064/1‑47 (MQ=255)
                                                            taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                            >  1:3661204/1‑53 (MQ=255)
                                                            taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta           >  1:4223154/1‑70 (MQ=255)
                                                              ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata         >  1:2060905/1‑70 (MQ=255)
                                                              ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata         <  1:4326404/70‑1 (MQ=255)
                                                                aaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                >  1:246516/1‑61 (MQ=255)
                                                                  gTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTc                             <  1:1266936/46‑1 (MQ=255)
                                                                  gTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc             <  1:153790/62‑1 (MQ=255)
                                                                     cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc  <  1:3511341/70‑1 (MQ=255)
                                                                       tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc  <  1:774512/69‑1 (MQ=255)
                                                                       |                                                                     
AATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTC  >  minE/2571276‑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: