Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,385	A→T	T127S (ACG→TCG)	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,385	0	A	T	100.0%	89.4 / NA	34	T127S (ACG→TCG)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/18);  total (16/18)

CCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGC  >  NC_000913/4636148‑4636617                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         cccGACGTCATGATTCTCGATGTTGGTTTTCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc                                                                                                                                                                                                                               <  1:97479/249‑1 (MQ=255) cccGACGTCATGATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTccc                                                                                                                                                                                                                               <  2:206207/249‑1 (MQ=255)            gATTCTCGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTg                                                                                                                                                                                                                  <  1:81760/251‑1 (MQ=255)                  cGATGTTGGTTTGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCAtt                                                                                                                                                                                                            >  2:94694/1‑251 (MQ=255)                                         gCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGCCGGCCCGAAGAGAAGCGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCGAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGATCTCGTCGCCGTCTCCCGTCATCCGTATTGGCCATTTTGAATTGAATGAACCCgcggcg                                                                                                                                                                                     >  2:204436/1‑251 (MQ=255)                                          cGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCAACCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAGGGAGCGCACCTTACTGCGTCGGGGGAAGAAGTTCTCGTCGCCGTCTACCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGCGgcgc                                                                                                                                                                                    >  2:63392/1‑251 (MQ=255)                                                           gccAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGCCGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGac                                                                                                                                                                   >  1:141856/1‑251 (MQ=255)                                                                        gccccATCCGGCGTTCCCTGTACTGTTCCGGACGGCCCGAAGTGAAGAGGTGGA‑CGCCCTGTTTGGG‑TGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCt                                                                                                                                                     <  1:13858/248‑1 (MQ=255)                                                                          tCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTg                                                                                                                                                    >  2:175470/1‑251 (MQ=255)                                                                                       aCCTGTTCTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAg                                                                                                                                        <  1:175470/250‑1 (MQ=255)                                                                                           ttACTTTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCCCTTACTGCGTCGGGGGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGtttt                                                                                                                                    <  1:157123/249‑1 (MQ=255)                                                                                           gTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCCCTTACTGCGTCGGTGGAAGAAGTTCTCGTCGCCGTCTCCCGTCACCCGTATCGGTCATTTTGAATTGAATGAACCCGCGGCGCAGTTCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGttt                                                                                                                                     <  1:104228/249‑1 (MQ=255)                                                                                                       ccGGCCCGAAGGGAAGAGGTCGAGCGCC‑TGCTTGGGCTGGAACTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCTGGGTGCGCACCTTACTGCTGCGGGTGAAGAAGTTCTCGTCGCCGTCTCCTGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGTTGGTTTGACACGCCATTAACGCGGACGCGGTATGAGTTTTTATTGCTGAAGAc                                                                                                                       <  2:106416/250‑1 (MQ=255)                                                                                                                  tGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGt                                                                                                            >  1:108134/1‑251 (MQ=255)                                                                                                                   gAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTc                                                                                                           <  1:23199/251‑1 (MQ=255)                                                                                                                   gAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGGGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTc                                                                                                           <  1:13100/251‑1 (MQ=255)                                                                                                                        gggcGATCGCC‑GGCTTGGGCGGGAAATTGGTGCGGACGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACGGCGTCGGGTGAAGAAGTCCTCGTCGCCGTCTCCGGTCATCCGTAGTGGTCATTTTGAATTGAATGAACCCGGGGGGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGGTACTCAAGtggccag                                                                                                      <  2:141856/248‑7 (MQ=255)                                                                                                                               cGCC‑TGCTTGGGCTGGACATTGGTGCTGACGCCTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCAAGGCCGCGt                                                                                               >  2:94127/1‑251 (MQ=255)                                                                                                                                         gggCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAGGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCg                                                                                      >  2:84893/1‑251 (MQ=255)                                                                                                                                           gCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGcc                                                                                    <  1:94694/251‑1 (MQ=255)                                                                                                                                               gAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGCAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGc                                                                                 >  2:219079/1‑250 (MQ=255)                                                                                                                                                           gacgacTATGTGGCTAAACCGTTTTCCCCCGGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAAAAGTTCTCGTCGCCGTCTCCCGTCACCCGTATTGGTCATTTTGAATTGAATGAACCCGGGGCGCAGAGCAGCTGTTTTGACACGCCATTAGCGCTGACGCGGTATGTGTTTTTATTGCTGAAGTCGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGAt                                                                      <  1:113839/249‑1 (MQ=255)                                                                                                                                                                  atgTGGCTAAACCGTTTTCACCCCGCGAAGTGGGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTACCGCTGACGCGGGATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTgg                                                             >  1:216935/3‑251 (MQ=255)                                                                                                                                                                     tGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGa                                                           >  1:77076/1‑250 (MQ=255)                                                                                                                                                                                    cACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTGCTCAAGTCGCCAGGCCGCGTCTGGGCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACAcc                                           <  1:160178/251‑1 (MQ=255)                                                                                                                                                                                            gaggtgggcgcCGGGGTCCGCCCCTTACTGCTGCGGGTGAAGAATTTCTCGTCGCCGTCCCCCGTCACCCTTATTGGTCATTTTGAATTGAATGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCGGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCg                                   <  2:207501/244‑1 (MQ=255)                                                                                                                                                                                              aGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGc                                  >  2:234355/1‑250 (MQ=255)                                                                                                                                                                                               gtgtGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAc                                >  2:126389/1‑251 (MQ=255)                                                                                                                                                                                                ggtgCGCCGGGGTCCGCACTTTACGGCGTGGGGGGAAGAAGTTCGCGTCGCCGTCCCCCGTCTTCCGTTTTGGTCATTTGGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACTTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                               <  1:179217/250‑1 (MQ=255)                                                                                                                                                                                                        aGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAc                       >  1:264928/1‑251 (MQ=255)                                                                                                                                                                                                             gctcACCTTACTGCGTCGGGGGAAGAAGTTCGCGTCGCCGTCCCCCGACACCCTTATTGGTCATTTGGAATTGAATGAACCCGGGGCGCAGATCAGCTGGTTTGCCACGCCATTAACGCTGACGCGGTATGGGGTTTTATTGCTGAAGACGTTCCTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCcaca                  <  1:84893/248‑1 (MQ=255)                                                                                                                                                                                                             gcgcACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGATCGCACCGTTGATAccccc                   >  2:215005/1‑248 (MQ=255)                                                                                                                                                                                                                     tACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACAAGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTCAAACg          >  1:114376/1‑251 (MQ=255)                                                                                                                                                                                                                             ggggTGAAGAAGTTCTCGTCGCCGCCCCCCGGAATCCGTATGGGTCATTTGGAATTGAATAACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGc  <  1:147149/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         CCCGACGTCATGATTCTCGATGTTGGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGC  >  NC_000913/4636148‑4636617

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