Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,371 T→C 6.6% intergenic (+139/‑1640) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3710TC 6.6% 351.6 / 12.9 121intergenic (+139/‑1640)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (54/59);  new base C (4/4);  total (58/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 135 total aligned reads displayed.

AGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTT  >  minE/2571301‑2571440
                                                                       |                                                                     
aGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAtttt                                                                       >  1:131548/1‑71 (MQ=255)
aGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAtttt                                                                       >  1:2605391/1‑71 (MQ=255)
 gCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAttttt                                                                      <  1:648562/71‑1 (MQ=255)
 gCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAttttt                                                                      <  1:493404/71‑1 (MQ=255)
 gCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAttttt                                                                      <  1:3307022/71‑1 (MQ=255)
       atgatTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATc                                                                  <  1:105756/69‑1 (MQ=255)
        tgatTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGt                                                               <  1:1924516/71‑1 (MQ=255)
          tttttcCTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTc                                                              <  1:682855/70‑1 (MQ=255)
           tttttCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCtt                                                            >  1:2318111/1‑71 (MQ=255)
            ttttCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                                                              <  1:580345/68‑1 (MQ=255)
            ttttCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                                           >  1:3513414/1‑71 (MQ=255)
            ttttCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                                           >  1:2079152/1‑71 (MQ=255)
             tttCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                                           <  1:2264728/70‑1 (MQ=255)
                cTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCt                                                       >  1:3596571/1‑71 (MQ=255)
                    aaTATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCg                                                   >  1:3934333/1‑71 (MQ=255)
                      tatTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAttttt                                                                      <  1:551851/50‑1 (MQ=255)
                      tatTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGAttttt                                                                      <  1:2094172/50‑1 (MQ=255)
                      tatTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTg                                                           <  1:2985924/61‑1 (MQ=255)
                        ttCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                <  1:1977310/70‑1 (MQ=255)
                        ttCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                <  1:2652559/70‑1 (MQ=255)
                        ttCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAg                                                <  1:3589264/70‑1 (MQ=255)
                         tCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGa                                              >  1:2149277/1‑71 (MQ=255)
                          cATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                             >  1:3206894/1‑71 (MQ=255)
                          cATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                             >  1:3991517/1‑71 (MQ=255)
                           aTTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGt                                                          >  1:322268/1‑57 (MQ=255)
                           aTTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                             <  1:4083275/70‑1 (MQ=255)
                            ttATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                           <  1:3833994/71‑1 (MQ=255)
                            ttATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                           <  1:338733/71‑1 (MQ=255)
                            ttATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                           <  1:3099649/71‑1 (MQ=255)
                              aTGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTc                                                              >  1:1517877/1‑50 (MQ=255)
                               tGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                                        <  1:592795/71‑1 (MQ=255)
                               tGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCt                                                       <  1:2919591/57‑1 (MQ=255)
                                 ttAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                                        <  1:4134059/69‑1 (MQ=255)
                                 ttAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      >  1:1290738/1‑71 (MQ=255)
                                 ttAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      <  1:3743354/71‑1 (MQ=255)
                                  taataaGTTC‑TGGATGAGTGTGCTCTCATTGGAttttt                                                                      >  1:3151480/1‑38 (MQ=255)
                                  taataaGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      >  1:791699/1‑70 (MQ=255)
                                  taataaGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                       >  1:2623489/1‑69 (MQ=255)
                                   aataaGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCAGTCTTGTTCTGCCGCAGGATGAGACtat                                     >  1:1415186/1‑70 (MQ=255)
                                     taaGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                  >  1:2255698/1‑71 (MQ=255)
                                       aGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                             >  1:3091884/1‑58 (MQ=255)
                                       aGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtc                                <  1:2170433/71‑1 (MQ=255)
                                        gTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                             <  1:1038152/57‑1 (MQ=255)
                                          tC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             >  1:1486724/1‑71 (MQ=255)
                                           tatGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgt                              >  1:2245535/3‑70 (MQ=255)
                                           c‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATgaga                                         >  1:4049697/1‑58 (MQ=255)
                                           c‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATgaga                                         >  1:314431/1‑58 (MQ=255)
                                           c‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             >  1:3026843/1‑70 (MQ=255)
                                           c‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             <  1:1431206/70‑1 (MQ=255)
                                           c‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             >  1:3342624/1‑70 (MQ=255)
                                             tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             <  1:1658026/69‑1 (MQ=255)
                                               gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                          <  1:3347870/70‑1 (MQ=255)
                                               gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                           <  1:860484/69‑1 (MQ=255)
                                                aTGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             >  1:1544499/1‑66 (MQ=255)
                                                aTGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGAt                        <  1:95336/71‑1 (MQ=255)
                                                 tGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATc                       <  1:1067535/71‑1 (MQ=255)
                                                   aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGc                     <  1:1993704/71‑1 (MQ=255)
                                                   aGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGc                     <  1:3025761/71‑1 (MQ=255)
                                                    gtgtGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATc                       >  1:3330941/1‑68 (MQ=255)
                                                      gtgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                  >  1:726925/1‑71 (MQ=255)
                                                       tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                  >  1:1159789/1‑70 (MQ=255)
                                                       tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                  <  1:918288/70‑1 (MQ=255)
                                                       tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                  <  1:3498040/70‑1 (MQ=255)
                                                         ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGa                                              >  1:3274310/1‑40 (MQ=255)
                                                         ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat               <  1:1811566/71‑1 (MQ=255)
                                                         ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat               <  1:353534/71‑1 (MQ=255)
                                                         ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                >  1:3285930/1‑70 (MQ=255)
                                                           ctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                  >  1:468776/1‑50 (MQ=255)
                                                           ctcATTGGATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga             >  1:375613/1‑71 (MQ=255)
                                                           ctcATTGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACtaatga             >  1:3027179/1‑71 (MQ=255)
                                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGa                                              >  1:1817655/1‑37 (MQ=255)
                                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                           <  1:3589774/40‑1 (MQ=255)
                                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      <  1:2521207/45‑1 (MQ=255)
                                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      <  1:224914/45‑1 (MQ=255)
                                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga             >  1:2338891/1‑70 (MQ=255)
                                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga             >  1:3018885/1‑70 (MQ=255)
                                                             cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAGTATACATGGTCGAGATCGCAatgat               <  1:3102297/67‑1 (MQ=255)
                                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                    >  1:3282848/1‑45 (MQ=255)
                                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga             <  1:3253065/68‑1 (MQ=255)
                                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga             <  1:2129203/68‑1 (MQ=255)
                                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat          <  1:4065114/71‑1 (MQ=255)
                                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat          <  1:2052191/71‑1 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                           <  1:3326977/37‑1 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                           <  1:2805248/37‑1 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      <  1:715103/42‑1 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                             <  1:1932209/51‑1 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat          >  1:3693318/1‑70 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg         <  1:3908495/71‑1 (MQ=255)
                                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATg         <  1:1627915/71‑1 (MQ=255)
                                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtcg                            >  1:2749450/1‑51 (MQ=255)
                                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                           >  1:3851271/1‑52 (MQ=255)
                                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat          >  1:4111812/1‑69 (MQ=255)
                                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata           >  1:1592109/1‑68 (MQ=255)
                                                                 ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACta                                      <  1:1387718/40‑1 (MQ=255)
                                                                 ggATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg       >  1:1950220/1‑71 (MQ=255)
                                                                 ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa        <  1:1875957/70‑1 (MQ=255)
                                                                    ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttg    <  1:2774117/71‑1 (MQ=255)
                                                                    ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttg    <  1:416103/71‑1 (MQ=255)
                                                                     tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttg    <  1:2425515/70‑1 (MQ=255)
                                                                       tttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt  <  1:455843/70‑1 (MQ=255)
                                                                       |                                                                     
AGCTTCGATGATTTTTCTTTAATATTCATTATGTTAATAAGTTC‑TGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTT  >  minE/2571301‑2571440

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: ‑