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_000913	1,501,912	0	T	A	100.0%	258.3 / ‑6.7	80	I117I (ATT→ATA)	tehB	tellurite, selenium methyltransferase, SAM‑dependent; tellurite, selenium resistance protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (36/42);  minor base G (1/0);  total (38/42)

TGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATTGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGTTT  >  NC_000913/1501816‑1502008                                                                                                 |                                                                                                 tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               <  2:1569442/100‑1 (MQ=255) tGATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               >  2:807498/1‑100 (MQ=255)  gATCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCa                                                                                              >  1:1312874/1‑100 (MQ=255)   aTCTGAATAACCTCATATTTGATAGACAGTACTATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGcc                                                                                               <  2:2377102/98‑1 (MQ=255)   aTCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                             <  2:123493/100‑1 (MQ=255)    tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCaa                                                                                             <  2:207228/99‑1 (MQ=255)    tCTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAat                                                                                            <  2:2484066/100‑1 (MQ=255)     cTGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAata                                                                                           <  2:2042527/100‑1 (MQ=255)      tGAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAatat                                                                                          >  1:2473409/1‑100 (MQ=255)       gAATAACCTCACATTTGATAGACAGTACGCTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGCTAGCCAata                                                                                           >  2:2467468/1‑98 (MQ=255)       gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         <  2:236065/100‑1 (MQ=255)       gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         >  2:1880855/1‑100 (MQ=255)       gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         >  2:2132142/1‑100 (MQ=255)       gAATAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         >  1:1121793/1‑100 (MQ=255)       gAATAACCTCACATTTGATAGACAGTACGATTGTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATg                                                                                         >  2:2318773/1‑100 (MQ=255)        aaTAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGc                                                                                        <  1:103546/100‑1 (MQ=255)          tAACCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCaa                                                                                      <  1:748922/100‑1 (MQ=255)            aCCTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACg                                                                                    >  1:304043/1‑100 (MQ=255)              cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGtt                                                                                  <  2:687144/100‑1 (MQ=255)              cTCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGt                                                                                   >  1:2329459/1‑99 (MQ=255)               tCACATTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTg                                                                                 <  1:639457/100‑1 (MQ=255)                cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  2:2283992/100‑1 (MQ=255)                cacaTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGc                                                                                <  2:692239/100‑1 (MQ=255)                  caTTTGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCan                                                                              <  2:1444840/100‑2 (MQ=255)                      tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAGCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTaaa                                                                          <  2:1761469/100‑1 (MQ=255)                      tGATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACAATCCCCGGGTTAATAGCCTATATGCAACGTTGTACTaa                                                                           >  2:1217890/1‑99 (MQ=255)                       gATAGACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAAc                                                                         >  1:2347119/1‑100 (MQ=255)                           gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                     >  1:727628/1‑100 (MQ=255)                           gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                     >  2:2221406/1‑100 (MQ=255)                           gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                     <  2:586756/100‑1 (MQ=255)                           gACAGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtgg                                                                     >  1:2127352/1‑100 (MQ=255)                              aGTACGATTTTATTCTTTCGACTGTGGTGCTGATGTTCCATGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCtggtgg                                                                  <  2:57514/100‑1 (MQ=255)                                   gATTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACa                                                             <  2:299603/100‑1 (MQ=255)                                    aTTTTATTCTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACaa                                                            >  1:1827706/1‑100 (MQ=255)                                            cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    >  1:1695408/1‑100 (MQ=255)                                            cTTTCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    >  1:2479380/1‑100 (MQ=255)                                            cTTTCGACTGTGGGGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtg                                                    <  2:2156094/100‑1 (MQ=255)                                             tttCGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATtgt                                                   >  2:1000213/1‑100 (MQ=255)                                                cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                                >  2:1851589/1‑100 (MQ=255)                                                cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                                <  1:964040/100‑1 (MQ=255)                                                cGACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAAAGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggc                                                <  2:1362210/100‑1 (MQ=255)                                                 gACTGTGGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcg                                               <  1:2129176/100‑1 (MQ=255)                                                 gACTGTGGAGAAGATGTTCCTTGAGGATAAAACCATAGCCAGGTTGAACGCCAATATACAACGTTGCACTAAACTTGGTGGTTACAAGCTGATTGAggat                                               >  2:522425/1‑98 (MQ=255)                                                  aCTGTGGTGCTGATGTTCCTTGAGGCCAAAACCATCCCCAGGTGGATAGCCAATATGCAACGTTGCACTAAACCTGCTGGTTACAACCTGATTGTggcgg                                              <  2:1992087/100‑1 (MQ=255)                                                    tgtgGTGCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             >  1:1523532/1‑99 (MQ=255)                                                     gtggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTggcggc                                             <  2:798349/98‑1 (MQ=255)                                                     gtggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGa                                           >  2:1715751/1‑100 (MQ=255)                                                       ggtgCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATg                                         <  2:1477311/100‑1 (MQ=255)                                                          gCTGATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGAt                                      >  1:991444/1‑100 (MQ=255)                                                             gATGTTCTTTGAGGCTAAAACCATCCCCTGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGCTGGTTACAACCTGATTGTGGCGGCGATGGATAcc                                   >  2:1345055/1‑100 (MQ=255)                                                             gATGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATAcc                                   >  1:867175/1‑100 (MQ=255)                                                               tGTTCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGc                                 <  1:1412212/100‑1 (MQ=255)                                                                 ttCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTg                               >  2:1778561/1‑100 (MQ=255)                                                                  tCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGa                              >  1:909630/1‑100 (MQ=255)                                                                  tCCTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGa                              >  1:292138/1‑100 (MQ=255)                                                                   ccTTGAGGCTAAAACCATCCCCGGGTTGATGGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             >  1:1051133/1‑100 (MQ=255)                                                                   ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             <  2:1854108/100‑1 (MQ=255)                                                                   ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             <  2:292134/100‑1 (MQ=255)                                                                   ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             <  1:916544/100‑1 (MQ=255)                                                                   ccTTGAGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGAt                             >  2:734790/1‑100 (MQ=255)                                                                        aGGCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATcc                        <  1:2176387/100‑1 (MQ=255)                                                                          tcTAAAACCATCCCCGGGTTGATAGCCAATATGGAACGTTGCGCTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATGATCCa                       <  2:2051301/98‑1 (MQ=255)                                                                          gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      >  2:395067/1‑100 (MQ=255)                                                                          gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      >  2:947895/1‑100 (MQ=255)                                                                          gCTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCAt                      >  1:1830922/1‑100 (MQ=255)                                                                           cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                     <  2:1444371/100‑1 (MQ=255)                                                                           cTAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATg                     >  1:515551/1‑100 (MQ=255)                                                                            tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                    >  1:674647/1‑100 (MQ=255)                                                                            tAAAACCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGt                    <  1:107831/100‑1 (MQ=255)                                                                              aaaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAc                  <  1:151772/100‑1 (MQ=255)                                                                               aaCCATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTAcc                 <  2:2501420/100‑1 (MQ=255)                                                                                 ccATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGt               <  2:1322354/100‑1 (MQ=255)                                                                                  cATCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTc              >  2:1353997/1‑100 (MQ=255)                                                                                   aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             <  2:1065026/100‑1 (MQ=255)                                                                                   aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             <  1:807505/100‑1 (MQ=255)                                                                                   aTCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCg             <  1:2425211/100‑1 (MQ=255)                                                                                    tCCCCGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCgg            <  1:947903/100‑1 (MQ=255)                                                                                     ccccGGGTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGc           <  1:1523989/100‑1 (MQ=255)                                                                                         gggTTGATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTcc       <  1:2455164/100‑1 (MQ=255)                                                                                              gATAGCCAATATGCAACGTTGCACTAAACCTGGTGGTTACAACCTGATTGTGGCGGCGATGGATACCGCTGATTATCCATGTACCGTCGGCTTCCCGttt  >  1:948359/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: ‑