Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,442 C→G 8.6% L792V (CTG→GTG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,4420CG 8.6% 237.1 / 10.5 70L792V (CTG→GTG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (3/3);  total (35/35)
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.76e-01

CTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCT  >  NC_000913/1660301‑1660581
                                                                                                                                             |                                                                                                                                           
cTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                                                                                       >  1:820027/1‑148 (MQ=255)
cTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                                                                                       <  4:128805/148‑1 (MQ=255)
  aCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCCGTGGTGAAg                                                                                                                                     >  1:335171/1‑148 (MQ=255)
   ccAGCCAAAGTAACAACACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCGCCAGTGGTGAAgg                                                                                                                                    <  2:335171/148‑1 (MQ=255)
        cAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAAc                                                                                                                               >  4:154332/1‑148 (MQ=255)
                 accacGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcac                                                                                                                      <  3:154332/148‑1 (MQ=255)
                    acGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACt                                                                                                                   >  3:236670/1‑148 (MQ=255)
                            tcCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGt                                                                                                           <  4:236670/148‑1 (MQ=255)
                             cccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaac                                                                                                                                <  1:533787/126‑2 (MQ=255)
                             cccTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaaa                                                                                                                                >  2:533787/1‑126 (MQ=255)
                                      tAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                     <  3:515503/112‑1 (MQ=255)
                                      tAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                     >  4:515503/1‑112 (MQ=255)
                                         cgcTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                                                                                              <  2:820027/148‑1 (MQ=255)
                                             aTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTg                                                                                                             >  1:269969/1‑129 (MQ=255)
                                             aTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTg                                                                                                             <  2:269969/129‑1 (MQ=255)
                                             aTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaa                                                                                                 <  2:793827/141‑1 (MQ=255)
                                             aTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaa                                                                                                 >  1:793827/1‑141 (MQ=255)
                                                   cAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCa                                                                                    >  2:80721/1‑148 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTg                                                                                                                                           >  1:1056699/1‑91 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTg                                                                                                                                           <  2:1056699/91‑1 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTg                                                                                                                                           >  3:327958/1‑91 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTg                                                                                                                                           <  4:327958/91‑1 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                                                                                              >  3:247288/1‑136 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                                                                                              <  4:247288/136‑1 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                           >  3:470153/1‑139 (MQ=255)
                                                     gggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                           <  4:470153/139‑1 (MQ=255)
                                                      ggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                                                                                                 >  2:735263/1‑100 (MQ=255)
                                                      ggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                                                                                                 <  1:735263/100‑1 (MQ=255)
                                                      ggAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGa                                                                                 >  2:584019/1‑148 (MQ=255)
                                                        aGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                     <  2:665045/110‑1 (MQ=255)
                                                        aGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGcaca                                                                                                                     >  1:665045/1‑110 (MQ=255)
                                                             ggCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGaa                                                                                                                                      <  4:622625/88‑1 (MQ=255)
                                                             ggCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGaa                                                                                                                                      >  3:622625/1‑88 (MQ=255)
                                                              gCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                     >  3:473933/1‑88 (MQ=255)
                                                              gCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAg                                                                                                                                     <  4:473933/88‑1 (MQ=255)
                                                                   aGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcagca                                                                                                                       >  2:826220/1‑97 (MQ=255)
                                                                   aGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCgcagca                                                                                                                       <  1:826220/97‑1 (MQ=255)
                                                                        aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCa                                                               >  2:257333/1‑148 (MQ=255)
                                                                        aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCa                                                               <  1:257333/148‑1 (MQ=255)
                                                                         atatGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCAt                                                              >  4:860586/1‑148 (MQ=255)
                                                                               cTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGa                                                                                                  <  3:722538/106‑1 (MQ=255)
                                                                               cTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGa                                                                                                  >  4:722538/1‑106 (MQ=255)
                                                                                   cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGaa                                                                                                                                      >  3:821593/1‑66 (MQ=255)
                                                                                   cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGaa                                                                                                                                      <  4:821593/66‑1 (MQ=255)
                                                                                   cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                           >  2:308588/1‑109 (MQ=255)
                                                                                   cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                                                                                           <  1:308588/109‑1 (MQ=255)
                                                                                   cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                    <  3:860586/148‑1 (MQ=255)
                                                                                       aaaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                <  1:80721/148‑1 (MQ=255)
                                                                                         aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaa                                              <  1:584019/148‑1 (MQ=255)
                                                                                                cATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                                                                                                >  2:397623/1‑75 (MQ=255)
                                                                                                cATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                                                                                                <  1:397623/75‑1 (MQ=255)
                                                                                                           tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAatg                            >  1:510189/1‑148 (MQ=255)
                                                                                                              gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                >  3:875002/1‑125 (MQ=255)
                                                                                                              gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGa                                                <  4:875002/125‑1 (MQ=255)
                                                                                                                aTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                    >  4:870871/1‑119 (MQ=255)
                                                                                                                aTACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATtttgt                                                    <  3:870871/119‑1 (MQ=255)
                                                                                                                  aCGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGt                                                                                                           >  3:1098324/1‑62 (MQ=255)
                                                                                                                  aCGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGt                                                                                                           <  4:1098324/62‑1 (MQ=255)
                                                                                                                     cTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATcc                  <  2:510189/148‑1 (MQ=255)
                                                                                                                           aCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATcagatcgg            >  1:563993/1‑141 (MQ=255)
                                                                                                                           aCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACca            >  3:731305/1‑148 (MQ=255)
                                                                                                                              ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaaa                                             <  1:418080/112‑1 (MQ=255)
                                                                                                                              ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGaaaa                                             >  2:418080/1‑112 (MQ=255)
                                                                                                                                  tcccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCa                                                                                    >  3:983273/2‑69 (MQ=38)
                                                                                                                                  tcccTTCGCAAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCa                                                                                    <  4:983273/68‑1 (MQ=255)
                                                                                                                                     cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCt  <  4:731305/148‑1 (MQ=255)
                                                                                                                                      ttCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                      <  2:328378/95‑1 (MQ=255)
                                                                                                                                      ttCGCCAGTGGTGAAGGAAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACAtttt                                                      >  1:328378/1‑95 (MQ=255)
                                                                                                                                      ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                   >  2:255621/1‑114 (MQ=255)
                                                                                                                                      ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTg                                   <  1:255621/114‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
CTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCT  >  NC_000913/1660301‑1660581

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