Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	231,803	C→T	P228S (CCG→TCG)	yafD →	endo/exonuclease/phosphatase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	231,803	0	C	T	100.0%	123.5 / NA	32	P228S (CCG→TCG)	yafD	endo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/19);  total (13/19)

CACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGG  >  NC_000913/231659‑231944                                                                                                                                                 |                                                                                                                                              caccacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGctcgctcg                                                                                                                                         <  1:560740/151‑1 (MQ=255)        gCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGt                                                                                                                                  <  2:477966/151‑1 (MQ=255)         cGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGtt                                                                                                                                 <  2:277570/151‑1 (MQ=255)                gTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTAcc                                                                                                                          <  2:25283/151‑1 (MQ=255)                   aTTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTg                                                                                                                       <  2:570700/151‑1 (MQ=255)                      aTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGt                                                                                                                     <  2:405063/150‑1 (MQ=255)                      aTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGt                                                                                                                     <  1:297823/150‑1 (MQ=255)                        ggCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTg                                                                                                                  >  2:582757/1‑151 (MQ=255)                         gCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGa                                                                                                                 <  2:70805/151‑1 (MQ=255)                                     aaTGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGa                                                                                                       >  1:546316/1‑149 (MQ=255)                                     aaTGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAg                                                                                                     <  1:262252/151‑1 (MQ=255)                                           tGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTg                                                                                               >  2:41290/1‑151 (MQ=255)                                            ggAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTg                                                                                               <  1:48680/150‑1 (MQ=255)                                             gAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTa                                                                                             <  2:500888/151‑1 (MQ=255)                                                   tAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGtt                                                                                       >  1:35675/1‑151 (MQ=255)                                                      aaGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcg                                                                                    <  2:306506/151‑1 (MQ=255)                                                       aGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgc                                                                                   <  1:106416/151‑1 (MQ=255)                                                               cgcgTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTcc                                                                           >  2:52592/1‑151 (MQ=255)                                                                  gTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGAt                                                                        >  2:333717/1‑151 (MQ=255)                                                                                cgcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACt                                                          <  1:565857/151‑1 (MQ=255)                                                                                 gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTCCCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTc                                                         <  2:80199/151‑1 (MQ=255)                                                                                   gggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGt                                                       >  2:302520/1‑151 (MQ=255)                                                                                    ggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGtt                                                      <  2:48900/151‑1 (MQ=255)                                                                                      aaaTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGa                                                    <  1:363389/151‑1 (MQ=255)                                                                                                    cAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGc                                       >  1:214950/1‑150 (MQ=255)                                                                                                         gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt                                 <  2:436970/151‑1 (MQ=255)                                                                                                           gTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGAGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGAAAGCCTg                                >  2:490520/1‑150 (MQ=255)                                                                                                                 ccGATGATCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaat                          >  1:3906/1‑150 (MQ=255)                                                                                                                     tgatCAGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAgg                      >  1:415398/1‑150 (MQ=255)                                                                                                                          aGCGCCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTc                >  1:188831/1‑151 (MQ=255)                                                                                                                           gcgcCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCa               <  2:602386/151‑1 (MQ=255)                                                                                                                             gcCGTCGGGCGTTTGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAgg             >  2:581874/1‑151 (MQ=255)                                                                                                                                         ttGGTCGCTCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAgg  <  1:404808/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              CACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGG  >  NC_000913/231659‑231944

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