Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,303	T→C	C376R (TGC→CGC)	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,303	0	T	C	91.7%	49.6 / ‑5.5	24	C376R (TGC→CGC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (7/15);  total (8/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 = 6.90e-01

GATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAAT  >  NC_000913/4639083‑4639541                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               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)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               GATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATGAACGGTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAAT  >  NC_000913/4639083‑4639541

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