Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,394 G→T 17.9% intergenic (+162/‑1617) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3940GT17.9% 397.8 / 50.5 150intergenic (+162/‑1617)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (54/69);  new base T (15/12);  total (69/81)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-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)
   aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTATATCCGTCTTGTTCTGCCGCAGGAt                                                                   <  1:5026670/70‑1 (MQ=255)
     tATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:5231582/1‑68 (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:2060905/1‑70 (MQ=255)
              aaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 >  1:246516/1‑61 (MQ=255)
                   cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:3511341/70‑1 (MQ=255)
                   cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:5097885/70‑1 (MQ=255)
                    tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:774512/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:1469927/1‑51 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:2566562/1‑51 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                                   >  1:1823475/1‑51 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                   <  1:5436525/67‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                 <  1:666272/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                 <  1:4143204/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                 <  1:4358569/69‑1 (MQ=255)
                      gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGACGCAGGATGAGACTATACATCGTCgag                                                <  1:3829152/70‑1 (MQ=255)
                        tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:3844992/46‑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: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:3201858/40‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:886221/40‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:5117893/40‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatg                                       <  1:4960781/71‑1 (MQ=255)
                              tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                        <  1:4984534/70‑1 (MQ=255)
                               gCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                <  1:4165716/61‑1 (MQ=255)
                                ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                     <  1:151787/71‑1 (MQ=255)
                                atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                            <  1:4986953/63‑1 (MQ=255)
                                atctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACg                                            <  1:4993226/63‑1 (MQ=255)
                                 tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:814179/1‑49 (MQ=255)
                                  ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                      <  1:1266480/36‑1 (MQ=255)
                                  ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                                   >  1:1943869/1‑71 (MQ=255)
                                   tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgaga                                               >  1:4175455/1‑58 (MQ=255)
                                   tcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCacta                                       <  1:4315466/66‑4 (MQ=255)
                                   tcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCacta                                       <  1:1135986/66‑4 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:1435367/1‑46 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                          >  1:1706652/1‑46 (MQ=255)
                                    cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                    <  1:3904393/68‑1 (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)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                                 <  1:4734283/70‑3 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:3556466/1‑71 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:2753796/1‑71 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:1715371/1‑71 (MQ=255)
                                     aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                >  1:1508719/1‑71 (MQ=255)
                                      ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                                 <  1:3939917/37‑1 (MQ=255)
                                      ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                          <  1:1107040/60‑1 (MQ=255)
                                      ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                        <  1:5173566/62‑1 (MQ=255)
                                       tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                >  1:5176006/1‑69 (MQ=255)
                                        ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                                     <  1:5399116/63‑1 (MQ=255)
                                        ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                      >  1:4112445/1‑62 (MQ=255)
                                        ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                             >  1:2577682/1‑71 (MQ=255)
                                        ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                             >  1:4664338/1‑71 (MQ=255)
                                        ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                              <  1:1182047/70‑1 (MQ=255)
                                        ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                              <  1:947627/70‑1 (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:4599096/1‑55 (MQ=255)
                                           tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                          >  1:1970997/1‑55 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtcg                                                  <  1:1363715/46‑1 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGc                                           <  1:4813418/53‑1 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                    >  1:2172702/1‑60 (MQ=255)
                                            tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                                    >  1:3515198/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)
                                                      tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                 <  1:5372862/69‑1 (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:4301107/46‑1 (MQ=255)
                                                           ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   >  1:3234548/1‑46 (MQ=255)
                                                           ttCTGCCCCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                   >  1:4292406/1‑46 (MQ=255)
                                                             cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc        >  1:4362226/1‑71 (MQ=255)
                                                             cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc        >  1:5671987/1‑71 (MQ=255)
                                                             cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct         >  1:3810355/1‑70 (MQ=255)
                                                              tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                        <  1:3935302/38‑1 (MQ=255)
                                                              tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                        <  1:64923/38‑1 (MQ=255)
                                                              tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACgg           >  1:5083692/1‑67 (MQ=255)
                                                                ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                      >  1:4033043/1‑38 (MQ=255)
                                                                 cGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCg               >  1:3029933/1‑60 (MQ=255)
                                                                    aGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc             <  1:3721742/59‑1 (MQ=255)
                                                                     |                                                                    
TTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCA  >  minE/2571325‑2571462

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: