Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,259,283	G→A	H123Y (CAT→TAT)	tdcG ←	L‑serine deaminase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,259,283	0	G	A	100.0%	121.0 / NA	35	H123Y (CAT→TAT)	tdcG	L‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

GCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGC  >  NC_000913/3259047‑3259529                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       gctgcGCCGCCCTAGCTCTTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACAttt                                                                                                                                                                                                                                          <  2:159164/251‑1 (MQ=255) gctgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACAtt                                                                                                                                                                                                                                           <  1:116060/250‑1 (MQ=255)   tgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCg                                                                                                                                                                                                                                        >  2:65736/1‑251 (MQ=255)             tAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTAAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGACGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGGAATAGGTTTTACTTAATCGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAAAGAGGCAACATTTCGGGATGGa                                                                                                                                                                                                                                 >  1:278367/1‑248 (MQ=255)                     tgtgCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGtt                                                                                                                                                                                                                       <  2:67969/250‑1 (MQ=255)                                         ataGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATcaa                                                                                                                                                                                                   <  1:253353/250‑1 (MQ=255)                                                                       agcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAgccg                                                                                                                                                                    >  1:247085/1‑251 (MQ=255)                                                                       agcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAgccg                                                                                                                                                                    >  2:291012/1‑251 (MQ=255)                                                                          agTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGccgccc                                                                                                                                                                 <  2:145462/251‑1 (MQ=255)                                                                               aCCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTgc                                                                                                                                                            <  1:87014/251‑1 (MQ=255)                                                                                ccTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCCGCCGCCCGCTgcc                                                                                                                                                           >  2:250424/1‑250 (MQ=255)                                                                                 ctgctgAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTgcgc                                                                                                                                                          >  1:257601/1‑251 (MQ=255)                                                                                 ctgctgAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTgcgc                                                                                                                                                          >  1:129856/1‑251 (MQ=255)                                                                                           gAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTGCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGAGCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAAGCAACAATATGCGCACCAGGTGCCACTGGCAGCCGCCCGCTGCGCGTTACTAAct                                                                                                                                                <  1:280754/251‑1 (MQ=255)                                                                                                   aaGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATaaa                                                                                                                                         <  1:110351/250‑1 (MQ=255)                                                                                                           gACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGAt                                                                                                                                >  2:23151/1‑251 (MQ=255)                                                                                                             cGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGAtct                                                                                                                              >  2:50328/1‑251 (MQ=255)                                                                                                                  cGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCtcatca                                                                                                                         <  1:101507/251‑1 (MQ=255)                                                                                                                  cGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCtcatc                                                                                                                          >  2:97202/1‑250 (MQ=255)                                                                                                                     cATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGACAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATg                                                                                                                      <  2:212370/251‑1 (MQ=255)                                                                                                                      aTCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGa                                                                                                                     <  2:257601/251‑1 (MQ=255)                                                                                                                            cGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACAt                                                                                                               <  1:216454/251‑1 (MQ=255)                                                                                                                                  gCCGAAGTGTTCTTCTTCGACAATAACCCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCg                                                                                                         <  1:168709/251‑1 (MQ=255)                                                                                                                                                         tAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCa                                                                                                                                                                        >  1:132650/1‑165 (MQ=255)                                                                                                                                                         tAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACAGTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCa                                                                                                                                                                        <  2:132650/165‑1 (MQ=255)                                                                                                                                                             cccGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCatgatg                                                                               >  1:51748/1‑250 (MQ=255)                                                                                                                                                                gccgccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATgg                                                                           <  2:57260/251‑1 (MQ=255)                                                                                                                                                                                             aTAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGACTTACGTTCTCATAGTGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGGTACCGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCAGACCGCTGCGCGTTACTAACTATATAACTGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAATCCCATGATGATGGCGACATCCGAGGCATGGCCTTTGCCCGTc                                              >  1:281533/1‑251 (MQ=255)                                                                                                                                                                                                   ctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTCCTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTg                                          <  1:70108/249‑1 (MQ=255)                                                                                                                                                                                                    ttCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGAc                                        >  2:286465/1‑250 (MQ=255)                                                                                                                                                                                                            ccTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGccgccc                                                                                                                                                                 >  1:93061/1‑121 (MQ=255)                                                                                                                                                                                                            ccTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGccgccc                                                                                                                                                                 <  2:93061/121‑1 (MQ=255)                                                                                                                                                                                                                          aTCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGCTCTTTGCTACAGGAAAATCAACAATATGCGCACAAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACcac                 >  2:224687/1‑251 (MQ=255)                                                                                                                                                                                                                            ccGCATTCCGTTCTCATAGCGAGCCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACACTATGCGCACCAGATGCCACTGGCAGCCGCACGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACcaca                >  2:200595/1‑250 (MQ=255)                                                                                                                                                                                                                                         tcATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCCTcgc  >  2:79804/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       GCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCGC  >  NC_000913/3259047‑3259529

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