Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,386,502	T→C	S215P (TCC→CCC)	yfbL →	putative M28A family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,502	0	T	C	100.0%	124.3 / NA	35	S215P (TCC→CCC)	yfbL	putative M28A family peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/25);  total (10/25)

AACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGTCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCAGGGTTTATTCCCGGT  >  NC_000913/2386356‑2386648                                                                                                                                                   |                                                                                                                                                   aaCCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTg                                                                                                                                                >  2:431844/1‑151 (MQ=255)      cccTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCttt                                                                                                                                           >  2:606541/1‑151 (MQ=255)       ccTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTa                                                                                                                                          <  2:450180/151‑1 (MQ=255)          tcttTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATcc                                                                                                                                       >  1:150442/1‑151 (MQ=255)               cGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGAt                                                                                                                                   <  1:588196/150‑1 (MQ=255)                gTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCg                                                                                                                                 <  2:704982/151‑1 (MQ=255)                    cGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCgggg                                                                                                                              >  1:725390/1‑150 (MQ=255)                     gATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGgag                                                                                                                            <  1:34472/151‑1 (MQ=255)                     gATGAAATGGGGAGCGCGGAGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGCGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCgggg                                                                                                                              >  1:280322/1‑149 (MQ=255)                         aaaTGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACtt                                                                                                                        >  1:535012/1‑151 (MQ=255)                          aaTGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACttt                                                                                                                       <  1:487732/151‑1 (MQ=255)                               ggAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGc                                                                                                                  >  1:679890/1‑151 (MQ=255)                                           atgcaGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAg                                                                                                      <  2:288286/151‑1 (MQ=255)                                                      cTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGgag                                                                                                                            >  1:22902/1‑118 (MQ=255)                                                      cTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGgag                                                                                                                            <  2:22902/118‑1 (MQ=255)                                                         gAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATc                                                                                         <  1:92374/150‑1 (MQ=255)                                                               ccAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCg                                                                                  <  2:269926/151‑1 (MQ=255)                                                                             gATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTaaa                                                                     <  2:75365/150‑1 (MQ=255)                                                                                     tGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCAt                                                             >  2:681970/1‑150 (MQ=255)                                                                                      ggAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCAtt                                                            <  2:685414/150‑1 (MQ=255)                                                                                       gagaTGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCAttgt                                                          <  1:306602/151‑1 (MQ=255)                                                                                               tGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCAtc                                                   <  1:632769/150‑1 (MQ=255)                                                                                                 gCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCa                                                <  2:444344/151‑1 (MQ=255)                                                                                                    aTTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAg                                               <  1:171486/149‑1 (MQ=255)                                                                                                       aCGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGAtt                                           <  1:376153/150‑1 (MQ=255)                                                                                                            tcagcgcCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGAATTTATTGTCGTGGTCGGAAGAATACAGGATATCAACGCCGTACGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGACTTa                                         <  2:743953/144‑1 (MQ=255)                                                                                                                cgcCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTa                                 <  2:761936/151‑1 (MQ=255)                                                                                                                 gcCTGGAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTAt                                <  1:646901/151‑1 (MQ=255)                                                                                                                     ggAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTa                            <  2:8091/151‑1 (MQ=255)                                                                                                                     ggAAGCCAGAATTACCCTTACCCGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTa                            <  1:50365/151‑1 (MQ=255)                                                                                                                                 tacccttacccGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCa                 <  2:679849/150‑1 (MQ=255)                                                                                                                                  acccttacccGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCAg                <  1:290068/150‑1 (MQ=255)                                                                                                                                    ccttacccGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCAGGGt             <  1:102501/151‑1 (MQ=255)                                                                                                                                     cttacccGGCAATGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCAGGGtt            >  2:578642/1‑151 (MQ=255)                                                                                                                                               aaTGCCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCAGGGTTTATTCCCGGt  <  2:739675/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGACTCTGCGCCTGGAAGCCAGAATTACCCTTACCCGGCAATGTCCTGGCTTTATCCCGATCGGGGAGACTTTATTGCCGTGGTCGGCAGAATACAGGATATCAACGCCGTTCGTCAGGTAAAAGCGGCATTGTTGTCATCTCAGGATTTATCTGTTTATTCTATGAATACCCCAGGGTTTATTCCCGGT  >  NC_000913/2386356‑2386648

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