Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,039,247	A→G	V166A (GTG→GCG)	xerD ←	site‑specific tyrosine recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,039,247	0	A	G	100.0%	196.4 / NA	56	V166A (GTG→GCG)	xerD	site‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/37);  total (19/37)

GCACGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCACGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGACC  >  NC_000913/3039100‑3039393                                                                                                                                                    |                                                                                                                                                   gcACGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCgcgc                                                                                                                                                 >  1:146794/1‑151 (MQ=255)    cGCTGGCTGGGAAACAACACTTAAATTGACACCCCATTCAGCAGCCACGGACCCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCt                                                                                                                                              <  1:382677/151‑1 (MQ=255)    cGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCt                                                                                                                                              >  1:751506/1‑151 (MQ=255)     gctggctgGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTg                                                                                                                                             <  2:198512/151‑1 (MQ=255)        ggctggGAAACAACACGTCAATTGACACCCCATTCAGCAGCCACGGACGCCCATGTTCC‑‑ATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACg                                                                                                                                          <  1:121605/149‑1 (MQ=255)        ggctggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACg                                                                                                                                          >  1:110344/1‑151 (MQ=255)        ggctggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACg                                                                                                                                          <  2:672607/151‑1 (MQ=255)           tggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGc                                                                                                                                         >  1:355966/1‑149 (MQ=255)            gggAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAg                                                                                                                                       <  1:326631/150‑1 (MQ=255)               aaaCAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCt                                                                                                                                    <  2:213/150‑1 (MQ=255)               aaaCAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTg                                                                                                                                   >  1:733630/1‑151 (MQ=255)                 acaacaCGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGat                                                                                                                                 <  2:389312/151‑1 (MQ=255)                   aacaCGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGatat                                                                                                                               >  2:654240/1‑151 (MQ=255)                    acacGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATc                                                                                                                              <  2:245439/151‑1 (MQ=255)                           aaTTGACACCCAATTCAGCAGCCACGGACCCCCTTGTCCAAGAGAGGTTCCCACCCAGAAAACCGCCTCTTCTCCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCACCCTGACGCAGGCTGATATCACTCAt                                                                                                                        <  1:606961/150‑1 (MQ=255)                           aaTTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCCGACGCAGGCTGATATCACTCAtt                                                                                                                       >  2:193631/1‑151 (MQ=255)                             ttGACACACCATTCAGCACCCACGGACGCCCATGTCCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTg                                                                                                                      <  1:513456/150‑1 (MQ=255)                                 cacacCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGt                                                                                                                 <  2:647696/151‑1 (MQ=255)                                  acacCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTc                                                                                                                >  1:566732/1‑151 (MQ=255)                                     ccATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCg                                                                                                              <  1:713354/150‑1 (MQ=255)                                                ccACGGACGCCCATGTTCCAGATAGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCag                                                                                                    <  1:451857/149‑1 (MQ=255)                                                         cccATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGc                                                                                         <  2:205026/151‑1 (MQ=255)                                                          ccAGGTCCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGcc                                                                                        <  1:573059/151‑1 (MQ=255)                                                               ttCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCgg                                                                                     <  1:34915/149‑1 (MQ=255)                                                                 ccAGATAGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAg                                                                                  <  2:240227/150‑1 (MQ=255)                                                                   aGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCAt                                                                               >  2:281507/1‑151 (MQ=255)                                                                   aGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCAt                                                                               >  1:516129/1‑151 (MQ=255)                                                                      aaGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAACATaca                                                                            <  2:184989/150‑1 (MQ=255)                                                                       aGGTTTCCACCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATaca                                                                            <  1:465984/150‑1 (MQ=255)                                                                         gTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATacaac                                                                          >  2:603463/1‑150 (MQ=255)                                                                                gccagTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCaa                                                                   <  2:415082/150‑1 (MQ=255)                                                                                 ccagTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAg                                                                  <  2:110376/150‑1 (MQ=255)                                                                                   agTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCAt                                                               >  2:514031/1‑151 (MQ=255)                                                                                        aCCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCtt                                                          >  1:433869/1‑151 (MQ=255)                                                                                                  cACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTa                                                 <  1:643921/150‑1 (MQ=255)                                                                                                         cGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGt                                         >  1:476026/1‑151 (MQ=255)                                                                                                                aGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGAt                                  >  2:95754/1‑151 (MQ=255)                                                                                                                aGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGAt                                  <  2:14858/151‑1 (MQ=255)                                                                                                                 gACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATc                                 <  2:349363/151‑1 (MQ=255)                                                                                                                     ctctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                             >  1:77227/1‑151 (MQ=255)                                                                                                                     ctctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                             >  2:648069/1‑151 (MQ=255)                                                                                                                      tctTTGTTGCCTTTACCAATGACCCGTACCGCCCCCTGACGCAGGCTGATATCACTCAATGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                             >  1:409409/1‑150 (MQ=255)                                                                                                                       ctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACTCTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                             <  1:296974/149‑1 (MQ=255)                                                                                                                       ctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                             <  2:391012/149‑1 (MQ=255)                                                                                                                       ctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                             <  1:761641/149‑1 (MQ=255)                                                                                                                         ttgttgCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAAt                          <  2:668575/150‑1 (MQ=255)                                                                                                                           gttgCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATgg                        <  2:733704/150‑1 (MQ=255)                                                                                                                                cTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTg                  >  1:243597/1‑151 (MQ=255)                                                                                                                                  ttACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTg                  <  2:7936/149‑1 (MQ=255)                                                                                                                                   tACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGta                <  1:661713/150‑1 (MQ=255)                                                                                                                                   tACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGta                <  2:561285/150‑1 (MQ=255)                                                                                                                                    aCCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGtaat              >  1:258488/1‑151 (MQ=255)                                                                                                                                    aCCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGtaa               <  2:692913/150‑1 (MQ=255)                                                                                                                                     ccAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGtaata             <  2:15890/151‑1 (MQ=255)                                                                                                                                        aTGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAAc          <  2:209973/151‑1 (MQ=255)                                                                                                                                        aTGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAAc          <  2:360782/151‑1 (MQ=255)                                                                                                                                                 aCCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGAcc  <  2:401070/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCACGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCACGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGACC  >  NC_000913/3039100‑3039393

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