Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,672,639	T→C	L227P (CTG→CCG)	ydgD →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,672,639	0	T	C	100.0%	191.2 / NA	53	L227P (CTG→CCG)	ydgD	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/26);  total (27/26)

AGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGT  >  NC_000913/1672510‑1672770                                                                                                                                  |                                                                                                                                    aggcaggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCa                                                                                                                            <  1:293661/139‑1 (MQ=255)                 aGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTgg                                                                                                            >  2:122491/1‑139 (MQ=255)                 aGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTgg                                                                                                            >  2:1155536/1‑139 (MQ=255)                 aGATCATCTCGATACGTTGTACAGTCATCAAAACTGAGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTgg                                                                                                            <  2:1012355/139‑1 (MQ=255)                     catcTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGc                                                                                                               >  2:671533/1‑132 (MQ=255)                     catcTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGc                                                                                                               <  1:671533/132‑1 (MQ=255)                     catcTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCaat                                                                                                        >  1:1447140/1‑139 (MQ=255)                        ctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGAt                                                                                                                                  >  1:1311459/1‑110 (MQ=255)                        ctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGAt                                                                                                                                  <  2:1311459/110‑1 (MQ=255)                        ctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCaattaa                                                                                                     <  1:1479821/139‑1 (MQ=255)                          cGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCaattaatt                                                                                                   <  2:1117593/139‑1 (MQ=255)                           gATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTg                                                                                                  <  2:1035319/139‑1 (MQ=255)                                  tGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATa                                                                                                                          >  1:1481764/1‑108 (MQ=255)                                  tGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATa                                                                                                                          <  2:1481764/108‑1 (MQ=255)                                  tGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTggg                                                                                                <  2:307382/134‑1 (MQ=255)                                  tGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTggg                                                                                                >  1:307382/1‑134 (MQ=255)                                           atcaAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCa                                                                                                          >  2:660244/1‑115 (MQ=255)                                           atcaAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCa                                                                                                          <  1:660244/115‑1 (MQ=255)                                           atcaAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc                                                                                  >  1:723509/1‑139 (MQ=255)                                           atcaAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc                                                                                  >  2:296830/1‑139 (MQ=255)                                           atcaAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc                                                                                  >  1:1110856/1‑139 (MQ=255)                                           atcaAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTGGCCGGGTGACAGCGGTTCGCCTCCGATGTTGGATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTAGGc                                                                                  >  1:1047762/1‑139 (MQ=255)                                                aaCTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTg                                                                             <  2:1110856/139‑1 (MQ=255)                                                   tgtgGAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCcg                                                                          >  1:1065351/1‑139 (MQ=255)                                                          tGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTg                                                                                                  >  2:284188/1‑108 (MQ=255)                                                          tGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTg                                                                                                  <  1:284188/108‑1 (MQ=255)                                                            aCTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATc                                                                 >  1:726701/1‑139 (MQ=255)                                                            aCTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATc                                                                 >  1:517050/1‑139 (MQ=255)                                                                 ctggGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATg                                                                                                                    <  2:435533/83‑1 (MQ=255)                                                                 ctggGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATg                                                                                                                    >  1:435533/1‑83 (MQ=255)                                                                  tggGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgc                                                           <  2:954613/139‑1 (MQ=255)                                                                  tggGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgc                                                           <  2:1447140/139‑1 (MQ=255)                                                                  tggGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgc                                                           <  1:604665/139‑1 (MQ=255)                                                                   gggCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTACGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgcg                                                          <  1:1155536/139‑1 (MQ=255)                                                                    ggCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCcgc                                                                         >  2:482329/1‑123 (MQ=255)                                                                    ggCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCcgc                                                                         <  1:482329/123‑1 (MQ=255)                                                                      cgcAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgcgcgc                                                       >  2:494744/1‑139 (MQ=255)                                                                                                tGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCg                             >  1:752282/1‑139 (MQ=255)                                                                                                 gCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCt                                                                                                              >  2:1591213/1‑57 (MQ=255)                                                                                                 gCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCt                                                                                                              <  1:1591213/57‑1 (MQ=255)                                                                                                 gCGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCgg                            >  1:563281/1‑139 (MQ=255)                                                                                                  cGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCg                                              >  2:260814/1‑120 (MQ=255)                                                                                                  cGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCg                                              <  1:260814/120‑1 (MQ=255)                                                                                                  cGATACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTAACGGt                           >  1:598669/1‑139 (MQ=255)                                                                                                     tACCTTGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGtttt                        <  1:296830/139‑1 (MQ=255)                                                                                                          tGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgcg                                                          >  1:860860/1‑100 (MQ=255)                                                                                                          tGCCGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgcg                                                          <  2:860860/100‑1 (MQ=255)                                                                                                             cGGGTGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACaa                <  1:682514/139‑1 (MQ=255)                                                                                                                 tGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCaa                                                                                          >  1:473655/1‑61 (MQ=255)                                                                                                                 tGACAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCaa                                                                                          <  2:473655/61‑1 (MQ=255)                                                                                                                   aCAGCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGa          <  1:1020793/139‑1 (MQ=255)                                                                                                                     aGCGGTTCGCCTCCGATTTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATc        <  2:1065351/139‑1 (MQ=255)                                                                                                                      gCGGTTCGCCTCCGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCa       >  1:1165673/1‑139 (MQ=255)                                                                                                                           tCGCCTACGATCTTGAATACCGATGACGGCTGGCAATTAATTGTGGTGCAAAGTTCGGCTCCTGCCGCGACAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGTTCAACTGt  <  2:1047762/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                    AGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGT  >  NC_000913/1672510‑1672770

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