Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,410 C→T 19.4% intergenic (+178/‑1601) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,4100CT19.4% 292.4 / 58.1 139intergenic (+178/‑1601)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (54/58);  new base T (15/12);  total (69/70)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 151 total aligned reads displayed.

CTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGC  >  minE/2571344‑2571479
                                                                  |                                                                     
cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                                    <  1:5004618/70‑1 (MQ=255)
cTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                                    <  1:3511341/70‑1 (MQ=255)
 tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                                    <  1:774512/69‑1 (MQ=255)
 tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                                    <  1:3932380/69‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                                                    <  1:5436525/67‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                                                 <  1:3413843/70‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                                  <  1:666272/69‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                                  <  1:4149641/69‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                                  <  1:4358569/69‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                                                  <  1:2161337/69‑1 (MQ=255)
   gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGACGCAGGATGAGACTATACATCGTCgag                                                                 <  1:3829152/70‑1 (MQ=255)
      gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                             >  1:3667859/1‑71 (MQ=255)
      gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                             >  1:1537116/1‑71 (MQ=255)
       gctgtgATCTCATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                               >  1:1551697/3‑67 (MQ=255)
       aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                                             <  1:1820505/70‑1 (MQ=255)
        gtgtGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCa                                                           <  1:227182/71‑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:4993226/63‑1 (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)
                 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)
                  aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                    <  1:552476/68‑1 (MQ=255)
                  aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                                                 <  1:3135749/71‑1 (MQ=255)
                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                               >  1:1906316/1‑56 (MQ=255)
                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTgagaa                                                               >  1:1460628/1‑56 (MQ=255)
                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                 >  1:3556466/1‑71 (MQ=255)
                  aTTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                                                 >  1:1715371/1‑71 (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)
                     ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                                              >  1:2577682/1‑71 (MQ=255)
                     ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                                              >  1:4664338/1‑71 (MQ=255)
                     ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                                               <  1:947627/70‑1 (MQ=255)
                      gATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGaat                                             >  1:3122013/1‑69 (MQ=255)
                        tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                           >  1:4599096/1‑55 (MQ=255)
                        tttcttATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCa                                                           >  1:5134465/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)
                          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)
                               tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcctcc                                    >  1:2529112/1‑71 (MQ=255)
                                 cGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTcctcct                                   <  1:5342248/70‑1 (MQ=255)
                                   tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTg                                  <  1:5372862/69‑1 (MQ=255)
                                     ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                              >  1:1955615/1‑71 (MQ=255)
                                     ttGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                              >  1:1260997/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)
                                         tCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                                          >  1:2466197/1‑39 (MQ=255)
                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                                                    <  1:1971366/44‑1 (MQ=255)
                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                         >  1:423132/1‑71 (MQ=255)
                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctc                         >  1:5671987/1‑71 (MQ=255)
                                          cTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGct                          >  1:3810355/1‑70 (MQ=255)
                                           tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                                                         <  1:64923/38‑1 (MQ=255)
                                           tGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACgg                            >  1:5083692/1‑67 (MQ=255)
                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                       >  1:4033043/1‑38 (MQ=255)
                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                       >  1:1664118/1‑38 (MQ=255)
                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAatga                                                       >  1:2298830/1‑38 (MQ=255)
                                             ccGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGctct                        <  1:5208430/69‑1 (MQ=255)
                                              cGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCg                                >  1:3029933/1‑60 (MQ=255)
                                                cAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                   <  1:4786608/71‑1 (MQ=255)
                                                 aGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGAc                              <  1:3721742/59‑1 (MQ=255)
                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgtt                                            <  1:4446608/43‑1 (MQ=255)
                                                   gATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaa                <  1:5505732/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)
                                                     tGAGACGATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTaa                      >  1:1459/1‑63 (MQ=255)
                                                       agaCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAag             <  1:3733606/70‑1 (MQ=255)
                                                         aCTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCaaaaaa              >  1:4354855/1‑67 (MQ=255)
                                                           tataCATCGTCGAGATCGCAATGATGTTATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg        >  1:3801726/1‑71 (MQ=255)
                                                             taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCTACGGCTCTAATCAAAAAAGAGCAgg       <  1:3670887/70‑1 (MQ=255)
                                                             taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAAt                     <  1:5045480/56‑1 (MQ=255)
                                                             taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAgg       <  1:2745143/70‑1 (MQ=255)
                                                             taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg        >  1:2680435/1‑69 (MQ=255)
                                                             taCATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg        >  1:4592541/1‑69 (MQ=255)
                                                               cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                   >  1:4450503/1‑56 (MQ=255)
                                                               cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCa                   >  1:3092649/1‑56 (MQ=255)
                                                               cATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGt    >  1:4219500/1‑71 (MQ=255)
                                                                aTCGTCGAGATCGCAATGATGATATGAGTTGTtcct                                      >  1:4080725/1‑36 (MQ=255)
                                                                aTCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAg        <  1:1990394/66‑1 (MQ=255)
                                                                 tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc  >  1:4684665/1‑71 (MQ=255)
                                                                 tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc  >  1:4018276/1‑71 (MQ=255)
                                                                 tcgtcgAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGc  >  1:1351916/1‑71 (MQ=255)
                                                                  |                                                                     
CTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGGCTCTAATCAAAAAAGAGCAGGAGTGC  >  minE/2571344‑2571479

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: