Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,849,128	G→T	Q199K (CAG→AAG)	selD ←	selenophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,849,128	0	G	T	100.0%	107.5 / NA	37	Q199K (CAG→AAG)	selD	selenophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/17);  total (20/17)

TTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGT  >  W3110S.gb/1849070‑1849194                                                           |                                                                   ttGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGttt                                                        >  1:2035297/1‑71 (MQ=255) ttGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGttt                                                        >  1:4198754/1‑71 (MQ=255)   ggCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTc                                                       <  1:1724413/70‑1 (MQ=255)    gCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCa                                                      <  1:3810942/70‑1 (MQ=255)    gCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCa                                                      <  1:3805510/70‑1 (MQ=255)          gACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGa                                               >  1:3841813/1‑71 (MQ=255)                      aTGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGc                                   >  1:828449/1‑71 (MQ=255)                        gTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGcc                                  >  1:1562747/1‑70 (MQ=255)                        gTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCg                                 >  1:2797522/1‑71 (MQ=255)                         ttCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTc                                                       <  1:3889386/48‑1 (MQ=255)                         ttCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTc                                                       <  1:2940122/48‑1 (MQ=255)                         ttCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTc                                                       <  1:1562407/48‑1 (MQ=255)                          tCATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg                               >  1:3193506/1‑71 (MQ=255)                           cATCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTg                               >  1:1654599/1‑70 (MQ=255)                            aTCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTgg                              <  1:2698699/70‑1 (MQ=255)                            aTCCGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTgg                              <  1:4136165/70‑1 (MQ=255)                              ccGGCACATCACTTCCGTCGCCAGTCCCTTATGtt                                                               >  1:832215/1‑35 (MQ=255)                              ccGGCACATCACTTCCGTCGCCAGTCCCTTATGtt                                                               >  1:4050237/1‑35 (MQ=255)                               cGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg                          <  1:2508444/71‑1 (MQ=255)                               cGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg                          <  1:3624488/71‑1 (MQ=255)                               cGGCACATCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAg                          <  1:1571586/71‑1 (MQ=255)                                    caTCACTTCCGTCGCCAGTCCCTTATGTTATGGTTTc                                                       <  1:3783085/37‑1 (MQ=37)                                     aTCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc                    <  1:357085/71‑1 (MQ=255)                                     aTCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc                    <  1:1344903/71‑1 (MQ=255)                                      tCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc                    >  1:4146303/1‑70 (MQ=255)                                      tCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc                    >  1:975507/1‑70 (MQ=255)                                      tCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc                    >  1:3306436/1‑70 (MQ=255)                                      tCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGcc                    >  1:3088426/1‑70 (MQ=255)                                      tCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCg                   <  1:684562/71‑1 (MQ=255)                                      tCACTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTAAGc                                   >  1:265687/1‑55 (MQ=255)                                        aCTTCCGTCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGAt                 >  1:256885/1‑71 (MQ=255)                                                       cccTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGt                             >  1:4332148/1‑44 (MQ=255)                                                       cccTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGt  >  1:305323/1‑71 (MQ=255)                                                       cccTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGT‑GTAAGAACGCCGATCCCCAGCGGTTTCGt  >  1:1149067/1‑70 (MQ=255)                                                        ccTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAg                                    >  1:3960630/1‑36 (MQ=255)                                                        ccTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTc    <  1:623847/68‑1 (MQ=255)                                                        ccTTATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTc    <  1:1769735/68‑1 (MQ=255)                                                           |                                                                   TTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTTTTCTCAGCCGTGGTAAGAACGCCGATCCCCAGCGGTTTCGT  >  W3110S.gb/1849070‑1849194

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