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%	177.3 / NA	51	V166A (GTG→GCG)	xerD	site‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/27);  total (24/27)

CGCACGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCACGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGAC  >  NC_000913/3039099‑3039392                                                                                                                                                     |                                                                                                                                                  cgcACGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACcgcg                                                                                                                                                 >  2:723603/1‑151 (MQ=255)     cGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCt                                                                                                                                             >  1:484784/1‑151 (MQ=255)         ggctggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTg                                                                                                                                            <  2:217287/148‑1 (MQ=255)         ggctggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACg                                                                                                                                         >  1:18738/1‑151 (MQ=255)         ggcgggAAAACAACCCGACAATTGACACACCATTCAGCAGCCACGGACGCCCAGGTTCCAGATAGTTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCt                                                                                                                                             <  2:682054/147‑1 (MQ=255)          gctggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGc                                                                                                                                        >  1:798922/1‑151 (MQ=255)            tggGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGc                                                                                                                                        >  1:589961/1‑149 (MQ=255)                aaaCAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCt                                                                                                                                   <  1:819524/150‑1 (MQ=255)                    aacaCGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGata                                                                                                                               >  2:651684/1‑150 (MQ=255)                            aaTTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTc                                                                                                                         <  1:400333/148‑1 (MQ=255)                                gACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTc                                                                                                                   <  1:72678/150‑1 (MQ=255)                                 acacacCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTc                                                                                                                   >  2:194673/1‑149 (MQ=255)                                      cctttCAGCAGCCACGGACGCCCATGTTCCAGATAGTTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCg                                                                                                             <  2:413579/147‑1 (MQ=255)                                        aTTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGAcc                                                                                                          >  2:109726/1‑151 (MQ=255)                                         tccagcagCAACGCACGCCCATGTTCAAAATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGAcc                                                                                                          <  2:244329/148‑1 (MQ=255)                                                gcCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTc                                                                                                     >  1:487131/1‑148 (MQ=255)                                                 ccACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCaga                                                                                                  <  1:830300/150‑1 (MQ=255)                                                     ggACGCCCATGTTCCAGATAGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGac                                                                                               <  1:417710/149‑1 (MQ=255)                                                                ttCCAGATAGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCgg                                                                                    <  2:643262/149‑1 (MQ=255)                                                                 tCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGt                                                                                   >  2:719938/1‑149 (MQ=255)                                                                 tCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGt                                                                                   >  2:829262/1‑149 (MQ=255)                                                                    aaaTAGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAg                                                                                 <  2:630136/146‑1 (MQ=255)                                                                      atgggTTTCCAGCCAGTAATCCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTCCCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACTCGCAGCCCGGTAGCATa                                                                             <  2:140617/147‑1 (MQ=255)                                                                      aaaGGTTCCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATa                                                                             <  2:503102/148‑1 (MQ=255)                                                                       tAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATac                                                                            <  1:696329/150‑1 (MQ=255)                                                                        aGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATaca                                                                           >  2:308852/1‑150 (MQ=255)                                                                         ggTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATacaac                                                                         >  1:635653/1‑151 (MQ=255)                                                                           tttCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATacaaca                                                                        <  1:264411/150‑1 (MQ=255)                                                                               cagccagTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTc                                                                    >  1:460373/1‑150 (MQ=255)                                                                                  caagTAAGCCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAg                                                                 <  2:44011/148‑1 (MQ=255)                                                                                       aaaCCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGc                                                           >  2:808519/1‑151 (MQ=255)                                                                                        aaCCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCt                                                          <  2:61440/151‑1 (MQ=255)                                                                                          ccGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCttt                                                        <  1:654397/151‑1 (MQ=255)                                                                                                   cACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGt                                                 >  2:664489/1‑149 (MQ=255)                                                                                                          cGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGt                                        <  1:280914/151‑1 (MQ=255)                                                                                                             cACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGc                                     <  1:646928/151‑1 (MQ=255)                                                                                                               ccAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCt                                    <  2:689642/150‑1 (MQ=255)                                                                                                                cAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTAAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTg                                   >  2:369783/1‑150 (MQ=255)                                                                                                                        ctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGAtt                            <  2:658469/149‑1 (MQ=255)                                                                                                                        ctTTGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTa                           >  2:788258/1‑150 (MQ=255)                                                                                                                         tttGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAAt                         >  1:309007/1‑151 (MQ=255)                                                                                                                         tttGTTGCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAAt                         >  1:798775/1‑151 (MQ=255)                                                                                                                            gttgCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATgg                       <  1:719987/150‑1 (MQ=255)                                                                                                                               gCCTTTACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGcc                   >  1:563924/1‑151 (MQ=255)                                                                                                                                   ttACCAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGt                >  2:549766/1‑150 (MQ=255)                                                                                                                                      ccAATGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGtaata            <  1:199495/151‑1 (MQ=255)                                                                                                                                        aaTGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGtaataa           >  1:72831/1‑150 (MQ=255)                                                                                                                                        aaTGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATaaa          <  1:666563/151‑1 (MQ=255)                                                                                                                                         aTGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAAc         <  1:871007/151‑1 (MQ=255)                                                                                                                                         aTGACCCGTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAAc         <  1:164153/151‑1 (MQ=255)                                                                                                                                                gTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGAc  >  2:570290/1‑151 (MQ=255)                                                                                                                                                gTACCGCGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGAc  >  1:561080/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGCACGCTGGCTGGGAAACAACACGTCAATTGACACACCATTCAGCAGCCACGGACGCCCATGTTCCAGATAGGTTTCCAGCCAGTAAACCGCCTCTTCACCTAACGGCACCAGACGCTCTTTGTTGCCTTTACCAATGACCCGTACCACGCCCTGACGCAGGCTGATATCACTCATTGTCAGTCCGACCAGTTCAGAGACACGCAGCCCGGTAGCATACAACACTTCAAGCATGGCTTTATCGCGTAGCTCCAGTGGCTGATCGATTAATGGTGCCTGTAATAAACGTTCGAC  >  NC_000913/3039099‑3039392

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