Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,392 C→G 16.8% intergenic (+160/‑1619) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3920CG16.8% 375.0 / 60.6 161intergenic (+160/‑1619)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (55/79);  new base G (15/12);  total (70/91)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 162 total aligned reads displayed.

TTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCA  >  minE/2571325‑2571462
                                                                   |                                                                      
ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAgt                                                                     <  1:5305193/71‑2 (MQ=255)
ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:1864496/70‑1 (MQ=255)
ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:3328060/70‑1 (MQ=255)
  cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:2803109/1‑71 (MQ=255)
  cATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:1549294/1‑71 (MQ=255)
   aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   <  1:3930766/70‑1 (MQ=255)
   aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATg                                                                  <  1:4985725/71‑1 (MQ=255)
       tGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                                                              <  1:4784700/71‑1 (MQ=255)
         ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                                            <  1:1616543/71‑1 (MQ=255)
         ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                                             <  1:3411049/70‑1 (MQ=255)
         ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                                             <  1:596909/70‑1 (MQ=255)
          taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                                            >  1:4223154/1‑70 (MQ=255)
            ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          <  1:4326404/70‑1 (MQ=255)
            ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:2060905/1‑70 (MQ=255)
              aaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 >  1:246516/1‑61 (MQ=255)
                   cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:5004618/70‑1 (MQ=255)
                   cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:3511341/70‑1 (MQ=255)
                   cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:5097885/70‑1 (MQ=255)
                    tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:3932380/69‑1 (MQ=255)
                    tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:774512/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:1823475/1‑51 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:1469927/1‑51 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:2566562/1‑51 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:5436525/67‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                 <  1:4149641/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                 <  1:4143204/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                 <  1:666272/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGACGCAGGATGAGACTATACATCGTCgag                                                <  1:3829152/70‑1 (MQ=255)
                        tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                        <  1:3179875/60‑1 (MQ=255)
                        tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                        <  1:1392556/60‑1 (MQ=255)
                        tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                        <  1:4685267/60‑1 (MQ=255)
                         gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                            >  1:463623/1‑71 (MQ=255)
                         gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                            >  1:3667859/1‑71 (MQ=255)
                         gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                            >  1:1537116/1‑71 (MQ=255)
                          aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 <  1:1982580/49‑1 (MQ=255)
                           gtgtGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                          <  1:227182/71‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:886221/40‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:3201858/40‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                       <  1:3810209/71‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                        <  1:4984534/70‑1 (MQ=255)
                               gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                <  1:417316/61‑1 (MQ=255)
                                ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                     <  1:151787/71‑1 (MQ=255)
                                atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                            <  1:4993226/63‑1 (MQ=255)
                                 tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:814179/1‑49 (MQ=255)
                                  ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:398488/36‑1 (MQ=255)
                                  ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:1266480/36‑1 (MQ=255)
                                  ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 >  1:3864690/1‑41 (MQ=255)
                                  ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                   >  1:1943869/1‑71 (MQ=255)
                                   tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   >  1:4314960/1‑70 (MQ=255)
                                   tcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCacta                                       <  1:1135986/66‑4 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:1706652/1‑46 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:510010/1‑46 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:1435367/1‑46 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                 >  1:1951030/1‑71 (MQ=255)
                                     aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   <  1:552476/68‑1 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                              >  1:1906316/1‑56 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                              >  1:1372611/1‑56 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCActaatg                                    <  1:4527043/67‑1 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:3556466/1‑71 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:1715371/1‑71 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:1508719/1‑71 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:2753796/1‑71 (MQ=255)
                                      ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 <  1:4703125/37‑1 (MQ=255)
                                      ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 <  1:3939917/37‑1 (MQ=255)
                                      ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                          <  1:1107040/60‑1 (MQ=255)
                                       tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                >  1:5176006/1‑69 (MQ=255)
                                        ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                     <  1:5399116/63‑1 (MQ=255)
                                        ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                             >  1:2577682/1‑71 (MQ=255)
                                        ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTTTGCATTGTTGAGAACGCACTAATGACATGa                              <  1:5040939/70‑1 (MQ=255)
                                        ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                              <  1:1182047/70‑1 (MQ=255)
                                        ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                              <  1:947627/70‑1 (MQ=255)
                                         gATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgt                            >  1:383292/1‑71 (MQ=255)
                                         gATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGaat                            >  1:3122013/1‑69 (MQ=255)
                                           ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                               <  1:2994731/66‑1 (MQ=255)
                                           ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                               <  1:231130/66‑1 (MQ=255)
                                           tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                          >  1:1970997/1‑55 (MQ=255)
                                           tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                          >  1:5134465/1‑55 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtcg                                                  <  1:1363715/46‑1 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                    >  1:3515198/1‑60 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                    >  1:2172702/1‑60 (MQ=255)
                                              tttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtc                       >  1:2823603/1‑71 (MQ=255)
                                                tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                      <  1:2406726/70‑1 (MQ=255)
                                                  tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   <  1:1538276/55‑1 (MQ=255)
                                                  tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcctcc                   >  1:2529112/1‑71 (MQ=255)
                                                    cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct                  <  1:5342248/70‑1 (MQ=255)
                                                      tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                 <  1:5372862/69‑1 (MQ=255)
                                                        ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                >  1:4034113/1‑52 (MQ=255)
                                                        ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg            >  1:1625714/1‑71 (MQ=255)
                                                        ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg            >  1:1619006/1‑71 (MQ=255)
                                                          gTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc             <  1:1230706/69‑1 (MQ=255)
                                                           ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   >  1:3234548/1‑46 (MQ=255)
                                                           ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   <  1:4301107/46‑1 (MQ=255)
                                                             cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc        >  1:423132/1‑71 (MQ=255)
                                                             cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc        >  1:5671987/1‑71 (MQ=255)
                                                             cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc        >  1:4362226/1‑71 (MQ=255)
                                                              tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                        <  1:64923/38‑1 (MQ=255)
                                                                ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct       <  1:5208430/69‑1 (MQ=255)
                                                                ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct       <  1:482278/69‑1 (MQ=255)
                                                                 cGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCg               >  1:3029933/1‑60 (MQ=255)
                                                                   cAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa  <  1:4786608/71‑1 (MQ=255)
                                                                   |                                                                      
TTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCA  >  minE/2571325‑2571462

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