Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,996,724	T→C	G146G (GGT→GGC)	xerC →	site‑specific tyrosine recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,996,724	0	T	C	100.0%	193.3 / NA	55	G146G (GGT→GGC)	xerC	site‑specific tyrosine recombinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/25);  total (30/25)

TTTCGGCACCGAAAGCGCCGCGTCATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGTCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGCAC  >  NC_000913/3996588‑3996859                                                                                                                                         |                                                                                                                                        tttCGGCACCGAAAGCGCCGCGTCATCTGCCGAAAAACATCGACGTCCAGGATATGAATCGGCTGCGGGATATTGATATCAATGATCCCCTCGCTGTACGCCACCGTGCCATGCTGGAAGTGATGTACGGCGCGGGCCtc                                                                                                                                      >  2:5439/1‑139 (MQ=255)   tCGGCACCGAAAGCGCCGCGTCATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCg                                                                                                                                    <  2:1000606/140‑1 (MQ=255)           gAAAGCGCCGCGTCATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTg                                                                                                                            <  1:1061741/140‑1 (MQ=255)             aaGCGCCGCGTCATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAg                                                                                                                          <  1:628285/140‑1 (MQ=255)                      gTCATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTgggg                                                                                                                 <  2:1003255/140‑1 (MQ=255)                        cATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCt                                                                                                               <  1:120710/140‑1 (MQ=255)                          tCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTgg                                                                                                             <  2:683839/140‑1 (MQ=255)                               cGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTgggg                                                                                                                 >  1:256193/1‑131 (MQ=255)                               cGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTgggg                                                                                                                 <  2:256193/131‑1 (MQ=255)                                      cATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGT‑ATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCt                                                                                                <  1:599236/140‑1 (MQ=255)                                              tcgacgATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTg                                                                                         <  1:685707/140‑1 (MQ=255)                                              tcgacgATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTg                                                                                         <  2:250079/140‑1 (MQ=255)                                                gacgaTATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGa                                                                                       <  1:293642/140‑1 (MQ=255)                                                               ctgctgGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTg                                                                               <  1:375248/133‑1 (MQ=255)                                                               ctgctgGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTg                                                                               >  2:375248/1‑133 (MQ=255)                                                               ctgctgGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGtgtg                                                                        >  2:253047/1‑140 (MQ=255)                                                               ctgctgGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGtgtg                                                                        >  1:307979/1‑140 (MQ=255)                                                                 gctgGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTggg                                                                      <  1:253047/140‑1 (MQ=255)                                                                          tGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGa                                                             >  1:653073/1‑140 (MQ=255)                                                                               tCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGgc                                                        >  2:766270/1‑140 (MQ=255)                                                                                    gATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCaa                                                   >  1:119316/1‑140 (MQ=255)                                                                                         ccTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATgggg                                                              >  1:497164/1‑124 (MQ=255)                                                                                         ccTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATgggg                                                              <  2:497164/124‑1 (MQ=255)                                                                                          cTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAgcg                                             >  2:206940/1‑140 (MQ=255)                                                                                           tCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAgcgc                                            <  2:686817/140‑1 (MQ=255)                                                                                           tCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAgcgc                                            >  2:440187/1‑140 (MQ=255)                                                                                           tCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAgcgc                                            >  1:531592/1‑140 (MQ=255)                                                                                              cTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGcgccgc                                         <  2:531592/140‑1 (MQ=255)                                                                                                 tACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTg                                      >  1:82058/1‑140 (MQ=255)                                                                                                   cgcgACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGcc                                    >  1:662220/1‑140 (MQ=255)                                                                                                         cGTGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGCGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTgg                              >  1:458953/1‑140 (MQ=255)                                                                                                           tGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTc                            >  1:11360/1‑140 (MQ=255)                                                                                                           tGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTc                            >  1:149777/1‑140 (MQ=255)                                                                                                           tGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTc                            >  1:1051866/1‑140 (MQ=255)                                                                                                           tGCAATGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGGGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTc                            >  1:1069312/1‑140 (MQ=255)                                                                                                                tGCTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAAc                       >  1:618820/1‑140 (MQ=255)                                                                                                                  cTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGc                     >  1:449200/1‑140 (MQ=255)                                                                                                                  cTGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGc                     >  2:224042/1‑140 (MQ=255)                                                                                                                   tGGAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCt                    <  2:618820/140‑1 (MQ=255)                                                                                                                     gAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACCCCTCGCCCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCtgt                  <  2:1069312/140‑1 (MQ=255)                                                                                                                     gAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCaa                        <  1:709398/134‑1 (MQ=255)                                                                                                                     gAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCaa                        >  2:709398/1‑134 (MQ=255)                                                                                                                     gAAGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCtgt                  >  1:12751/1‑140 (MQ=255)                                                                                                                      aaGTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCtgtg                 <  2:307979/140‑1 (MQ=255)                                                                                                                        gTGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGc               >  2:144090/1‑140 (MQ=255)                                                                                                                         tGATGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCg              <  1:440187/140‑1 (MQ=255)                                                                                                                           aTGTACGGCGCGGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTg            >  2:733631/1‑140 (MQ=255)                                                                                                                                   cgcgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCt                    <  2:688346/124‑1 (MQ=255)                                                                                                                                   cgcgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCt                    >  1:688346/1‑124 (MQ=255)                                                                                                                                   cgcgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGc    >  2:971821/1‑140 (MQ=255)                                                                                                                                   cgcgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGc    >  1:951618/1‑140 (MQ=255)                                                                                                                                    gcgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGCa   <  2:449200/140‑1 (MQ=255)                                                                                                                                    gcgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGCa   <  2:82058/140‑1 (MQ=255)                                                                                                                                     cgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCt                    <  2:300080/122‑1 (MQ=255)                                                                                                                                     cgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCt                    >  1:300080/1‑122 (MQ=255)                                                                                                                                     cgGGCCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGCAc  >  2:172259/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        TTTCGGCACCGAAAGCGCCGCGTCATCTGCCGAAAAACATCGACGTCGACGATATGAATCGGCTGCTGGATATTGATATCAATGATCCCCTCGCTGTACGCGACCGTGCAATGCTGGAAGTGATGTACGGCGCGGGTCTGCGTCTTTCTGAGCTGGTGGGGCTGGATATTAAACACCTCGACCTGGAGTCTGGTGAAGTGTGGGTTATGGGGAAAGGCAGCAAAGAGCGCCGCCTGCCGATTGGTCGCAACGCTGTGGCGTGGATTGAGCAC  >  NC_000913/3996588‑3996859

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