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

GTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACT  >  NC_000913/4636268‑4636736                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgccagg                                                                                                                                                                                                                            >  1:240282/1‑245 (MQ=255)                                    cgacTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCg                                                                                                                                                                                         >  1:208060/1‑251 (MQ=255)                                                      gTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTCCTCAAGTCGCCAGGCCGAGTCTTGTACCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAgg                                                                                                                                                                       >  2:83627/1‑251 (MQ=255)                                                             ccccGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGCTAGCGTGTGGGGAGATGCGCAGGacacat                                                                                                                                                                 >  2:352250/1‑248 (MQ=255)                                                                    aaGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGc                                                                                                                                                         >  1:207666/1‑251 (MQ=255)                                                                       tgtgCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCAcc                                                                                                                                                      >  1:87669/1‑251 (MQ=255)                                                                        gtgCGCCAGGGTGCCCACCTTACGGCGTCGGGTGAAGATGTTCTCGTCGCCGTCTCCCGTCATCCGTTTTGGTCATTTTGAATTGAATAACCCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGAGGCGCAGGACACCTACGATCGCACCg                                                                                                                                                     <  1:105500/251‑1 (MQ=255)                                                                                gggTGCGCACCTTACTGCGGCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTTTTGGTCATTTTGAATGGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAcc                                                                                                                                             <  1:352250/251‑1 (MQ=255)                                                                                          cTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATTGATTGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAACGCTGACGCGGCATGAGTTTTTATTGATGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTACTGGATGGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACAtaaat                                                                                                                                    >  1:224293/1‑246 (MQ=255)                                                                                                    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)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGA‑CACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACT  >  NC_000913/4636268‑4636736

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