Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,657,898 A→T 24.7% intergenic (+28/‑171) ynfD → / → ynfE DUF1161 family periplasmic protein/putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,657,8980AT24.7% 146.1 / 54.5 85intergenic (+28/‑171)ynfD/ynfEDUF1161 family periplasmic protein/putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base T (10/11);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

CCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACT  >  NC_000913/1657803‑1657996
                                                                                               |                                                                                                  
ccaccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACaaa                                                                                                <  1:1981839/100‑1 (MQ=255)
   ccaGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGCAACACACTATCGACTGAACGCCGGATATGACAAATcc                                                                                             >  1:2194944/1‑100 (MQ=255)
    caGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                                                            <  1:1016038/100‑1 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGTTGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          <  1:1127528/100‑1 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                          <  1:105191/100‑1 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                          <  1:922805/100‑1 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGc                                                                                          <  2:2561282/100‑1 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          >  2:1094316/1‑100 (MQ=255)
      gCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGc                                                                                          <  1:2417611/100‑1 (MQ=255)
       cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                         <  2:2074210/100‑1 (MQ=255)
       cGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCg                                                                                         <  2:2105591/100‑1 (MQ=255)
        ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                        <  1:1449607/100‑2 (MQ=255)
        ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                        <  1:1751393/100‑2 (MQ=255)
        ggTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgt                                                                                        <  1:903687/100‑2 (MQ=255)
             cGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTACCACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagt                                                                                   <  2:1073641/100‑7 (MQ=37)
             cGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagt                                                                                   <  2:1761173/100‑7 (MQ=38)
              gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCGGCgttcagtc                                                                                  <  2:1045060/100‑8 (MQ=37)
              gTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGaa                                                                                  >  2:749899/1‑100 (MQ=255)
               tctcTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                 <  1:895122/100‑1 (MQ=255)
               tcgctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATCATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                                                 <  1:1986014/97‑1 (MQ=255)
                  ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa                                                                              <  2:2495847/100‑1 (MQ=255)
                  ctGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATa                                                                              >  1:921602/1‑100 (MQ=255)
                    gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                            >  2:769809/1‑100 (MQ=255)
                    gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                            <  2:324441/100‑1 (MQ=255)
                    gCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATaca                                                                            <  1:1932608/100‑1 (MQ=255)
                     cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                           >  2:136885/1‑100 (MQ=255)
                     cTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaa                                                                           >  2:1545490/1‑100 (MQ=255)
                        ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                                                        >  2:183234/1‑100 (MQ=255)
                               gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                 >  2:367666/1‑100 (MQ=255)
                               gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                 <  2:1094029/100‑1 (MQ=255)
                               gTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCa                                                                 >  2:770858/1‑100 (MQ=255)
                                tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATCATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCAc                                                                <  1:626069/100‑1 (MQ=255)
                                    gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTCAACATACAACATAATACCACCtta                                                            >  1:207042/1‑100 (MQ=255)
                                      cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                                                          >  2:600472/1‑100 (MQ=255)
                                      cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                                                          >  1:347402/1‑100 (MQ=255)
                                      cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCAACCttagt                                                          >  1:1566596/1‑98 (MQ=255)
                                        ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTAc                                                        <  2:1850711/100‑1 (MQ=255)
                                         gATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                                                       <  1:1545481/100‑1 (MQ=255)
                                             gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                                                    <  1:656282/99‑1 (MQ=255)
                                               gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAcc                                                 <  1:600469/100‑1 (MQ=255)
                                                 gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCt                                               <  1:2347601/100‑1 (MQ=255)
                                                  ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                                              >  2:185564/1‑100 (MQ=255)
                                                  ccTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                                              >  1:2158787/1‑100 (MQ=255)
                                                     gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACCAATCCGGCGATTTCAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa                                           >  1:611210/1‑100 (MQ=255)
                                                       cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                         <  2:704208/100‑1 (MQ=255)
                                                       cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                         >  2:94994/1‑100 (MQ=255)
                                                       cGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt                                         >  1:1911811/1‑100 (MQ=255)
                                                         aaCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGa                                       <  2:1112100/100‑1 (MQ=255)
                                                          aCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat                                      <  1:1888856/100‑1 (MQ=255)
                                                             gCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg                                   <  1:1932318/100‑1 (MQ=255)
                                                                gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACCCATACCCTTCTATTGATATGGAt                                >  2:1225347/1‑100 (MQ=255)
                                                                   acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaa                             >  1:703048/1‑100 (MQ=255)
                                                                    cacaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat                            <  2:2187283/100‑1 (MQ=255)
                                                                    cacaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat                            >  2:1971384/1‑100 (MQ=255)
                                                                     tccgatctgACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata                           <  1:582861/92‑1 (MQ=255)
                                                                     acaTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaata                           >  1:1563576/1‑100 (MQ=255)
                                                                      caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaat                            >  1:1068428/1‑98 (MQ=255)
                                                                          aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                      <  2:1147966/100‑1 (MQ=255)
                                                                          aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                      >  1:2445232/1‑100 (MQ=255)
                                                                          aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCt                      <  2:1459353/100‑1 (MQ=255)
                                                                           tCGACTGAACGCCGGATATGACAAATCCGGCGCTTTGAACATACTACATAATCCCACCTTATTACTCATACCCTTCTATTGCTCTGGATTAATAATTCt                      <  1:2064472/99‑1 (MQ=255)
                                                                             ttcaaatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTaa                   >  1:1274162/8‑100 (MQ=37)
                                                                             gACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATtaataagatcgga                   >  2:582863/1‑93 (MQ=255)
                                                                              aCTGAACGCCGGATATGACAAATCCGGCGATTTGAGCATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAc                  <  1:109281/100‑1 (MQ=255)
                                                                               cTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAAcc                 <  2:1390109/100‑1 (MQ=255)
                                                                                 gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAGCCTTATATTGATATGGCTTAATATTTCTTAACCCa               >  1:2151422/1‑100 (MQ=255)
                                                                                 gAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa               >  1:2432254/1‑100 (MQ=255)
                                                                                 gAACGCCGGATATGACAAATCCGGCGATGTGAACATACAACATAATTCCACTTTATTACTCATACCCTTCTATAGCGATGGATTAATAATTCTTAAcaca               >  1:1936160/1‑100 (MQ=255)
                                                                                 aatcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCa               >  1:1851047/4‑100 (MQ=38)
                                                                                  atcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATACTTCTTAACCCaa              >  1:472460/3‑100 (MQ=255)
                                                                                   tcGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaa             >  2:414998/2‑100 (MQ=255)
                                                                                    cGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATGCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa            >  1:400101/1‑100 (MQ=255)
                                                                                    cGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa            >  1:34191/1‑100 (MQ=255)
                                                                                    cGCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCaaaa            >  2:460739/1‑100 (MQ=255)
                                                                                     gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt           >  1:954254/1‑100 (MQ=255)
                                                                                     gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt           >  1:939767/1‑100 (MQ=255)
                                                                                     gCCGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCAACTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt           >  1:1438346/1‑100 (MQ=255)
                                                                                     gCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAAt           >  2:1711627/1‑100 (MQ=255)
                                                                                         gATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATggct       >  2:849198/1‑98 (MQ=255)
                                                                                             tGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGAc   >  2:882281/1‑100 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTa      >  1:2023675/1‑96 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  >  2:1963865/1‑100 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  >  2:1589962/1‑100 (MQ=255)
                                                                                              gACAAATCCGGCGATTTGAACATACAACATAATCCCACCATATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  <  1:1539732/100‑1 (MQ=255)
                                                                                              gACAAATCAGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACt  >  1:2445986/1‑100 (MQ=255)
                                                                                               |                                                                                                  
CCACCAGCGGTAACGTCTCTGCTCCCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACT  >  NC_000913/1657803‑1657996

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑