Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,299	T→C	G374G (GGT→GGC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,299	0	T	C	92.0%	54.1 / ‑1.4	25	G374G (GGT→GGC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (7/16);  total (8/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAAT  >  NC_000913/4639052‑4639541                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cggACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTgg                                                                                                                                                                                                                                                 <  1:72277/249‑1 (MQ=255)                                gATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCtgt                                                                                                                                                                                                                  <  1:130521/251‑1 (MQ=255)                                                aTGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGgtgg                                                                                                                                                                                                  >  1:75335/1‑250 (MQ=255)                                                  ggCGTTCTTTTTTTTGGAAACGTTCACCGCGCATCGTTTACTCCCAATGCAATATTTGCTGGTGGGGCTTTCATTTGTGATGTTTTATTTGCTCTTGATGGCGCTTTCTGAACATATCGGTGTTACCGTGGCATGTATAATCGCCAGTCTGATTGGGGCGTGAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGGTTACCCTCGCGCTGTTGTTGCTGGATGGtgtg                                                                                                                                                                                                 <  2:235578/250‑1 (MQ=255)                                                           tGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATTTTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGa                                                                                                                                                                                        <  1:54217/250‑1 (MQ=255)                                                               tttGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCt                                                                                                                                                                                   >  1:168824/1‑251 (MQ=255)                                                                            ccGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGGAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGCTGGTGTGATGTGGGGACTGCTCAACTCTGCCGAt                                                                                                                                                                      >  1:308999/1‑251 (MQ=255)                                                                                    cGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTTCTGGATGGTGTGATGTGGGGACTGCTCAATTCTGCCGATAGCGCGCt                                                                                                                                                              <  1:104411/251‑1 (MQ=255)                                                                                     gTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGGTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGTATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtg                                                                                                                                                             <  1:44838/251‑1 (MQ=255)                                                                                        tACACCCAATGCAATATTTGCTGGTGGGTCTTTCATTGGGTATGTGTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGTATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCtgtt                                                                                                                                                           <  2:175734/250‑1 (MQ=255)                                                                                          cacCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTgttgtt                                                                                                                                                        <  1:93950/251‑1 (MQ=255)                                                                                                aTGCAATCTTTGCTGGTGGGGCTTTCATTGGCGAAGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTTAAAGGCTGGCGCAACAGCCTGTTGTTTACCCTCGCGCTGTTGTTGCTGGGTGGTGTGATGTGGGGACTGCTCAACTCTGCCGATTACGGGCTGTTGTTgggta                                                                                                                                                   >  1:337906/1‑248 (MQ=255)                                                                                                   cAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAg                                                                                                                                               <  1:7001/251‑1 (MQ=255)                                                                                                                               tGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGGATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGAGGGTGTGGGGTGGGGGCTGTACAACTCTGTCGATAGCGCGCTGTTGTTGGGGACAAGTGTGCtaggggggggtctgggcggactg                                                                                                                   >  2:299231/1‑229 (MQ=255)                                                                                                                                                      gcgcTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGTGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTtgt                                                                                                           <  1:41039/236‑1 (MQ=255)                                                                                                                                                      gcgcTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGa                                                                                            <  1:218435/251‑1 (MQ=255)                                                                                                                                                                         gTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGAAGGCGGTATTGAAAGGCTTGCGCGACAGCACGTTGTTTAACCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAACGCGCGGGTTTTGGGAACCAGTGTGCTGGGGGTGGCGCTGGCCGGCATGGTGTTTTGGACTCGTAATATCGACTGGTATGGGTTTTCACTg                                                                         >  1:103121/1‑251 (MQ=255)                                                                                                                                                                             tacgggtgCATGTATAAACGCTAGTCTGATAGGGGCGTTAATTAACGGTATTCATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGTTTGTTGTTTACCCGCGCGCGTTTGTTGCTGGATGGTTTGATGTGGGGATTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGTTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACGCGTAATATCGACTTGTATGCGTTTTTACTGCCg                                                                      <  2:343827/243‑1 (MQ=255)                                                                                                                                                                                                gCCAGTCTGATTGGGGCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAaaga                                                  <  1:111018/251‑1 (MQ=255)                                                                                                                                                                                                        taTTGTGGCGTTAATGAACGGTATTTATTTGTAGGCGGTATGGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGTATGGTGTGATGTGGTGTCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATTATGGTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTATAGAAGTTa                                             <  1:42803/247‑1 (MQ=255)                                                                                                                                                                                                         agttgggcattAATGACCGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTCCAACAGCATGTTGTTTACCCTCGCGCTGTGGTTGCTGGATGGTGTGATGTGGGGATTGCTCATCTCCTCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTTGCGCTGGCCGGCATGATGTTTGTGACCCGTAATAACGACTGTTATGCGTTTTCCCTGCCGAAAATGAAAGCCAGTAAAGAAGATACaa                                          <  2:93819/240‑1 (MQ=255)                                                                                                                                                                                                            ggggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGTCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGTTGTGATGTGGGTACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATAACGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGa                                      <  1:222128/251‑1 (MQ=255)                                                                                                                                                                                                             gggCGTTAATGAACGGTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGAc                                     <  1:74858/251‑1 (MQ=255)                                                                                                                                                                                                                       gcaCGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAAAAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGAGATGTGGGGACTGCACACCTCTGCCGATAGCGCGCTGTTGTTGGGAAACAGTGTGCTGGGGGTGGCGATGGCCGGCATGATGATTGGGACCCGTAATATCGACGGGTGTGCGTTTTAACTTCCGAACTCGAAAGCCAGAAAAGAAGGTAAAACGGaaggttctgtgc                           >  2:101945/3‑240 (MQ=255)                                                                                                                                                                                                                                                 tGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGGGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCTGCATGATGTTTGTGACTCGTAATATCGACTGGGATGCGTTTTCACTGCCGAAAATGAAAGCCAGGAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAAt  >  1:260180/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAAT  >  NC_000913/4639052‑4639541

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