Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,610	C→T	L202L (CTG→TTG)	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,610	0	C	T	96.9%	94.8 / ‑4.7	34	L202L (CTG→TTG)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/13);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCG  >  NC_000913/4636367‑4636856                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        cGGGAGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGAATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCCGGGCACCTACCATCGCACCGTTGATACCCACATTAAAACGTTGCGTg                                                                                                                                                                                                                                                   >  2:52470/1‑250 (MQ=255)                                            gtcATTTGGAATTGAATGAACCCGCGGCGAATACAAGCTG‑TTTGAACACGCCATTAACGCGGACGCGGTATGAGTTTTTATTGCTAAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACt                                                                                                                                                                                                       <  1:181600/249‑1 (MQ=255)                                                          aatgaaCCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATggg                                                                                                                                                                                         >  1:11875/1‑251 (MQ=255)                                                              aaCCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGatat                                                                                                                                                                                     >  1:134914/1‑251 (MQ=255)                                                                      gcggaaaTCAGCTGGTTGGA‑CACGCCATTAACGCTGACGCGTTATGAGTTTTTATTGCTGAAGACGTTACTCAAGCCGCCAGGCCGCGTCTGGGCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAg                                                                                                                                                                             <  1:376864/245‑1 (MQ=255)                                                                        gcAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAggg                                                                                                                                                                           >  1:302516/1‑251 (MQ=255)                                                                        gcAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAggg                                                                                                                                                                           >  1:265283/1‑251 (MQ=255)                                                                         cAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAgggg                                                                                                                                                                          >  1:9772/1‑251 (MQ=255)                                                                               aGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                                                                                                                                                                     >  1:67123/1‑250 (MQ=255)                                                                                       tGA‑CATGCCATTAACGCGGCCGCGGTATGATTGTTTATTGCTGAGGACGTTACTCAAGGCGCCAGGCCGCTTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTa                                                                                                                                                             <  1:142919/250‑1 (MQ=255)                                                                                       tGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTAt                                                                                                                                                            >  1:25096/1‑251 (MQ=255)                                                                                                             ggggTATGATTTTTTATTGCTGAAACCGTTTCTCAAGTCGCCAGGCCGCGTTTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTggg                                                                                                                                       <  1:58798/249‑1 (MQ=255)                                                                                                                gTATGAGTTTTTATTGCTGAAGCCGTTCCTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTa                                                                                                                                    <  1:83627/251‑1 (MQ=255)                                                                                                                         tttATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACt                                                                                                                           >  1:271535/1‑251 (MQ=255)                                                                                                                                   aggacgttcttcaagcggccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcggtg                                                                                                                 <  1:163376/230‑1 (MQ=255)                                                                                                                                        gTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcggtggcggc                                                                                                            >  1:83872/1‑248 (MQ=255)                                                                                                                                         ttACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGcc                                                                                                           >  1:30696/1‑251 (MQ=255)                                                                                                                                           attcaaatcgccaggctgcgTCTGGCCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGTCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTg                                                                                                         <  1:263395/234‑1 (MQ=255)                                                                                                                                                     ccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTg                                                                                               >  1:173635/4‑251 (MQ=255)                                                                                                                                                      caggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGGCTGGTTCGGACTgg                                                                                              >  2:80363/3‑251 (MQ=255)                                                                                                                                                      caggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTgg                                                                                              <  1:231802/249‑1 (MQ=255)                                                                                                                                                                           aGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTa                                                                        <  1:80363/251‑1 (MQ=255)                                                                                                                                                                             cAACTGATGTATAGCATATGGGAAGACGCGCAGTACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGGCAAGCTGCGCGCCAGCAGCCCCGATCTTTCCCCGAATAATACTCAACGCTGCATGGGATATATTCTGAGGGGCCTGTAATGGTTATCCGCAGGCGGTTGTTGCTGGGCTAGTTGTTACTGGTGGCGGGGGCAGCCGGGTTCTTACTGGCCATTTTTGTCAAAG‑AAGTTaaa                                                                      <  1:15879/251‑1 (MQ=255)                                                                                                                                                                                     ggATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCgg                                                              >  1:137866/1‑250 (MQ=255)                                                                                                                                                                                            gtgtGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGcc                                                                                                                                                                        <  1:313070/136‑1 (MQ=255)                                                                                                                                                                                            gtgtGGGAAGAGGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGGTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCTCGGCATGGGCTATAGCCTGAGGGGcc                                                                                                                                                                        >  2:313070/4‑139 (MQ=255)                                                                                                                                                                                               tGGGAAGATGCGCAGGACAACTACGATCGCACCGATGATACCCCCATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGAATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTTAAAGTAAGT‑AAACCGGGCGTGCGAGGAGc                                                     >  2:253636/1‑250 (MQ=255)                                                                                                                                                                                                         cgcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAggggg                                          >  1:250591/1‑250 (MQ=255)                                                                                                                                                                                                          gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttttt                                                                                                                              >  1:253586/1‑167 (MQ=255)                                                                                                                                                                                                          gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttttt                                                                                                                              <  2:253586/167‑1 (MQ=255)                                                                                                                                                                                                                 ccctccgaccgccccgtgtaTACCCACATTAAAACGTGGCGTGCCAAGTTGCGCGCCATCAACCCCGATCTGTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGTGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGTGGACGTTAATc                                  <  2:250591/232‑1 (MQ=255)                                                                                                                                                                                                                                  tgATACCCACATTAAAACG‑TTCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTg                >  1:74948/2‑251 (MQ=255)                                                                                                                                                                                                                                              tAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTGTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTgg      <  2:67123/250‑1 (MQ=255)                                                                                                                                                                                                                                                 aaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGcgcg  >  1:245171/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        CGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCG  >  NC_000913/4636367‑4636856

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