Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,398 G→C 18.2% intergenic (+166/‑1613) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3980GC18.2% 335.1 / 60.8 148intergenic (+166/‑1613)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (50/71);  new base C (15/12);  total (65/83)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 155 total aligned reads displayed.

ATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAA  >  minE/2571328‑2571465
                                                                      |                                                                   
aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATg                                                                     <  1:4985725/71‑1 (MQ=255)
    tGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                                                                 <  1:4784700/71‑1 (MQ=255)
       taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                                               >  1:4223154/1‑70 (MQ=255)
         ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                             <  1:4326404/70‑1 (MQ=255)
                cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                      <  1:5097885/70‑1 (MQ=255)
                cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                      <  1:5004618/70‑1 (MQ=255)
                cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                      <  1:3511341/70‑1 (MQ=255)
                 tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                      <  1:3932380/69‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                      <  1:5436525/67‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgaga                                                  <  1:5619635/71‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                   <  1:3413843/70‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                    <  1:4143204/69‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                    <  1:2161337/69‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                    <  1:4358569/69‑1 (MQ=255)
                   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGACGCAGGATGAGACTATACATCGTCgag                                                   <  1:3829152/70‑1 (MQ=255)
                     tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                                           <  1:3179875/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)
                       gctgtgATCTCATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                 >  1:1551697/3‑67 (MQ=255)
                       aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                    <  1:1982580/49‑1 (MQ=255)
                       aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                               <  1:1820505/70‑1 (MQ=255)
                        gtgtGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                             <  1:227182/71‑1 (MQ=255)
                           tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                          <  1:3810209/71‑1 (MQ=255)
                           tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                          <  1:4960781/71‑1 (MQ=255)
                            gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                   <  1:1624458/61‑1 (MQ=255)
                            gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                   <  1:4165716/61‑1 (MQ=255)
                             ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                        <  1:151787/71‑1 (MQ=255)
                             atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                               <  1:1855778/63‑1 (MQ=255)
                             atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                               <  1:4986953/63‑1 (MQ=255)
                              tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                             >  1:814179/1‑49 (MQ=255)
                               ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                    >  1:2578423/1‑41 (MQ=255)
                               ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                      >  1:2317238/1‑71 (MQ=255)
                                tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgaga                                                  >  1:4175455/1‑58 (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:1435367/1‑46 (MQ=255)
                                 cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                       <  1:1937277/68‑1 (MQ=255)
                                 cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                       <  1:3904393/68‑1 (MQ=255)
                                 cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                    >  1:1951030/1‑71 (MQ=255)
                                 cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                                    >  1:2811266/1‑71 (MQ=255)
                                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                 >  1:1372611/1‑56 (MQ=255)
                                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCActaatg                                       <  1:4527043/67‑1 (MQ=255)
                                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                    <  1:2542554/70‑3 (MQ=255)
                                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                   >  1:3556466/1‑71 (MQ=255)
                                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                   >  1:1508719/1‑71 (MQ=255)
                                   ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                    <  1:3939917/37‑1 (MQ=255)
                                   ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                    <  1:4703125/37‑1 (MQ=255)
                                   ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                             <  1:1107040/60‑1 (MQ=255)
                                   ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                           <  1:5173566/62‑1 (MQ=255)
                                     ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                        <  1:5399116/63‑1 (MQ=255)
                                     ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                         >  1:4112445/1‑62 (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)
                                        tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                             >  1:5134465/1‑55 (MQ=255)
                                         tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtcg                                                     <  1:1363715/46‑1 (MQ=255)
                                          ttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt                           >  1:622813/1‑71 (MQ=255)
                                           tttATCGGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                            >  1:731348/1‑53 (MQ=255)
                                           tttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtc                          >  1:2823603/1‑71 (MQ=255)
                                             tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                         <  1:2406726/70‑1 (MQ=255)
                                               tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                      <  1:3553576/55‑1 (MQ=255)
                                               tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                      <  1:1538276/55‑1 (MQ=255)
                                                 cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct                     <  1:5342248/70‑1 (MQ=255)
                                                 cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                    >  1:3330444/1‑71 (MQ=255)
                                                   tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                    <  1:5372862/69‑1 (MQ=255)
                                                   tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                    <  1:3229031/69‑1 (MQ=255)
                                                     ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                >  1:1955615/1‑71 (MQ=255)
                                                     ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                >  1:1260997/1‑71 (MQ=255)
                                                     ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGAATGATGATATGAGTTGTTCCTCCTGCGACg               >  1:1619006/1‑71 (MQ=255)
                                                       gTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                <  1:1230706/69‑1 (MQ=255)
                                                        ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                      <  1:4301107/46‑1 (MQ=255)
                                                         tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                            >  1:2466197/1‑39 (MQ=255)
                                                         tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                             >  1:813073/1‑38 (MQ=255)
                                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                      <  1:1971366/44‑1 (MQ=255)
                                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc           >  1:5671987/1‑71 (MQ=255)
                                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc           >  1:423132/1‑71 (MQ=255)
                                                           tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                           <  1:64923/38‑1 (MQ=255)
                                                           tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACgg              >  1:5083692/1‑67 (MQ=255)
                                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                         >  1:1664118/1‑38 (MQ=255)
                                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                         >  1:4033043/1‑38 (MQ=255)
                                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                         >  1:2298830/1‑38 (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)
                                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                              <  1:4446608/43‑1 (MQ=255)
                                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                              <  1:2142377/43‑1 (MQ=255)
                                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                              <  1:3708620/43‑1 (MQ=255)
                                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa  <  1:5505732/71‑1 (MQ=255)
                                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa  <  1:258067/71‑1 (MQ=255)
                                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa  <  1:3290580/71‑1 (MQ=255)
                                                                     tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa        >  1:5210211/1‑63 (MQ=255)
                                                                     tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa  >  1:2546351/1‑69 (MQ=255)
                                                                     tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa     >  1:4208551/1‑66 (MQ=255)
                                                                     tGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa     >  1:711638/1‑66 (MQ=255)
                                                                     tGAGACGATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa        >  1:1459/1‑63 (MQ=255)
                                                                      |                                                                   
ATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAA  >  minE/2571328‑2571465

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