Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,466	G→A	A154T (GCG→ACG)	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,466	0	G	A	93.5%	86.4 / NA	32	A154T (GCG→ACG)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (15/14);  minor base C (1/1);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAAT  >  NC_000913/4636219‑4636697                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        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)                                                                                                                                                                                                    atttggatttgaacaaaCCCGCGGCTCATATTAGTTGTTGTGACACGCCATTACCGCTTACGCGTTTTGATTTTTTATTGCTGAACACGTTACTCAGGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCCTTGATCCCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATTCTCAt                                    <  2:69122/236‑1 (MQ=255)                                                                                                                                                                                                            ttgaatGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCAt                            >  1:44931/1‑251 (MQ=255)                                                                                                                                                                                                                              cAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCATCTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAgggg          >  2:70297/1‑251 (MQ=255)                                                                                                                                                                                                                               aGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGc         >  2:246974/1‑251 (MQ=255)                                                                                                                                                                                                                                       tGGTTTGACACGCCATAAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGCTGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGAGGCCTTTAAt  >  2:65066/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        GCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCC‑TGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAAT  >  NC_000913/4636219‑4636697

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