Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,735	T→C	P186P (CCT→CCC)	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,638,735	0	T	C	100.0%	53.6 / NA	20	P186P (CCT→CCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/12);  total (8/12)

ATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTT  >  NC_000913/4638500‑4638969                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           ataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCgg                                                                                                                                                                                                                             >  2:215768/1‑251 (MQ=255)     cGGTCAGCTCTGGCACCGTGATTTAACGCTAAAAGCCGATTTAGCTGTTTCGCGTCTTCGCGCGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTCTTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCTGGCAGGGCGTGCATATCCCTTTAACCGATGGGGCATGGCGGAAGc                                                                                                                                                                                                                         >  2:166443/1‑251 (MQ=255)            gTCTGGCGCAGTGAGTTAACGCTAAAAGCCGATGTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCCTTACCCGAAGGGTACTTGCGGAAGCATAACCt                                                                                                                                                                                                                  <  2:232786/251‑1 (MQ=255)                          tttAACGCTAAAAGCAGATTTCAATGTTTCGCGTCTTAGCGAGCTCAACGCGCAAAATATCATCTTATGCAAGCCCTTTATTGTGATTAGCGCCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCCGGGTAGAGCTAGGCGGGCAGGGCGTGCATATACCTCTACCCGACGGGTACTGGCGGatgcgagacctaaagttgaaaaagc                                                                                                                                                                                                    >  2:120191/1‑227 (MQ=255)                                    aaaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAAttt                                                                                                                                                                                          <  1:7205/251‑1 (MQ=255)                                      aGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTaa                                                                                                                                                                                        <  1:35420/251‑1 (MQ=255)                                       gCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAg                                                                                                                                                                                       <  1:61258/251‑1 (MQ=255)                                               cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCg                                                                                                                                                                               <  2:60806/251‑1 (MQ=255)                                                gATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGGAGACGAAGATGGCCCCGAATTTAATGGGTACCgg                                                                                                                                                                              >  2:244545/1‑251 (MQ=255)                                                  tGTTTCGCGTCTTAGCGAGCTCAACGAGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGCAGTTAACGGAACGGCGCTGCACATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGAGTGCATATACCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATAATGCCCTACATTTAAGAGGTACCGGCg                                                                                                                                                                            >  1:142473/1‑251 (MQ=255)                                                                               cAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAATGGGACTGGCGGAAGCAGAACCTGGAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTTTTTGGTGCCTGCCGGGCGTact                                                                                                                                               >  2:65453/1‑249 (MQ=255)                                                                                      caggtgaggcaagagctttATTGTGAGGAGGGGCGGGGATGCGCGTGGTAGTGGTGTGGTGAAAGCGCCTGAAGGTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGAAGGGAGAGCATATCCCTGGACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAGGGACCTGAATTTAAGCGGTACAGGCGATCATTCTGTGGTGCCTGCCGGGCGGAATAGCGaaa                                                                                                                                        <  2:72179/235‑1 (MQ=255)                                                                                                        tATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGTTGTGGTGAAAGCGCCTTAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAAc                                                                                                                       <  1:116803/250‑1 (MQ=255)                                                                                                                                               ggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGTGTAGAGCAAGGCGTGCAGGGCGTGAATATCCCTATACCCGAAGGGGACTGGCGGAAGCAGAACATGAAGCTGAAGATGGCCCTGAAGTTAAGCGGTACCGGCGATCTTTCGGTGGTGCCTGCCGGGCGTTATAGTGAAATGTCTTTATCCCGCAACTGGTCGCATCCCAGTTTTTTAGGTGAtttgttgcaagcta                                                                               <  2:186689/251‑12 (MQ=255)                                                                                                                                                   aaaGCGCCTGAAGTTAACGGAACGGCGCTGTCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACg                                                                           <  1:167075/251‑1 (MQ=255)                                                                                                                                                                       aacggCGCTGACCATTGAACCCGGCAGCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGtt                                                       <  1:120191/251‑1 (MQ=255)                                                                                                                                                                          ggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGCAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCCTCCCAGTTTTTTAGGGGATTTTTTGCCAGCCAAACGGGAAGTTAGCGCGTAAGGTTTTc                                                    >  1:234414/1‑251 (MQ=255)                                                                                                                                                                                                    ggTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGttt                          <  1:135482/251‑1 (MQ=255)                                                                                                                                                                                                       tAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGAACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGGAATAGCGAAATGACCTTAACCAGCCACTGGCCGCATTCCAGTTTTTTAGGGGATTTTTTGCCAGCACAACGGGAAGTTAGCGGGTCAGGTTTTCAGGCGCAATGGCACAACAGCTGGGTTgtt                       >  2:152266/1‑249 (MQ=255)                                                                                                                                                                                                                             gCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGtttt  <  1:244545/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           ATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTT  >  NC_000913/4638500‑4638969

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