Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,872,450 +G coding (338/387 nt) yidW ← regulator protein for dgo operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,872,4501.G99.3% 365.2 / ‑3.8 140coding (338/387 nt)yidWregulator protein for dgo operon
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (72/67);  total (72/68)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 146 total aligned reads displayed.

GTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGT  >  NC_000913/3872375‑3872519
                                                                             |                                                                     
gTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCt                                                                        <  1:2824116/76‑1 (MQ=255)
 tCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtg                                                                       <  1:3143339/76‑2 (MQ=255)
  ccTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCtgc                                                                      <  1:2135880/76‑3 (MQ=255)
  ccTGATTGTTGGAAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAACTCGCGCGCTCCGCtgc                                                                      >  1:240706/1‑74 (MQ=255)
      aTTGTTGGCAATCATCTCGTTCGGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCaa                                                                  >  1:122269/1‑76 (MQ=255)
        tgttgtCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACg                                                                <  2:3299075/76‑1 (MQ=255)
         gttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCCACgg                                                               >  2:1479642/1‑76 (MQ=255)
         gttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACgg                                                               <  2:1301232/76‑1 (MQ=255)
         gttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACgg                                                               >  2:4679487/1‑76 (MQ=255)
         gttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACgg                                                               >  1:2369217/1‑76 (MQ=255)
          ttgGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACggg                                                              <  1:1210332/76‑1 (MQ=255)
            ggCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCa                                                            >  1:987429/1‑76 (MQ=255)
              cAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAAt                                                          <  1:422430/76‑1 (MQ=255)
               aaTCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATc                                                         <  2:2115654/76‑1 (MQ=255)
               aaTCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATc                                                         <  1:2085416/76‑1 (MQ=255)
               aaTCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTACTGCCCAACGGGCAATc                                                         <  2:446239/76‑1 (MQ=255)
                atcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCCGCCCAACGGGCAATCg                                                        >  1:3229656/1‑76 (MQ=255)
                 tcatcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGc                                                       <  2:3575142/76‑1 (MQ=255)
                  catcTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGcc                                                      <  2:3683219/76‑1 (MQ=255)
                        gTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCCACGGGCAATCCCCCGTTcc                                                >  2:3101687/1‑76 (MQ=255)
                         ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCACTCGCCGGTTccc                                               >  2:4868381/1‑75 (MQ=255)
                         ttCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTcca                                               >  1:1348178/1‑76 (MQ=255)
                          tCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccac                                              <  1:2569595/76‑1 (MQ=255)
                          tCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccac                                              >  1:4701705/1‑76 (MQ=255)
                           cAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacc                                             >  1:3137328/1‑76 (MQ=255)
                            aGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacca                                            >  1:4156599/1‑76 (MQ=255)
                            aGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTccacca                                            >  1:3366102/1‑76 (MQ=255)
                             gcgcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAg                                           >  2:2193324/1‑76 (MQ=255)
                               gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAt                                         <  1:642293/76‑1 (MQ=255)
                               gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAt                                         >  2:1220849/1‑76 (MQ=255)
                               gcCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAt                                         <  1:2631047/76‑1 (MQ=255)
                                ccGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAtt                                        <  1:78062/76‑1 (MQ=255)
                                 cGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                       <  2:3228004/76‑1 (MQ=255)
                                  gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTc                                      >  2:823085/1‑76 (MQ=255)
                                  gATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGCTCCACCAGATTTc                                      >  2:155521/1‑76 (MQ=255)
                                    ttCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGc                                    <  1:395285/76‑1 (MQ=255)
                                    ttCAAGCTGCGCCAGATCGCTGGAAGTCGACCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGc                                    >  2:3317209/1‑76 (MQ=255)
                                      cAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCAc                                  >  2:2996132/1‑76 (MQ=255)
                                       aaccTGCGCCAGATCGCTGGAAGCCGCGCGCTCTGCTGCCCAACCGCCACTCGCCCTCTCCACCCCATCTCGCACt                                 >  2:35680/4‑76 (MQ=255)
                                        aTCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACtt                                <  1:3063305/76‑1 (MQ=255)
                                         tCTGCGCCAGCTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                               <  1:4744317/76‑1 (MQ=255)
                                         tCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                               <  1:1384238/76‑1 (MQ=255)
                                         tCTGCGCCAGATCGCTGGAAGCCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTc                               >  2:495518/1‑76 (MQ=255)
                                          cTGCGCCAGATCGCTGGAAGTCTCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                              <  1:3431720/76‑1 (MQ=255)
                                          cTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                              <  1:2940339/76‑1 (MQ=255)
                                          cTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCg                              >  1:3039528/1‑76 (MQ=255)
                                           tGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGc                             <  1:4528553/76‑1 (MQ=255)
                                            gcgcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCt                            >  2:4300156/1‑76 (MQ=255)
                                             cgcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTc                           >  1:4124139/1‑76 (MQ=255)
                                              gcCAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCa                          <  1:4509002/76‑1 (MQ=255)
                                               ccAGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCAt                         <  2:1873603/76‑1 (MQ=255)
                                                 aGATCGGGGGAAGGCGGGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATgg                       <  2:2362520/76‑1 (MQ=255)
                                                 aGATCGCTGGAGTGCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATgg                       <  1:2409032/76‑1 (MQ=255)
                                                 aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                       <  2:2292077/60‑1 (MQ=255)
                                                 aGATCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGAttt                                       >  1:2292078/1‑60 (MQ=255)
                                                   aTCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCa                     >  2:3737842/1‑76 (MQ=255)
                                                    tCGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                    <  1:834300/76‑1 (MQ=255)
                                                    tCGCTGGAAGTCGCGCGCCCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCAc                    <  1:247392/76‑1 (MQ=255)
                                                     ctcTGTAAGGCGCGCGCTCGGCCGCGCAACGGGCAATCGCCGGTTCCACCAGATTGCGCACTTCGCTCATGGCACt                   <  2:4477425/74‑1 (MQ=255)
                                                     cGCTGTAAGTCGCTCGCTCTGCTGCCCAACGGGCAGTCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                   <  1:706245/76‑1 (MQ=255)
                                                     cGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                   <  1:4148298/76‑1 (MQ=255)
                                                     cGCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACt                   <  2:4283096/76‑1 (MQ=255)
                                                      gCTGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTg                  >  1:2926424/1‑76 (MQ=255)
                                                        tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCCCTTCGCTCATGGCACTgct                >  2:1570337/1‑74 (MQ=255)
                                                        tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCACGGCACTGAt                >  2:3810292/1‑76 (MQ=255)
                                                        tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGGTCCATCAGATTTCGCACTTCGCTCATGGCACTGAt                >  1:3133888/1‑76 (MQ=255)
                                                        tGGAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCAGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGAt                >  1:853024/1‑76 (MQ=255)
                                                        tGGAAGTCGCGCGCGCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGAt                <  2:4375409/76‑1 (MQ=255)
                                                        gggAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGAt                <  2:92561/75‑1 (MQ=255)
                                                         ggAAGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATa               >  1:2824679/1‑76 (MQ=255)
                                                         ggAAGTCGCGCGCTCTGCTGCCCAACGGGACATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATa               >  1:3897710/1‑76 (MQ=255)
                                                          gAAGTCGCGAGCTCTGCTGCACAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATaa              >  2:611658/1‑76 (MQ=255)
                                                           aaGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCTCCGATAAg             >  2:55241/1‑76 (MQ=255)
                                                           aaGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATaac             >  2:4836926/1‑75 (MQ=255)
                                                           aaGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAg             <  1:1692437/76‑1 (MQ=255)
                                                           aaGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAg             >  2:1673209/1‑76 (MQ=255)
                                                           aaGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAg             <  2:2695036/76‑1 (MQ=255)
                                                            aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc            <  2:1262382/76‑1 (MQ=255)
                                                            aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc            <  1:2619394/76‑1 (MQ=255)
                                                            aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc            <  2:3545203/76‑1 (MQ=255)
                                                            aGTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGc            >  1:3261363/1‑76 (MQ=255)
                                                             gTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGcc           <  2:4156599/76‑1 (MQ=255)
                                                             gTCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGcc           <  2:2615210/76‑1 (MQ=255)
                                                              tCGCGCGCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCg          >  2:4255360/1‑76 (MQ=255)
                                                              gcgcgggCGCTGCTGCCCAACGGGCAATCGCCGGTTCCGCCAGATTGCGCACTTCGCTCCTGGCACTGATAAACCg          <  2:4530359/75‑1 (MQ=255)
                                                               cgcgcgCTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt         >  2:3449553/1‑76 (MQ=255)
                                                               cgcgcgCTCTGCTGCCCAACGGGCAATCGCCGGCCCCACCAGATTTCGCACTTCGCTCATGGCCCTGATAAGCCGt         >  2:112962/1‑76 (MQ=255)
                                                               cgcgcgCTCTGCTGCCCAACGGGCAATCGCCGGATCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGt         <  1:60573/76‑1 (MQ=255)
                                                                gcgcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTg        <  2:3238035/76‑1 (MQ=255)
                                                                gcgcgcTCTGCTGCCCAACGGGAAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTg        <  2:861358/76‑1 (MQ=255)
                                                                  gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg      >  1:1954124/1‑76 (MQ=255)
                                                                  gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg      >  1:4594712/1‑76 (MQ=255)
                                                                  gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg      >  1:2669332/1‑76 (MQ=255)
                                                                  gcgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTggg      >  2:3454112/1‑76 (MQ=255)
                                                                   cgcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGt     <  2:3916612/76‑1 (MQ=255)
                                                                    gcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc    >  2:4610545/1‑76 (MQ=255)
                                                                    gcTCTGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTc    >  2:2161237/1‑76 (MQ=255)
                                                                     ctctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCg   <  1:1114218/76‑1 (MQ=255)
                                                                      tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCgc  >  2:4733721/1‑75 (MQ=255)
                                                                      tctGCTGCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGt  >  1:4656214/1‑76 (MQ=255)
                                                                             |                                                                     
GTCCTGATTGTTGGCAATCATCTCGTTCAGCGCCGATTCAATCTGCGCCAGATCGCTGGAAGTCGCGCGCTCTGCTCCCAACGGGCAATCGCCGGTTCCACCAGATTTCGCACTTCGCTCATGGCACTGATAAGCCGTGGGTCGT  >  NC_000913/3872375‑3872519

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