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

AGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCT  >  NC_000913/3014153‑3014402                                                                                                                                       |                                                                                                                    aGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAAtt                                                                                                                >  2:533383/1‑139 (MQ=255)   tCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTca                                                                                                              >  1:728664/1‑139 (MQ=255)          tCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCt                                                                                                       >  2:1043042/1‑139 (MQ=255)               tATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                  <  2:843724/139‑1 (MQ=255)                    acaaTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGAtt                                                                                             <  2:1044650/139‑1 (MQ=255)                          tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                       <  1:446864/139‑1 (MQ=255)                           gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                       <  2:1062930/138‑1 (MQ=255)                           gcgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                       >  1:1062930/1‑138 (MQ=255)                            cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     <  2:728664/139‑1 (MQ=255)                                     aaaCGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                                                            <  2:1083134/139‑1 (MQ=255)                                               ttttGCTCTTCGAGTTTTGCAGTCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                                  <  2:152667/139‑1 (MQ=255)                                                    ctctTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     <  2:919563/115‑1 (MQ=255)                                                    ctctTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     >  1:919563/1‑115 (MQ=255)                                                               ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     <  1:1189518/104‑1 (MQ=255)                                                               ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     <  2:423092/104‑1 (MQ=255)                                                               ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     >  2:1189518/1‑104 (MQ=255)                                                               ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     >  1:423092/1‑104 (MQ=255)                                                                 gCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                >  1:977719/1‑139 (MQ=255)                                                                 gCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                >  2:136718/1‑139 (MQ=255)                                                                      ctctGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCctct                                                      >  2:734716/1‑128 (MQ=255)                                                                      ctctGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCctct                                                      <  1:734716/128‑1 (MQ=255)                                                                         tgtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                  <  1:920091/81‑1 (MQ=255)                                                                         tgtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgcgcg                                                                                                  >  2:920091/1‑81 (MQ=255)                                                                         tgtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                        >  1:1234106/1‑139 (MQ=255)                                                                           tgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGTCTCTACGTCGTGTTGTTGCa                                      <  1:51661/139‑1 (MQ=255)                                                                            gAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAg                                     <  1:136718/139‑1 (MQ=255)                                                                               aaTATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                        <  1:1056567/101‑1 (MQ=255)                                                                               aaTATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCg                                                                        >  2:1056567/1‑101 (MQ=255)                                                                                  atTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     >  1:576337/1‑85 (MQ=255)                                                                                  atTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     <  2:576337/85‑1 (MQ=255)                                                                                   ttCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                        >  2:31508/1‑81 (MQ=255)                                                                                   ttCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                        <  1:31508/81‑1 (MQ=255)                                                                                        aCCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                   >  1:1052560/1‑49 (MQ=255)                                                                                        aCCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                   <  2:1052560/49‑1 (MQ=255)                                                                                           aaCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                   >  1:664017/1‑46 (MQ=38)                                                                                           aaCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                   <  2:664017/46‑1 (MQ=255)                                                                                                  aTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     <  2:160274/69‑1 (MQ=255)                                                                                                  aTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                     >  1:160274/1‑69 (MQ=255)                                                                                                  aTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                             <  2:1051245/93‑1 (MQ=255)                                                                                                  aTGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                             >  1:1051245/1‑93 (MQ=255)                                                                                                        tatcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAAt          <  1:737544/136‑1 (MQ=255)                                                                                                          tcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg        <  2:127577/138‑1 (MQ=255)                                                                                                            tttCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCa      <  1:533383/139‑1 (MQ=255)                                                                                                                aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                   >  1:297531/1‑90 (MQ=255)                                                                                                                aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcg                                                   <  2:297531/90‑1 (MQ=255)                                                                                                                aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct  <  2:977719/139‑1 (MQ=255)                                                                                                                aGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATtct  >  2:351160/1‑139 (MQ=255)                                                                                                                 gTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                >  2:244269/1‑92 (MQ=255)                                                                                                                 gTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                <  1:244269/92‑1 (MQ=255)                                                                                                                     cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGa                              <  1:106359/106‑1 (MQ=255)                                                                                                                     cccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGa                              >  2:106359/1‑106 (MQ=255)                                                                                                                      ccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg        <  2:212536/127‑1 (MQ=255)                                                                                                                      ccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg        >  1:212536/1‑127 (MQ=255)                                                                                                                                cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAg                                                                                         >  1:947880/1‑36 (MQ=255)                                                                                                                                cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAg                                                                                         <  2:947880/36‑1 (MQ=255)                                                                                                                                cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                             <  2:256530/64‑1 (MQ=255)                                                                                                                                cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAt                                                             >  1:256530/1‑64 (MQ=255)                                                                                                                                   gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc         <  1:168047/113‑1 (MQ=255)                                                                                                                                   gCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATc         >  2:168047/1‑113 (MQ=255)                                                                                                                                       |                                                                                                                    AGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCT  >  NC_000913/3014153‑3014402

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