Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,964,837	G→A	R107C (CGT→TGT)	thyA ←	thymidylate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,964,837	0	G	A	100.0%	189.5 / NA	56	R107C (CGT→TGT)	thyA	thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/33);  total (23/33)

GAAGAATGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCA  >  NC_000913/2964706‑2964968                                                                                                                                    |                                                                                                                                    gaagaaTGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTg                                                                                                                             >  2:818848/1‑140 (MQ=255) gaagaaTGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTg                                                                                                                             >  1:778223/1‑140 (MQ=255) gaagaaTGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTg                                                                                                                             >  2:535924/1‑140 (MQ=255)  aagaaTGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTgg                                                                                                                            >  1:413857/1‑140 (MQ=255)   agaaTGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGc                                                                                                                           >  1:159608/1‑140 (MQ=255)      aTGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGtt                                                                                                                        >  1:419894/1‑140 (MQ=255)           tGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTGCAGTACCGTAGGGATCCGGTCAATACGACAACCATCTGGCGTTGGCCa                                                                                                                   >  2:305596/1‑140 (MQ=255)                  gTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGgcgcgc                                                                                                            <  1:747260/140‑1 (MQ=255)                   tGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGcc                                                                                                           >  1:658167/1‑140 (MQ=255)                   tGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGcc                                                                                                           >  1:197945/1‑140 (MQ=255)                   tGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGcc                                                                                                           <  1:45065/140‑1 (MQ=255)                   tGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGcc                                                                                                           <  2:38551/140‑1 (MQ=255)                   tGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGcc                                                                                                           >  1:585020/1‑140 (MQ=255)                   tGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATCACACCCATCTGGCGTTGGCCAGGCGCGcc                                                                                                           >  1:286232/1‑140 (MQ=255)                    gCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCa                                                                                                          <  1:414801/140‑1 (MQ=255)                      cAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACt                                                                                                        <  2:158700/140‑1 (MQ=255)                      cAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACt                                                                                                        <  1:655101/140‑1 (MQ=255)                      cAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTCTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACt                                                                                                        <  1:8790/140‑1 (MQ=255)                        gcgcCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGt                                                                                                      <  1:818848/140‑1 (MQ=255)                        gcgcCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGt                                                                                                      <  2:735096/140‑1 (MQ=255)                        gcgcCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGt                                                                                                      <  2:413857/140‑1 (MQ=255)                        gcgcCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGt                                                                                                      <  2:419894/140‑1 (MQ=255)                        gcgcCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGt                                                                                                      <  2:778223/140‑1 (MQ=255)                         cgcCATTTTATCCAGTTCGCCTCCGTTCCTCGCTGAATCAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGtt                                                                                                     <  1:371949/140‑1 (MQ=255)                         cgcCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGtt                                                                                                     <  1:719171/140‑1 (MQ=255)                           ccATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGCGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGTGCGCCACTGTTTa                                                                                                   <  2:475945/140‑1 (MQ=255)                             aTTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAGCAAACTGGCGTTGGCCAAGCGCGCCACTGTTTAcc                                                                                                 >  2:371949/1‑140 (MQ=255)                                                     cACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGtt                                                                         >  1:154131/1‑140 (MQ=255)                                                            aaaCAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCATGCGCGCCACTGTTTACAATACACTGTCCCGGGGTCGCCGTTTTCATCg                                                                  >  2:332779/1‑140 (MQ=255)                                                            aaaCAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCg                                                                  >  2:6873/1‑140 (MQ=255)                                                                 ataatGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCa                                                             <  2:34018/140‑1 (MQ=255)                                                                      gcggcgCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGt                                                        <  2:154131/140‑1 (MQ=255)                                                                      gcggcgCGAATCCGGGTCGTTTTTCAGCTGGTTCAGCACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGt                                                        >  2:730394/1‑140 (MQ=255)                                                                        ggcgCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTcc                                                      <  2:585020/140‑1 (MQ=255)                                                                        ggcgCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTcc                                                      <  1:6873/140‑1 (MQ=255)                                                                        ggcgCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTcc                                                      <  2:197945/140‑1 (MQ=255)                                                                        ggcgCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCACCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTcc                                                      <  1:753336/140‑1 (MQ=255)                                                                           gcgAATCCGGGTCGTTTTTCGGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTcc                                                      <  1:100517/137‑1 (MQ=255)                                                                           gcgAATCCGGGTCGTTTTTCGGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTcc                                                      >  2:100517/1‑137 (MQ=255)                                                                            cgAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGa                                                  <  1:652994/140‑1 (MQ=255)                                                                            cgAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGACGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGa                                                  <  2:159608/140‑1 (MQ=255)                                                                                    ggTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACa                                          >  2:809851/1‑140 (MQ=255)                                                                                          ttttCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGttt                                    <  2:658167/140‑1 (MQ=255)                                                                                               aGCTGGTTCAGCACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCgtg                               <  1:730394/140‑1 (MQ=255)                                                                                                   ggTTCAGTACCGTAGCGATCTGATCAATATGACAACCCTCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGa                           <  1:305596/140‑1 (MQ=255)                                                                                                    gTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTa                             <  1:64798/137‑1 (MQ=255)                                                                                                    gTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTa                             >  2:64798/1‑137 (MQ=255)                                                                                                    gTTCAGTACCGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGAt                          <  1:19154/140‑1 (MQ=255)                                                                                                    gTTCAGTACCGTAGTGAACTGGTCAATATGACAACCATCTGGCGTTGGGCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGAt                          <  2:286232/140‑1 (MQ=255)                                                                                                             cGTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGt                 <  1:194263/140‑1 (MQ=255)                                                                                                              gTAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGtt                >  1:458882/1‑140 (MQ=255)                                                                                                               tAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCGGATGGTGCCATTGTGTTCGTGGAGATAAGCAATGTTa               >  1:333015/1‑140 (MQ=255)                                                                                                               tAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGa                           >  1:88877/1‑128 (MQ=255)                                                                                                               tAGTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGa                           <  2:88877/128‑1 (MQ=255)                                                                                                                 gTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTATATAAGc                      <  2:23426/131‑1 (MQ=255)                                                                                                                 gTGATCTGGTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGc                      >  1:23426/1‑131 (MQ=255)                                                                                                                         gTCAATATGACAACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCt     >  1:658812/1‑140 (MQ=255)                                                                                                                            aaTATGACAACCATCTGGCGTTAGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCa  <  1:332779/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    GAAGAATGCATGGCACGGTGCCAGCGCCATTTTATCCAGTTCGCCTACGTTCCACGCTGAAACAATAATGCGGCGCGAATCCGGGTCGTTTTTCAGCTGGTTCAGTACCGTAGTGATCTGGTCAATATGACGACCATCTGGCGTTGGCCAGGCGCGCCACTGTTTACCATACACTGGCCCGAGGTCGCCGTTTTCATCGGCCCATTCGTCCCAGATGGTGACATTGTTTTCGTGTAGATAAGCAATGTTAGTGTCGCCCTGCA  >  NC_000913/2964706‑2964968

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