Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,501,912 T→A 100% I117I (ATT→ATA tehB → tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,501,9120TA97.2% 224.5 / ‑2.0 71I117I (ATT→ATAtehBtellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (36/33);  total (37/34)
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.36e-01

TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008
                                                                                                |                                                                                                
tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               >  1:1262025/1‑100 (MQ=255)
  aTCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                             >  2:1141061/1‑100 (MQ=255)
   tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                             >  1:1184759/1‑99 (MQ=255)
    cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGTGTTGATAGCCAat                                                                                            >  1:727474/1‑99 (MQ=255)
     tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAatat                                                                                          >  2:1071247/1‑100 (MQ=255)
      gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         <  1:2065088/100‑1 (MQ=255)
       aaTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGc                                                                                        <  1:845711/100‑1 (MQ=255)
        aTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAATTATCCCCGGGTTGATAGCCAATATGCa                                                                                       <  2:1805549/100‑1 (MQ=255)
             cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTCAAACCATCCCCGGGTTGATAGCCAATATGCAACg                                                                                    >  2:2089370/1‑98 (MQ=255)
               cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  1:21628/100‑1 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                              <  2:1674308/100‑1 (MQ=255)
                 caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCAc                                                                              <  2:1181149/100‑1 (MQ=255)
                  aTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCATGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACt                                                                             <  2:707794/100‑1 (MQ=255)
                     tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGATGATGGTACTTGAGGCTAACACTATCCCAGGGGTGTTAGAAAAGATGCAAGTTTGAACTaaa                                                                          >  2:639639/1‑100 (MQ=255)
                      gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGAAACGTTACACTAAAc                                                                         >  1:158055/1‑100 (MQ=255)
                       aTAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAGCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAcc                                                                        >  2:1987199/1‑100 (MQ=255)
                         aGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCt                                                                       >  1:629180/1‑99 (MQ=255)
                           aCAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggt                                                                    >  1:1759483/1‑100 (MQ=255)
                            cAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtg                                                                   <  1:926591/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 >  2:1716807/1‑100 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 <  1:295046/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 <  1:2268708/100‑1 (MQ=255)
                              gTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTggtggt                                                                 <  1:700291/100‑1 (MQ=255)
                                aCGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTa                                                               <  2:667709/100‑1 (MQ=255)
                                aCGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTa                                                               <  2:634646/100‑1 (MQ=255)
                                 cGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTAc                                                              <  1:2310839/100‑1 (MQ=255)
                                      ttATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCt                                                         <  2:733085/100‑1 (MQ=255)
                                       tATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTg                                                        <  2:1270099/100‑1 (MQ=255)
                                       tATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTg                                                        >  2:128432/1‑100 (MQ=255)
                                            tttCGACTGTTGTGCTGATGTCCCTTGAGGCAATAACCATCCCCGGGTTTATAGCCAATATGCAACGTTGCACTACACCTGGTCGTTACAACCTAATtgt                                                   >  2:2326264/1‑100 (MQ=255)
                                             ttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtgtg                                                  >  1:510703/1‑100 (MQ=255)
                                               cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCGCTAAACCTGGTGGTTACTACCTGATtgtgtc                                                <  1:1613562/100‑3 (MQ=255)
                                                gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                               >  1:319004/1‑100 (MQ=255)
                                                  cTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             >  2:433847/1‑100 (MQ=255)
                                                  cTGTGGTGCTGATGTTCCCTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             <  2:1722113/100‑1 (MQ=255)
                                                    gtggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGa                                           >  2:2014848/1‑100 (MQ=255)
                                                       gtgCTTATGTTCCTTGAGGCTAAAACCAACCACGGGTTGATAGCCAATATGAAACGTTGCACTAAACCTGGTGGTTACAACCTGATTCTGGCGGCGATg                                         <  1:632700/99‑1 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        >  1:284336/1‑100 (MQ=255)
                                                       gtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATgg                                        <  2:949645/100‑1 (MQ=255)
                                                        tgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGa                                       >  2:538710/1‑100 (MQ=255)
                                                        tgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGa                                       >  2:141237/1‑100 (MQ=255)
                                                         gCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGAt                                      >  1:314166/1‑100 (MQ=255)
                                                          cTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATa                                     >  1:381089/1‑100 (MQ=255)
                                                           tGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAc                                    >  1:463865/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:1531428/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:596093/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  1:2026420/1‑100 (MQ=255)
                                                              tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 >  2:1545359/1‑100 (MQ=255)
                                                               gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCt                                <  1:1271595/100‑1 (MQ=255)
                                                               gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCt                                <  2:1897323/100‑1 (MQ=255)
                                                               gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCt                                <  2:1324740/100‑1 (MQ=255)
                                                               gTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCt                                <  1:1362262/100‑1 (MQ=255)
                                                                ttCCTTGAGGCTAAAGCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                               <  1:1987221/100‑1 (MQ=255)
                                                                  ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGa                              <  2:68561/99‑1 (MQ=255)
                                                                  ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             <  1:1139858/100‑1 (MQ=255)
                                                                   cTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             <  2:1262445/99‑1 (MQ=255)
                                                                       aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                        >  1:2284211/1‑100 (MQ=255)
                                                                        ggCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCGGATTATcc                        >  2:2334899/1‑99 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCCGATTATCCAt                      <  2:1099407/100‑1 (MQ=255)
                                                                         gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      >  1:2184876/1‑100 (MQ=255)
                                                                          cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                     >  1:482185/1‑100 (MQ=255)
                                                                          cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                     <  1:2376765/100‑1 (MQ=255)
                                                                            aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                   >  2:172612/1‑100 (MQ=255)
                                                                            aaaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTa                   <  1:348730/100‑1 (MQ=255)
                                                                                   tCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCgg            >  1:977673/1‑100 (MQ=255)
                                                                                   tCCCCGGGTTGATAGCCAATATGCAAAGTTACACTAAACATGGTGGTTACAAACTGATTGTGGCGGCGATGGATACCGCTAATTATCCATGTACCGTCgg            >  2:686841/1‑100 (MQ=255)
                                                                                    ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           <  2:2284185/100‑1 (MQ=255)
                                                                                      ccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           <  1:48813/98‑1 (MQ=255)
                                                                                           ttGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGt    >  1:2102260/1‑100 (MQ=255)
                                                                                            tGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGtt   <  2:463862/100‑1 (MQ=255)
                                                                                             gATAGCCAATATGCAACGTTGCACAAAACCAGGTGGTTAAAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGttt  >  2:928543/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008

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: