Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,391 G→C 17.4% intergenic (+159/‑1620) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3910GC17.4% 373.0 / 62.9 161intergenic (+159/‑1620)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (54/79);  new base C (16/12);  total (70/91)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

ATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCT  >  minE/2571322‑2571457
                                                                     |                                                                  
atatTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCg                                                                    <  1:1877765/70‑1 (MQ=255)
   ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAgt                                                                <  1:5305193/71‑2 (MQ=255)
   ttCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                 <  1:1864496/70‑1 (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:596909/70‑1 (MQ=255)
             taataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                                       >  1:4223154/1‑70 (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:774512/69‑1 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                              >  1:2566562/1‑51 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                              >  1:1469927/1‑51 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                              >  1:1823475/1‑51 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                              <  1:5436525/67‑1 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                           <  1:3413843/70‑1 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                            <  1:4358569/69‑1 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                            <  1:666272/69‑1 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                            <  1:2161337/69‑1 (MQ=255)
                         gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                            <  1:4143204/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)
                                 tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                 <  1:886221/40‑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:4993226/63‑1 (MQ=255)
                                    tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                     >  1:814179/1‑49 (MQ=255)
                                    tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                     >  1:2847201/1‑49 (MQ=255)
                                     ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                                 <  1:1266480/36‑1 (MQ=255)
                                     ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                            >  1:2578423/1‑41 (MQ=255)
                                     ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                              >  1:1943869/1‑71 (MQ=255)
                                     ctcATTAGATTTCTCATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga                              >  1:2317238/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)
                                       cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                     >  1:510010/1‑46 (MQ=255)
                                       cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                               <  1:3904393/68‑1 (MQ=255)
                                       cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                               <  1:1937277/68‑1 (MQ=255)
                                       cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                            >  1:2811266/1‑71 (MQ=255)
                                        aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                              <  1:552476/68‑1 (MQ=255)
                                        aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                              <  1:3581123/68‑1 (MQ=255)
                                        aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                           <  1:3135749/71‑1 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                         >  1:1372611/1‑56 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                            <  1:2542554/70‑3 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                            <  1:4734283/70‑3 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatgaca                            <  1:3093088/70‑3 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                           >  1:2753796/1‑71 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                           >  1:1508719/1‑71 (MQ=255)
                                        aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                           >  1:1715371/1‑71 (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)
                                           ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                 >  1:4112445/1‑62 (MQ=255)
                                           ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                        >  1:4664338/1‑71 (MQ=255)
                                           ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTTTGCATTGTTGAGAACGCACTAATGACATGa                         <  1:5040939/70‑1 (MQ=255)
                                           ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                         <  1:947627/70‑1 (MQ=255)
                                           ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                         <  1:1182047/70‑1 (MQ=255)
                                              ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg                          <  1:2994731/66‑1 (MQ=255)
                                              tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                     >  1:4599096/1‑55 (MQ=255)
                                              tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                     >  1:5134465/1‑55 (MQ=255)
                                               tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtcg                                             <  1:1363715/46‑1 (MQ=255)
                                               tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                               >  1:2172702/1‑60 (MQ=255)
                                                   tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc                 <  1:2406726/70‑1 (MQ=255)
                                                     tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                              <  1:1538276/55‑1 (MQ=255)
                                                     tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                              <  1:3553576/55‑1 (MQ=255)
                                                       cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct             <  1:5342248/70‑1 (MQ=255)
                                                       cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg            >  1:3330444/1‑71 (MQ=255)
                                                         tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg            <  1:3229031/69‑1 (MQ=255)
                                                         tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg            <  1:5372862/69‑1 (MQ=255)
                                                           ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc        >  1:1955615/1‑71 (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)
                                                              ttCTGCCCCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                              >  1:4292406/1‑46 (MQ=255)
                                                               tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                    >  1:2466197/1‑39 (MQ=255)
                                                                cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                              <  1:1971366/44‑1 (MQ=255)
                                                                cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc   >  1:5671987/1‑71 (MQ=255)
                                                                cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc   >  1:4362226/1‑71 (MQ=255)
                                                                cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct    >  1:3810355/1‑70 (MQ=255)
                                                                 tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                   <  1:64923/38‑1 (MQ=255)
                                                                 tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                   <  1:3935302/38‑1 (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)
                                                                     |                                                                  
ATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCT  >  minE/2571322‑2571457

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: