Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,662,154 A→G A533A (GCA→GCG ynfF → putative selenate reductase YnfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,662,1540AG98.0% 171.6 / ‑4.1 50A533A (GCA→GCGynfFputative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/25);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCA  >  NC_000913/1662017‑1662290
                                                                                                                                         |                                                                                                                                        
gTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTTGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCggg                                                                                                                                        <  2:984206/140‑1 (MQ=255)
gTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCggg                                                                                                                                        >  2:604310/1‑140 (MQ=255)
 tGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCgggg                                                                                                                                       >  1:734840/1‑140 (MQ=255)
 tGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCgggg                                                                                                                                       <  1:494756/140‑1 (MQ=255)
  gCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGa                                                                                                                                      <  2:734840/140‑1 (MQ=255)
   cTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGaa                                                                                                                                     <  1:391929/140‑1 (MQ=255)
   cTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCACTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGaa                                                                                                                                     >  2:521245/1‑140 (MQ=255)
     tCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAata                                                                                                                                   <  1:513852/140‑1 (MQ=255)
       aGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATg                                                                                                                                 <  1:451896/140‑1 (MQ=255)
          aTGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGc                                                                                                                              <  1:1047821/140‑1 (MQ=255)
             aCAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTAt                                                                                                                           >  2:100034/1‑140 (MQ=255)
              cAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTAtg                                                                                                                          >  1:845805/1‑140 (MQ=255)
               aGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTAtgt                                                                                                                         >  1:179614/1‑140 (MQ=255)
                         agaTGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGcc                                                                                                               <  1:969663/140‑1 (MQ=255)
                          gatgatTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGccc                                                                                                              <  1:797846/140‑1 (MQ=255)
                          gatgatTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGccc                                                                                                              <  1:506023/140‑1 (MQ=255)
                              atTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAAcc                                                                                                          >  2:337483/1‑140 (MQ=255)
                               ttgttgGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAAccc                                                                                                         >  1:903612/1‑140 (MQ=255)
                                     gCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAAcc                                                                                                   <  2:179614/140‑1 (MQ=255)
                                       aTTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTc                                                                                                 <  1:337483/140‑1 (MQ=255)
                                         tGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAg                                                                                               >  2:180916/1‑140 (MQ=255)
                                           aCCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCa                                                                                             >  1:357396/1‑140 (MQ=255)
                                                 tCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAAttt                                                                                       <  1:92130/140‑1 (MQ=255)
                                                        ggCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCa                                                                                <  2:903612/140‑1 (MQ=255)
                                                          ccTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCaaa                                                                              <  1:180916/140‑1 (MQ=255)
                                                          ccTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCaaa                                                                              >  1:454331/1‑140 (MQ=255)
                                                          ccTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCaaa                                                                              >  2:704477/1‑140 (MQ=255)
                                                          ccTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCaaa                                                                              >  2:404596/1‑140 (MQ=255)
                                                             cGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAccc                                                                           <  1:404596/140‑1 (MQ=255)
                                                               gCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCAt                                                                         <  1:347155/140‑1 (MQ=255)
                                                               gCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCAt                                                                         <  2:403365/140‑1 (MQ=255)
                                                                    gTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACt                                                                    <  2:845805/140‑1 (MQ=255)
                                                                                 cTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACt                                                                    <  2:588322/127‑1 (MQ=255)
                                                                                 cTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACt                                                                    >  1:588322/1‑127 (MQ=255)
                                                                                         cGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCa                                               >  2:610254/1‑140 (MQ=255)
                                                                                         cGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCa                                               <  1:610254/140‑1 (MQ=255)
                                                                                         cGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCa                                               >  2:98327/1‑140 (MQ=255)
                                                                                              tGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTg                                                        >  2:786934/1‑126 (MQ=255)
                                                                                              tGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTg                                                        <  1:786934/126‑1 (MQ=255)
                                                                                                          agCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGAc                              >  2:782462/1‑140 (MQ=255)
                                                                                                          agCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGAc                              <  1:356803/140‑1 (MQ=255)
                                                                                                           gCAAGAAGACCTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACg                             <  1:704477/140‑1 (MQ=255)
                                                                                                                    ccTTATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCaaa                    >  1:79576/1‑140 (MQ=255)
                                                                                                                       tATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCa                      >  2:1033029/1‑135 (MQ=255)
                                                                                                                       tATCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCa                      <  1:1033029/135‑1 (MQ=255)
                                                                                                                        aTCTCTCATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCtt                >  2:421598/1‑140 (MQ=255)
                                                                                                                         tctctcATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCttt               >  2:585195/1‑140 (MQ=255)
                                                                                                                             tcATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTg           >  1:115090/1‑140 (MQ=255)
                                                                                                                              cATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGa          >  2:145788/1‑140 (MQ=255)
                                                                                                                              cATGAATCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGa          <  1:100034/140‑1 (MQ=255)
                                                                                                                              cATGAATCTGCGGGGAATATGGGCTATGTGACCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGa          >  1:863815/1‑140 (MQ=255)
                                                                                                                                    tCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCg    >  2:442798/1‑140 (MQ=255)
                                                                                                                                    tCTGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATACTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCg    <  2:454331/140‑1 (MQ=255)
                                                                                                                                      tGCGGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCa  <  1:782462/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCA  >  NC_000913/1662017‑1662290

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