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

GGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAA  >  NC_000913/3014151‑3014393                                                                                                                                         |                                                                                                           ggAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTaa                                                                                                         >  2:82680/1‑139 (MQ=255)       gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                   >  2:623295/1‑139 (MQ=255)       gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                   >  2:233256/1‑139 (MQ=255)             cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                             >  1:1048359/1‑139 (MQ=255)               gCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                           <  2:449735/139‑1 (MQ=255)               gCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcg                                                                                           >  1:340017/1‑139 (MQ=255)                 tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                         <  2:45338/139‑1 (MQ=255)                 tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                         <  2:1237844/139‑1 (MQ=255)                       caaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                   <  2:1211338/139‑1 (MQ=255)                           ttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                               <  1:464555/139‑1 (MQ=255)                            tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                              <  2:1048359/139‑1 (MQ=255)                            tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                              <  2:1133024/139‑1 (MQ=255)                             gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGAt                                                                             >  1:1036175/1‑139 (MQ=255)                              cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            <  2:139191/139‑1 (MQ=255)                              cgcCATAGCAAACGGACAATTGTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCatgatg                                                                            <  1:1175041/139‑1 (MQ=255)                                         aCGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTatat                                                                 <  2:1036175/139‑1 (MQ=255)                                               aaTTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCcca                                                           <  1:233256/139‑1 (MQ=255)                                                             aGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCctct                                             <  1:82680/139‑1 (MQ=255)                                                                           tgtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                               >  2:1330076/1‑139 (MQ=255)                                                                             tgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                             <  1:623295/139‑1 (MQ=255)                                                                                     ttCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                <  2:243799/112‑1 (MQ=255)                                                                                     ttCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                                                >  1:243799/1‑112 (MQ=255)                                                                                       cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                         <  2:936842/69‑1 (MQ=255)                                                                                       cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                         >  1:936842/1‑69 (MQ=255)                                                                                       cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                             >  1:8357/1‑129 (MQ=255)                                                                                       cATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                             <  2:8357/129‑1 (MQ=255)                                                                                        aTACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAgg                  >  1:334817/1‑139 (MQ=255)                                                                                         tACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                              >  1:575449/1‑78 (MQ=255)                                                                                         tACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCAAAAATCTGCGCGCGATTGCAGGa                                                                              <  2:575449/78‑1 (MQ=255)                                                                                              aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGAGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCg            >  2:914611/1‑139 (MQ=255)                                                                                               ccTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcgc                                                                                        >  2:573829/1‑62 (MQ=255)                                                                                               ccTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcgc                                                                                        <  1:573829/62‑1 (MQ=255)                                                                                               ccTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGc           >  1:863981/1‑139 (MQ=255)                                                                                                  cAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                          >  1:46098/1‑41 (MQ=255)                                                                                                  cAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                          <  2:46098/41‑1 (MQ=255)                                                                                                     tGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCACCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCt              <  1:452010/130‑1 (MQ=255)                                                                                                     tGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCACCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCt              >  2:452010/1‑130 (MQ=255)                                                                                                         atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                             <  1:510690/108‑1 (MQ=255)                                                                                                         atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCa                             >  2:510690/5‑112 (MQ=255)                                                                                                         atatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCaa  >  1:699157/5‑139 (MQ=255)                                                                                                          tatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGAtt                                                                                    <  1:1351132/53‑1 (MQ=255)                                                                                                          tatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGAtt                                                                                    >  2:1351132/4‑56 (MQ=255)                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                   <  1:1234670/54‑1 (MQ=255)                                                                                                           atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                   >  2:1234670/3‑56 (MQ=255)                                                                                                              tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                         <  2:1226282/47‑1 (MQ=255)                                                                                                              tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                         >  1:1226282/1‑47 (MQ=255)                                                                                                                      aCCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            <  2:443398/52‑1 (MQ=255)                                                                                                                      aCCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            >  1:443398/1‑52 (MQ=255)                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            <  2:605219/51‑1 (MQ=255)                                                                                                                       cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            >  1:605219/1‑51 (MQ=255)                                                                                                                         cGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            <  2:1212631/49‑1 (MQ=255)                                                                                                                         cGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                            >  1:1212631/1‑49 (MQ=255)                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                               >  2:741732/1‑37 (MQ=255)                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                               <  1:741732/37‑1 (MQ=38)                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc               >  2:797272/1‑101 (MQ=255)                                                                                                                                  cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc               <  1:797272/101‑1 (MQ=255)                                                                                                                                         |                                                                                                           GGAGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAA  >  NC_000913/3014151‑3014393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑