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%	219.6 / NA	68	*257* (TAG→TAA)	yqeC	putative selenium‑dependent hydroxylase accessory protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/35);  total (33/35)

GGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCG  >  NC_000913/3014150‑3014396                                                                                                                                          |                                                                                                              gggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACTATTATGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCGCATTGGATAT‑CTTTCAGTAACCCGGAGTACTCGTGCCGTTa                                                                                                             <  1:259424/139‑1 (MQ=255)        gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                      >  2:594521/1‑139 (MQ=255)                gCTATCAACAATTTGCGCCGTAGCAAACGGACAATCTTGCTCTTCGGGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTCACCCGGAATACCCGGGCCGTTCATTCATAAATCTgca                                                                                              >  1:285679/1‑138 (MQ=255)                gCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                              <  1:845169/139‑1 (MQ=255)                            ttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                  <  1:439591/139‑1 (MQ=255)                             tGCGCCATAGCAAATGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                 <  1:124465/139‑1 (MQ=255)                             tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                  <  2:1531484/138‑1 (MQ=255)                             tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                  >  1:1531484/1‑138 (MQ=255)                             tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                 <  2:1531283/139‑1 (MQ=255)                               cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                               <  1:773554/139‑1 (MQ=255)                               cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                               <  1:166700/139‑1 (MQ=255)                                gcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                            <  1:852882/125‑1 (MQ=255)                                gcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                            >  2:852882/1‑125 (MQ=255)                                gcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGt                                                                              >  2:311827/1‑139 (MQ=255)                                     aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                 >  2:495986/1‑131 (MQ=255)                                     aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                 <  1:495986/131‑1 (MQ=255)                                     aGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTg                                                                         >  1:920531/1‑139 (MQ=255)                                       cAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGta                                                                       >  1:1085568/1‑139 (MQ=255)                                            ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                  <  2:911994/139‑1 (MQ=255)                                            ggACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                  <  1:582596/139‑1 (MQ=255)                                              aCAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCgac                                                                >  1:324601/1‑137 (MQ=255)                                                aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCcca                                                              <  2:527158/139‑1 (MQ=255)                                                           tCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                   <  1:311827/139‑1 (MQ=255)                                                            cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                      <  2:216939/103‑1 (MQ=255)                                                            cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                      >  1:216939/1‑103 (MQ=255)                                                            cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                  >  1:14981/1‑139 (MQ=255)                                                            cGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                  >  1:1066675/1‑139 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTACTTCACAAATCTGCGCGCGATTGCa                                                                                    >  1:932667/1‑77 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                    <  2:932667/77‑1 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                    <  1:897597/77‑1 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCa                                                                                    >  2:897597/1‑77 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                       <  1:1237860/106‑1 (MQ=255)                                                                                        cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                       >  2:1237860/1‑106 (MQ=255)                                                                                           aCCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                   <  2:920531/139‑1 (MQ=255)                                                                                            ccAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                            <  1:241768/65‑1 (MQ=255)                                                                                            ccAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                            >  2:241768/1‑65 (MQ=255)                                                                                              aaCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                               <  1:906562/76‑1 (MQ=255)                                                                                              aaCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                               >  2:906562/1‑76 (MQ=255)                                                                                               aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg               >  2:604708/1‑139 (MQ=255)                                                                                               aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg               <  2:14981/139‑1 (MQ=255)                                                                                                ccTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGc              >  2:1154937/1‑139 (MQ=255)                                                                                                    aaTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                              >  1:88723/1‑55 (MQ=255)                                                                                                    aaTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                              <  2:88723/55‑1 (MQ=255)                                                                                                     aTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGa                        <  2:311586/124‑1 (MQ=255)                                                                                                     aTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGa                        >  1:311586/1‑124 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                  <  1:848170/90‑1 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCct                                                  >  2:848170/3‑92 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc   >  2:167568/3‑139 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc   >  1:335059/3‑139 (MQ=255)                                                                                                            atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc   >  2:1255371/3‑139 (MQ=255)                                                                                                             tcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg  <  1:167568/138‑1 (MQ=255)                                                                                                                   aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                            <  1:561938/43‑1 (MQ=255)                                                                                                                   aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                            >  2:561938/1‑43 (MQ=255)                                                                                                                        cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                       >  2:1563316/1‑75 (MQ=255)                                                                                                                        cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                       <  1:1563316/75‑1 (MQ=255)                                                                                                                             aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                      <  1:473715/39‑1 (MQ=38)                                                                                                                             aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                      >  2:473715/1‑39 (MQ=255)                                                                                                                             aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg  <  2:552213/123‑1 (MQ=255)                                                                                                                             aaTACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg  >  1:552213/1‑123 (MQ=255)                                                                                                                                     ggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                  <  2:353174/51‑1 (MQ=255)                                                                                                                                     ggCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                  >  1:353174/1‑51 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                  >  1:1465703/1‑50 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                  <  2:1465703/50‑1 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                   <  2:683250/97‑1 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAg                   >  1:683250/1‑97 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg               >  2:1409401/1‑101 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg               <  1:1409401/101‑1 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAAt    >  2:1414538/1‑112 (MQ=255)                                                                                                                                      gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAAt    <  1:1414538/112‑1 (MQ=255)                                                                                                                                          |                                                                                                              GGGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCG  >  NC_000913/3014150‑3014396

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