Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,329 T→G 29.2% intergenic (+97/‑1682) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3290TG29.2% 154.6 / 120.5 137intergenic (+97/‑1682)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/51);  new base G (21/19);  total (67/70)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 144 total aligned reads displayed.

TGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATG  >  minE/2571259‑2571398
                                                                       |                                                                       
  cTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCAtt                                                                        <  1:3430721/70‑1 (MQ=255)
     cTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGt                                                                    >  1:1917431/1‑71 (MQ=255)
       tGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGt                                                                    >  1:3885188/1‑69 (MQ=255)
       tGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGt                                                                    >  1:4428264/1‑69 (MQ=255)
        gTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGtt                                                                   <  1:5048383/69‑1 (MQ=255)
               cTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                          >  1:901165/1‑71 (MQ=255)
               cTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                          <  1:4524869/71‑1 (MQ=255)
                  attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATg                                                                     >  1:486374/1‑57 (MQ=255)
                  attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                      >  1:3597404/1‑70 (MQ=255)
                  attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                      <  1:4295718/70‑1 (MQ=255)
                  attattTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                      <  1:1810774/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)
                              tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                                         >  1:2722097/1‑71 (MQ=255)
                              tatTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctc                                         >  1:2352113/1‑71 (MQ=255)
                               atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                      <  1:3553326/57‑1 (MQ=255)
                               atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGa                                                    >  1:4377557/1‑59 (MQ=255)
                               atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                        >  1:5669681/1‑71 (MQ=255)
                               atTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGctct                                        >  1:4021442/1‑71 (MQ=255)
                                ttACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTtaataa                                                              >  1:4959213/1‑50 (MQ=255)
                                   cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTTTGTTAATAAGTTc                                                          <  1:687733/51‑1 (MQ=255)
                                   cTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTgca                                                    >  1:5408570/1‑53 (MQ=255)
                                    ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGt                                                            <  1:5675227/48‑1 (MQ=255)
                                    ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                      <  1:284437/52‑1 (MQ=255)
                                    ttCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTg                                                      <  1:1216093/52‑1 (MQ=255)
                                      cATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGa                                                 <  1:1926925/55‑1 (MQ=255)
                                       aTTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAt                                <  1:327776/71‑1 (MQ=255)
                                        ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                               >  1:4836689/1‑71 (MQ=255)
                                        ttAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGAtt                               >  1:598891/1‑71 (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)
                                               cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTc                                                          <  1:4416152/39‑1 (MQ=255)
                                               cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                                        >  1:3069491/1‑40 (MQ=255)
                                               cGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCt                                                        >  1:4491163/1‑40 (MQ=255)
                                                gatgatTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATcc                       <  1:1022779/71‑1 (MQ=255)
                                                 atgatTTTTCCTTTAATATCCAGCATGTTaagaag                                                             >  1:1241167/1‑35 (MQ=9)
                                                 atgatTTTTCCTTTAATATCCAGCATGTTaagaag                                                             >  1:4082362/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)
                                                   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)
                                                     tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGtgtg                                            <  1:1745162/45‑1 (MQ=255)
                                                     tttttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt                  <  1:119029/71‑1 (MQ=255)
                                                      tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCtt                  >  1:2110719/1‑71 (MQ=255)
                                                      tttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCt                   >  1:1847902/1‑70 (MQ=255)
                                                       tttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                 <  1:4212067/70‑1 (MQ=255)
                                                        ttccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCtt                  <  1:1754179/67‑1 (MQ=255)
                                                        ttCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTa                          >  1:2233766/1‑60 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                 <  1:2033084/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                 <  1:4726710/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                 <  1:2353435/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                 <  1:529762/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                 <  1:327799/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                 <  1:4309757/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                 <  1:4254032/68‑1 (MQ=255)
                                                         tccTTTAATATCCAGCATGTTAAGAAGTTCTGGGATAGCTGTAATCTCATTAGATTTCTCATCCGTCTTg                 <  1:1522883/68‑1 (MQ=255)
                                                         tCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGtt               <  1:5695919/70‑1 (MQ=255)
                                                           cTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt             >  1:4191709/1‑71 (MQ=255)
                                                           cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc              >  1:1524999/1‑71 (MQ=255)
                                                           cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc              <  1:116806/71‑1 (MQ=255)
                                                           cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc              >  1:4842888/1‑71 (MQ=255)
                                                            tttAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTc              <  1:475903/69‑1 (MQ=255)
                                                            tttAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTg                 <  1:5480766/67‑1 (MQ=255)
                                                              tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGcc          >  1:109743/1‑71 (MQ=255)
                                                              tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGcc          >  1:837254/1‑71 (MQ=255)
                                                              tAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGGTCTGcc          >  1:5014420/1‑71 (MQ=255)
                                                              tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                          >  1:4814237/1‑56 (MQ=255)
                                                              tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTa                          >  1:2221760/1‑56 (MQ=255)
                                                              tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCt             >  1:3759709/1‑69 (MQ=255)
                                                              tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc           >  1:3672754/1‑71 (MQ=255)
                                                              tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc           >  1:5528958/1‑71 (MQ=255)
                                                              tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTGc           >  1:2674080/1‑71 (MQ=255)
                                                               aaTATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt             <  1:2164482/67‑1 (MQ=255)
                                                               aaTATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTgat               >  1:3672967/1‑64 (MQ=255)
                                                               aaTATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTgat               >  1:2205858/1‑64 (MQ=255)
                                                                   ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATc                        >  1:5095556/1‑52 (MQ=255)
                                                                   ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATc                        >  1:4038816/1‑52 (MQ=255)
                                                                   ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt             <  1:1611576/63‑1 (MQ=255)
                                                                   ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt             <  1:2512567/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:2803109/1‑71 (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)
                                                                       |                                                                       
TGCTCCTTGTAAGGACTAATTATTTTTGCATATTACTTCATTAGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATG  >  minE/2571259‑2571398

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: