Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,760,776	T→G	L158F (TTA→TTC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,776	0	T	G	100.0%	214.7 / NA	56	L158F (TTA→TTC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/26);  total (30/26)

CGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGC  >  NC_000913/3760632‑3760923                                                                                                                                                 |                                                                                                                                                    cgtgcgGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACAcc                                                                                                                                               <  1:159056/151‑1 (MQ=255) cgtccgGTTGGTGGTCCCTCCTTCGAGCAGGCAGTTGCCTGCCTGACGCATAACATCGGGAATACTCAACGCGCCGCCAATCTCCAACAGCTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGTAGTTCCCGCCAACATCCAGAacac                                                                                                                                                >  2:763921/1‑150 (MQ=255)  gtgcgGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACAcc                                                                                                                                               >  2:858901/1‑150 (MQ=255)       gTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACAccgccgcc                                                                                                                                         <  2:1820/151‑1 (MQ=255)          ggtggtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACcgccgccgc                                                                                                                                       >  2:895476/1‑150 (MQ=255)                cccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGtt                                                                                                                                 >  1:840814/1‑150 (MQ=255)                cccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTa                                                                                                                                <  1:890294/151‑1 (MQ=255)                 ccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTAt                                                                                                                               <  2:530916/151‑1 (MQ=255)                    aCTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGa                                                                                                                            <  1:9448/151‑1 (MQ=255)                       tCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACg                                                                                                                          >  2:479106/1‑150 (MQ=255)                        cGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGa                                                                                                                         <  1:749737/150‑1 (MQ=255)                             aGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCTCCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAg                                                                                                                    <  1:506515/150‑1 (MQ=255)                                     gcctgcctGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTcc                                                                                                           >  2:186090/1‑151 (MQ=255)                                      cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCg                                                                                                          <  1:612331/151‑1 (MQ=255)                                      cctgcctgACGCATAACATCGGGAATACGAAACGCGACGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCg                                                                                                          <  1:509393/151‑1 (MQ=255)                                       ctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGc                                                                                                         >  1:742652/1‑151 (MQ=255)                                          cctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGcccc                                                                                                      >  2:818409/1‑151 (MQ=255)                                            tgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCg                                                                                                    >  2:45763/1‑151 (MQ=255)                                             gACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGccccc                                                                                                     <  2:727483/149‑1 (MQ=255)                                                gCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAAt                                                                                                >  2:63167/1‑151 (MQ=255)                                                  aTAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATAc                                                                                              >  2:899887/1‑151 (MQ=255)                                                  aTAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATAc                                                                                              >  2:389845/1‑151 (MQ=255)                                                         cGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGcaca                                                                                        <  2:732560/150‑1 (MQ=255)                                                         cGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCAcag                                                                                       <  2:775603/151‑1 (MQ=255)                                                          gggAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCAcagc                                                                                      >  1:650972/1‑151 (MQ=255)                                                               tACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTg                                                                                 >  2:751463/1‑151 (MQ=255)                                                                aCGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTgc                                                                                >  1:616607/1‑151 (MQ=255)                                                                    aaCGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGcc                                                                             >  1:691240/1‑150 (MQ=255)                                                                        ccccgccAATCGCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAgcg                                                                         <  2:635883/148‑1 (MQ=255)                                                                         gccgccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGcc                                                                       >  1:625093/1‑151 (MQ=255)                                                                           cgccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCg                                                                     <  2:233418/151‑1 (MQ=255)                                                                           cgccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCg                                                                     >  1:239211/1‑151 (MQ=255)                                                                             ccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGa                                                                    <  2:430198/150‑1 (MQ=255)                                                                              cAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATc                                                                  <  1:573104/151‑1 (MQ=255)                                                                                  ctcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATcgcg                                                               >  1:592865/1‑150 (MQ=255)                                                                                     caccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGt                                                           >  2:286187/1‑151 (MQ=255)                                                                                     caccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGt                                                           <  2:804506/151‑1 (MQ=255)                                                                                       ccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTcc                                                         <  1:915150/151‑1 (MQ=255)                                                                                             cgccgcGAGATACCCCCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTc                                                    <  2:236336/150‑1 (MQ=255)                                                                                                cgcgAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGc                                                <  1:95837/151‑1 (MQ=255)                                                                                                    agaTACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgtcg                                            >  2:617544/1‑151 (MQ=255)                                                                                                      aTACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCgtcgtc                                          >  1:494187/1‑151 (MQ=255)                                                                                                       tACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                         >  2:898310/1‑151 (MQ=255)                                                                                                       tACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                         >  1:637299/1‑151 (MQ=255)                                                                                                          caccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                         >  2:666647/1‑148 (MQ=255)                                                                                                            ccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                         >  2:635965/1‑146 (MQ=255)                                                                                                            ccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGAt                                     <  1:451221/150‑1 (MQ=255)                                                                                                                  ctTTTCCGCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATAc                                <  1:45763/149‑1 (MQ=255)                                                                                                                         gCTGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTc                       >  2:795013/1‑151 (MQ=255)                                                                                                                           tGGCAGTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                     >  1:627346/1‑151 (MQ=255)                                                                                                                                gTGGCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                 <  2:831890/150‑1 (MQ=255)                                                                                                                                  ggCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGaa              <  2:415987/151‑1 (MQ=255)                                                                                                                                  ggCCGCCAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGaa              <  1:172207/151‑1 (MQ=255)                                                                                                                                      gccAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa          <  1:895496/151‑1 (MQ=255)                                                                                                                                        cAACATCAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCAt         >  2:824749/1‑150 (MQ=255)                                                                                                                                              cAAGAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc  >  2:34616/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGC  >  NC_000913/3760632‑3760923

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