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

AGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGA  >  NC_000913/3014153‑3014410                                                                                                                                       |                                                                                                                            aGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAAtt                                                                                                                        >  2:530348/1‑139 (MQ=255)     gAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGATAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTcaca                                                                                                                    >  1:663032/1‑139 (MQ=255)          tCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCt                                                                                                               >  2:183626/1‑139 (MQ=255)           cGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTg                                                                                                              <  1:306253/139‑1 (MQ=255)            ggCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTgc                                                                                                             >  2:614300/1‑139 (MQ=255)                       aTTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGTGCGCGATTGCa                                                                                                  <  2:644512/139‑1 (MQ=255)                         ttGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAgg                                                                                                <  2:722666/139‑1 (MQ=255)                          tGCGCCATAGCAAACGGATAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                               <  2:663032/139‑1 (MQ=255)                          tGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGa                                                                                               <  2:35690/139‑1 (MQ=255)                            cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                           >  2:762727/1‑109 (MQ=255)                            cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTa                                                                                                                           <  1:762727/109‑1 (MQ=255)                            cgcCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                             <  2:709274/139‑1 (MQ=255)                            cgcCATAGCAAAAGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGACAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                             <  1:488600/139‑1 (MQ=255)                                               ttttGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                                                                          <  1:183626/139‑1 (MQ=255)                                                               ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgt                                                          <  1:162264/139‑1 (MQ=255)                                                               ttGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgt                                                          <  2:337457/139‑1 (MQ=255)                                                                 gCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTAcgtcg                                                        <  1:188653/139‑1 (MQ=255)                                                                         tgtgAAAATATTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGttgttg                                                <  1:256052/139‑1 (MQ=255)                                                                                 tatTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                             <  1:743078/86‑1 (MQ=255)                                                                                 tatTCATACCAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                             >  2:743078/1‑86 (MQ=255)                                                                                         ccAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                <  2:497185/139‑1 (MQ=255)                                                                                         ccAACCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGc                                >  1:533957/1‑139 (MQ=255)                                                                                            aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                             <  2:74479/75‑1 (MQ=255)                                                                                            aCCTCAATGGATAT‑CTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATg                                                                                             >  1:74479/1‑75 (MQ=255)                                                                                                         atcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTGTGTGCAATc                 >  2:454960/3‑139 (MQ=255)                                                                                                          tcTTTCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                <  1:405700/138‑1 (MQ=255)                                                                                                              tCAGTAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCAtt            >  2:51991/1‑139 (MQ=255)                                                                                                                  tAACCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTca        >  2:386708/1‑139 (MQ=255)                                                                                                                    aCCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                                    >  2:438888/1‑45 (MQ=255)                                                                                                                    aCCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTg                                                                                                    <  1:438888/45‑1 (MQ=255)                                                                                                                    aCCCGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTcaca      >  2:557088/1‑139 (MQ=255)                                                                                                                      ccGGAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACAtt    >  2:224691/1‑139 (MQ=255)                                                                                                                        ggAATACCCGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTga  >  1:425199/1‑139 (MQ=255)                                                                                                                                cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                <  1:231957/117‑1 (MQ=255)                                                                                                                                cGGGCCGTTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCg                >  2:231957/1‑117 (MQ=255)                                                                                                                                       |                                                                                                                            AGTCGAGCCTCGGCTATCAACAATTTGCGCCATAGCAAACGGACAATTTTGCTCTTCGAGTTTTGCAGCCTCTGTGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTTGCAGCGGTTGAAGGAGCTCGCTTTGTGCAATCGCATTCTCACATTGA  >  NC_000913/3014153‑3014410

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