Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,847,303	G→A	H236Y (CAC→TAC)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,847,303	0	G	A	100.0%	182.0 / NA	54	H236Y (CAC→TAC)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/21);  total (33/21)

TAGCTGGCAAAGTTACGTTCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTT  >  NC_000913/1847167‑1847439                                                                                                                                         |                                                                                                                                         tAGCTGGCAAAGTTACGTTCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGcc                                                                                                                                       >  1:390517/1‑140 (MQ=255)    cTGGCAAAGTTACGTTCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTCCTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCcag                                                                                                                                    <  1:304630/140‑1 (MQ=255)    cTGGCAAAGTTACGTTCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCcag                                                                                                                                    <  1:1045881/140‑1 (MQ=255)    cTGGCAAAGTTACGTTCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCcag                                                                                                                                    <  2:390517/140‑1 (MQ=255)                  ttCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCg                                                                                                                      <  1:1426672/140‑1 (MQ=255)                            ccAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTTAAGTAGCCCAGCAGACCAAAGCCCg                                                                                                                      >  1:766956/1‑130 (MQ=255)                            ccAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTTAAGTAGCCCAGCAGACCAAAGCCCg                                                                                                                      <  2:766956/130‑1 (MQ=255)                            ccAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGt                                                                                                            >  2:1460721/1‑140 (MQ=255)                                     gcgcCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATc                                                                                                        >  2:105886/1‑135 (MQ=255)                                     gcgcCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATc                                                                                                        <  1:105886/135‑1 (MQ=255)                                      cgcCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCTCAGCAGATCAACGCGCGTAACATCGGTCATCGCtttt                                                                                                  >  1:395950/1‑140 (MQ=255)                                       gcCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTa                                                                                                 <  1:1193800/140‑1 (MQ=255)                                           aaCTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGGCCACAGCCCGTAGCGTCAGTCATCGCTTTTACGcc                                                                                             >  2:452689/1‑140 (MQ=255)                                            agttAATGTACTCTTCACCACCCGGGAGTTTCGGGTTCGTTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCt                                                                                            <  2:537439/138‑1 (MQ=255)                                                       tctTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCGCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGc                                                                                 >  2:1378132/1‑140 (MQ=255)                                                       tctTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGc                                                                                 >  2:904591/1‑140 (MQ=255)                                                       tctTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGc                                                                                 >  1:297184/1‑140 (MQ=255)                                                       tctTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGc                                                                                 >  2:121525/1‑140 (MQ=255)                                                        ctTCAACACCGGGGAGTTTCGGGATCACTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGCAACGTCGTTCATCGGTTTTACACCTTCGATGTTGGCa                                                                                >  1:549155/1‑140 (MQ=255)                                                         ttCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTTACGTCGGACATGGCATTTACGCCTTCGATGTTGGcct                                                                               >  1:537439/1‑138 (MQ=255)                                                              cacCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGAc                                                                          >  1:110439/1‑140 (MQ=255)                                                              cacCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGAc                                                                          >  1:166820/1‑140 (MQ=255)                                                               acCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGAcg                                                                         <  2:900129/140‑1 (MQ=255)                                                                 cGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGAcgcg                                                                       >  2:1542573/1‑140 (MQ=255)                                                                 cGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGAcgcg                                                                       >  1:861368/1‑140 (MQ=255)                                                                  ggggAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACgcgc                                                                      <  2:861368/140‑1 (MQ=255)                                                                  ggggAGTTTCGGGATCGATTCAGAGTCGACGCGTGACTGTACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACgcgc                                                                      <  2:742927/140‑1 (MQ=255)                                                                   gggAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGcc                                                                     >  1:1010564/1‑140 (MQ=255)                                                                   gggAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGcc                                                                     >  2:1527376/1‑140 (MQ=255)                                                                   gggAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGCAACGTCGGGCATGGATTTTACGCCTGCGATGCTGGCAAACGACGCGcc                                                                     >  2:694348/1‑140 (MQ=255)                                                                        tttCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGa                                                                >  2:1218290/1‑140 (MQ=255)                                                                         ttCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGAt                                                               <  1:1460721/140‑1 (MQ=255)                                                                         ttCGGGATCGCATCATAGTCGAACCGGGCCTCCAAACCAACCCCCTGAAACATTTCGCTAAAGTAGCCCAGCAAAACAAAGCCCGTAACGTCGGTCATCGCTTTTACACCTTCGATGTTGGCAAAGGACGCGCCTGCGAt                                                               <  2:549155/140‑1 (MQ=255)                                                                              gATCGCTTCATAGTCAACGCGTTCCGGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCa                                                          <  1:694348/140‑1 (MQ=255)                                                                                   cTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTGCGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCgg                                                     <  2:1010564/140‑1 (MQ=255)                                                                                   cTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCgg                                                     >  2:364303/1‑140 (MQ=255)                                                                                   cTTCATAGTCGACGCCTGCCTGCACCCCAGCCCCCTGCCACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCgg                                                     <  1:452689/140‑1 (MQ=255)                                                                                        tAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCgg                                                     <  1:923530/135‑1 (MQ=255)                                                                                        tAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCgg                                                     >  2:923530/1‑135 (MQ=255)                                                                                          gTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCa                                              <  1:121525/140‑1 (MQ=255)                                                                                              aCGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTAAAGTAGCCCAGCAGACCAAA‑TCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGAACATCACTTcc                                         <  2:1537015/140‑1 (MQ=255)                                                                                                        gCACGCCAGCCCCCTGAGACATATCGCTCAACTAGCCCACCAGACCAATGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGt                                <  2:395950/140‑1 (MQ=255)                                                                                                             ccAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCGCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTg                           <  1:1378132/140‑1 (MQ=255)                                                                                                              cAGCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGa                          >  2:225074/1‑140 (MQ=255)                                                                                                                gCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATg                        >  2:1470027/1‑140 (MQ=255)                                                                                                                gCCCCCTGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGATCATCCGTCACATCCCTTCCGACGCCAGTCCCTGATg                        >  1:603822/1‑140 (MQ=255)                                                                                                                      tGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGtt                      >  2:363950/1‑136 (MQ=255)                                                                                                                      tGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGtt                      <  1:363950/136‑1 (MQ=255)                                                                                                                      tGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTgg                  >  1:1391811/1‑140 (MQ=255)                                                                                                                      tGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTgg                  >  1:156612/1‑140 (MQ=255)                                                                                                                      tGACACATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTgg                  >  1:1234348/1‑140 (MQ=255)                                                                                                                        acacATTTCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGtt                >  2:941453/1‑140 (MQ=255)                                                                                                                              ttCGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACa          >  1:63585/1‑140 (MQ=255)                                                                                                                                cGCTCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGt        >  1:1337371/1‑140 (MQ=255)                                                                                                                                   tCAAGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGAt     <  1:1290984/140‑1 (MQ=255)                                                                                                                                      aGTAGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGCCGCCAGTCCCTTATGTTCTGGTTTCAACAGTGAtttt  >  1:579219/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                         TAGCTGGCAAAGTTACGTTCAGTGCCGCCAGGTACTGCGCCCAACTTAATGTACTCTTCAACACCGGGGAGTTTCGGGATCGCTTCATAGTCGACGCGTGCCTGCACACCAGCCCCCTGACACATTTCGCTCAAGTGGCCCAGCAGACCAAAGCCCGTAACGTCGGTCATCGCTTTTACGCCTTCGATGTTGGCAAAGGACGCGCCTGCGATGTTCATCCGGCACATCACTTCCGTCGCCAGTCCCTGATGTTCTGGTTTCAACAGTGATTTT  >  NC_000913/1847167‑1847439

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