Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,201	G→A	L65L (TTG→TTA)	creB →	DNA‑binding transcriptional regulator CreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,201	0	G	A	100.0%	103.1 / NA	33	L65L (TTG→TTA)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/13);  total (20/13)

TAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGC  >  NC_000913/4635961‑4636439                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               tAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGGCTACC‑GACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAg                                                                                                                                                                                                                                       <  1:128804/250‑1 (MQ=255)   gcgcGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTa                                                                                                                                                                                                                                    <  1:70398/251‑1 (MQ=255)      cgGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTg                                                                                                                                                                                                                                 >  1:106416/1‑251 (MQ=255)            cTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTATGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTc                                                                                                                                                                                                                           >  1:25444/1‑251 (MQ=255)             tGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTcc                                                                                                                                                                                                                          >  1:49035/1‑251 (MQ=255)              gTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGCTGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGAGCTcc                                                                                                                                                                                                                          >  2:335/1‑250 (MQ=255)                   gCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCgg                                                                                                                                                                                                                    <  1:149482/251‑1 (MQ=255)                            gcAATAACAGAGGCGGTTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATAAGCTAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGAAGG‑CTCCCGACGTCATGATTCTCGATGTTGTTTTGCCTGATATCGGCGGCTTTGAATTATGCCGCAAGTTACTGGCGCTCCATCCGGCGTTACCTg                                                                                                                                                                                                           >  2:62404/1‑251 (MQ=255)                            gcAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGCGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAACGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTg                                                                                                                                                                                                           >  1:252618/1‑251 (MQ=255)                                                   aaCGGGAAAAGGTCTGGTTAGAGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAgg                                                                                                                                                                                    >  2:40168/1‑250 (MQ=255)                                                                cTGGTTAGTGTAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGCCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGCTGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTACTGACTGCCCGGAGTGAAGAGGTCGATc                                                                                                                                                                       >  2:179217/1‑251 (MQ=255)                                                                             gATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCt                                                                                                                                                          <  1:215874/251‑1 (MQ=255)                                                                                       ggATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGt                                                                                                                                                <  2:49035/251‑1 (MQ=255)                                                                                               gACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGCAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgacg                                                                                                                                         >  2:13100/1‑250 (MQ=255)                                                                                               gACACGCTGGACTACATGTTGCAGCAGGAGGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTCCCGGCGCTGGATAAAGCTCGGCAGCAGG‑CTCCCGACGACATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGGCGCCAGGTACTGGCGCTCCATCCGGCGTTCCCTGTACTGTTCCTGACGGCCCGGCGGGAAGAGGGCGATCGCCTGATTGGTCTGGAAACTGGCGCTgacg                                                                                                                                         >  2:143680/1‑250 (MQ=255)                                                                                                  acGCTGGTCTACATTTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTa                                                                                                                                     <  1:143816/251‑1 (MQ=255)                                                                                                   cGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTat                                                                                                                                    >  1:102045/1‑250 (MQ=255)                                                                                                      tGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTAtgtg                                                                                                                                 >  2:147149/1‑251 (MQ=255)                                                                                                                              gTTTTGCCGTCGAGGGCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCATG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATGATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTTTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGaa                                                                                                         <  1:258271/251‑1 (MQ=255)                                                                                                                                tttGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAgt                                                                                                       >  1:143676/1‑251 (MQ=255)                                                                                                                                         gAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGTGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAg                                                                                              <  1:212466/251‑1 (MQ=255)                                                                                                                                                         gCTTGCCGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTg                                                                              >  1:207501/1‑251 (MQ=255)                                                                                                                                                               gggTGCTGGATAAAGCGCGCCAGTAGG‑CTCCCGACGTCATGATTCTTGATTTTGGTTTGCATGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCTCTACATGCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGACCGCCTGCTTGGGCTGGAAAGTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGTTGCGCCCCTTACTGCGTCgg                                                                        <  2:98944/250‑1 (MQ=255)                                                                                                                                                               cGGTGCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCCGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGCAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTCAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGtagg                                                                        >  1:160967/1‑248 (MQ=255)                                                                                                                                                                   gCTGGATAAAGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGATTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGCCGACTATGTGGCTAAACCGTTTTCACCCCGCGTAGAGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTg                                                                     >  2:104228/1‑250 (MQ=255)                                                                                                                                                                        aTAAAGCTCGCCAGCAGG‑CTCCCGACGACATGATTCTCGATGTTGGCTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTgaaga                                                                 >  2:23199/1‑249 (MQ=255)                                                                                                                                                                            aGCTCGCCAGCAGG‑CTCCCGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTtct                                                           <  1:335/251‑1 (MQ=255)                                                                                                                                                                                             cccGACGTCATGATTCTCGATGTTGGTTTTCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc                                             <  1:97479/249‑1 (MQ=255)                                                                                                                                                                                             cccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc                                             <  2:206207/249‑1 (MQ=255)                                                                                                                                                                                                        gATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTg                                <  1:81760/251‑1 (MQ=255)                                                                                                                                                                                                              cGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCAtt                          >  2:94694/1‑251 (MQ=255)                                                                                                                                                                                                                                     gCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGCCGGCCCGAAGAGAAGCGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCGAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGATCTCGTCGCCGTCTCCCGTCATCCGTATTGGCCATTTTGAATTGAATGAACCCgcggcg   >  2:204436/1‑251 (MQ=255)                                                                                                                                                                                                                                      cGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCAACCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAGGGAGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCGTCGCCGTCTACCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGCGgcgc  >  2:63392/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               TAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATAAAGCTCGCAAGCAGG‑TACCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGC  >  NC_000913/4635961‑4636439

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