Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,014,287	C→T	*257* (TAG→TAA)	yqeC ←	putative selenium‑dependent hydroxylase accessory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,014,287	0	C	T	100.0%	235.4 / NA	73	*257* (TAG→TAA)	yqeC	putative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/37);  total (36/37)

GGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGTT  >  NC_000913/3014150‑3014421                                                                                                                                          |                                                                                                                                       gggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                      <  1:211479/139‑1 (MQ=255)        gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                                               >  2:558191/1‑139 (MQ=255)        gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                                               >  1:1550268/1‑139 (MQ=255)              cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                                                         <  2:1550268/139‑1 (MQ=255)                gCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                                                       <  2:799864/139‑1 (MQ=255)                           tttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGNNNACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAg                                                                                                            >  2:542867/1‑139 (MQ=255)                             tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          <  2:80030/139‑1 (MQ=255)                              gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                                                         <  2:464557/139‑1 (MQ=255)                               cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                        <  2:991097/139‑1 (MQ=255)                               cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                        <  1:1636109/139‑1 (MQ=255)                                gcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGt                                                                                                       <  1:558191/139‑1 (MQ=255)                                   aTAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCt                                                                                                    >  2:652170/1‑139 (MQ=255)                                    tAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                                                                   <  1:652170/139‑1 (MQ=255)                                        aaaCGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                                               <  2:1682847/139‑1 (MQ=255)                                               cAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                           <  1:1747843/120‑1 (MQ=255)                                               cAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                           >  2:1747843/1‑120 (MQ=255)                                                                    gCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                                   <  1:1316623/139‑1 (MQ=255)                                                                          tctGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgttgt                                                             >  1:995827/1‑139 (MQ=255)                                                                            tgtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                                           >  1:1504875/1‑139 (MQ=255)                                                                             gtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGc                                                          >  2:1703349/1‑139 (MQ=255)                                                                                aaaaTATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGc                                                       <  2:543783/139‑1 (MQ=255)                                                                                      ttCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          >  1:554180/1‑82 (MQ=255)                                                                                      ttCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          <  2:554180/82‑1 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                      >  1:1612672/1‑52 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                      <  1:806547/52‑1 (MQ=255)                                                                                              aaCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                      <  2:579972/46‑1 (MQ=255)                                                                                              aaCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                                      >  1:579972/1‑46 (MQ=38)                                                                                               aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                                           <  2:414300/120‑1 (MQ=255)                                                                                               aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg                                        >  1:930535/1‑139 (MQ=255)                                                                                               aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGACAGATTGCCTCTACGTCGTAttgttg                                                           >  1:414300/1‑120 (MQ=255)                                                                                                ccTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGtt                                                                                                      <  2:921166/76‑1 (MQ=255)                                                                                                ccTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGtt                                                                                                      >  1:921166/1‑76 (MQ=255)                                                                                                          atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          >  1:429923/5‑63 (MQ=255)                                                                                                          atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          <  2:429923/59‑1 (MQ=255)                                                                                                           tatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                        <  1:1613425/61‑1 (MQ=255)                                                                                                           tatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                        >  2:1613425/4‑64 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGTGCCGTTAATTCACAAATCTGCGCGCGATTGCATGATGTTCTTGTATATCGCCCAGCCAGATTACCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgtg                                 <  1:1774260/132‑1 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                        >  1:1551911/3‑63 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                                        <  2:1551911/61‑1 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCGCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                         <  1:676596/108‑1 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                         >  2:676596/3‑110 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                           <  2:1430523/122‑1 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                           >  1:1430523/3‑124 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTtgtg                                 >  2:1774260/3‑134 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc                            >  2:768513/3‑139 (MQ=255)                                                                                                             tcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                           <  2:995827/138‑1 (MQ=255)                                                                                                               tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                                      <  2:1084055/94‑1 (MQ=255)                                                                                                               tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                                      >  1:1084055/1‑94 (MQ=255)                                                                                                               tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                         >  1:53800/1‑139 (MQ=255)                                                                                                               tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                         >  2:552247/1‑139 (MQ=255)                                                                                                                   aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                                                                   >  2:1586022/1‑61 (MQ=255)                                                                                                                   aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCtt                                                                                                   <  1:1586022/61‑1 (MQ=255)                                                                                                                        cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGAt                                                                                                                 >  2:78140/1‑42 (MQ=255)                                                                                                                        cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGAt                                                                                                                 <  1:78140/42‑1 (MQ=255)                                                                                                                        cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGaa                                                <  1:436522/107‑1 (MQ=255)                                                                                                                        cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGaa                                                >  2:436522/1‑107 (MQ=255)                                                                                                                             aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                      <  2:541419/112‑1 (MQ=255)                                                                                                                             aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCt                                      >  1:541419/1‑112 (MQ=255)                                                                                                                                  ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGt                                                    >  2:1300012/1‑93 (MQ=255)                                                                                                                                  ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGt                                                    <  1:1300012/93‑1 (MQ=255)                                                                                                                                  ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                                            <  1:247889/101‑1 (MQ=255)                                                                                                                                  ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                                            >  2:247889/1‑101 (MQ=255)                                                                                                                                  ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                         >  2:1274393/1‑120 (MQ=255)                                                                                                                                  ccGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa                         <  1:1274393/120‑1 (MQ=255)                                                                                                                                   cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          <  1:959819/38‑1 (MQ=255)                                                                                                                                   cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                                          >  2:959819/1‑38 (MQ=38)                                                                                                                                   cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCgtg                                                                 <  1:382203/79‑1 (MQ=255)                                                                                                                                   cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCgtg                                                                 >  2:382203/1‑79 (MQ=255)                                                                                                                                   cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                         <  1:948357/87‑1 (MQ=255)                                                                                                                                   cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                                                         >  2:948357/1‑87 (MQ=255)                                                                                                                                    gggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCTGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGaaaaa       <  1:1283637/136‑1 (MQ=255)                                                                                                                                    gggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGaaaaa       >  2:1283637/1‑136 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGtt  >  1:536358/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                       GGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGAGAAAAACGGTT  >  NC_000913/3014150‑3014421

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: ‑