Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,017,078	A→G	F408L (TTC→CTC)	ycbZ ←	putative ATP‑dependent protease YcbZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,017,078	0	A	G	100.0%	193.3 / NA	54	F408L (TTC→CTC)	ycbZ	putative ATP‑dependent protease YcbZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/29);  total (25/29)

GAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGG  >  NC_000913/1016947‑1017210                                                                                                                                    |                                                                                                                                     gAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCa                                                                                                                              <  1:753104/140‑1 (MQ=255)  aGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCAt                                                                                                                             <  1:750590/140‑1 (MQ=255)  aGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCAt                                                                                                                             <  1:674225/140‑1 (MQ=255)  aGGGCGCTTATCAGGGCCCAGAATTCAGCCATCGACGCACTATCTCCACCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCAt                                                                                                                             <  2:220980/140‑1 (MQ=255)    ggCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATg                                                                                                                            >  2:517485/1‑139 (MQ=255)    ggCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATg                                                                                                                            <  1:517485/139‑1 (MQ=255)               gggCGCAGAGTTCAGCCATCGACGCTCTACCTCCATCCACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCt                                                                                                                <  1:220132/140‑1 (MQ=255)                ggCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGCCATCAGGAGCGCTTGCATGATCATCATCCCtt                                                                                                               <  2:564005/140‑1 (MQ=255)                  cgcAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGatcatcatc                                                                                                                   <  2:423130/134‑1 (MQ=255)                  cgcAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGatcatcatc                                                                                                                   >  1:423130/1‑134 (MQ=255)                    cAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGc                                                                                                           >  1:74470/1‑140 (MQ=255)                     agagTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCa                                                                                                          <  2:804602/140‑1 (MQ=255)                      gagTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCAt                                                                                                         >  2:66837/1‑140 (MQ=255)                                 tCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCttt                                                                                                              <  1:502761/124‑1 (MQ=255)                                 tCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCttt                                                                                                              >  2:502761/1‑124 (MQ=255)                                 tCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTgccg                                                                                              <  2:524243/140‑1 (MQ=255)                                        aCTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCt                                                                                       >  2:549825/1‑140 (MQ=255)                                             ctcCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCt                                                                                  <  2:509724/140‑1 (MQ=255)                                                  tCAACTTCACTGTATAACTGCTCAACTGTCAGCGATGATGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCg                                                                             <  1:210465/140‑1 (MQ=255)                                                    gagacagcTGTATGATTGCTCACATCTCAACGATGCTCAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGtt                                                                           <  1:221273/133‑1 (MQ=255)                                                         cACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATg                                                                      <  2:54439/140‑1 (MQ=255)                                                          aCTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGt                                                                     >  1:757667/1‑140 (MQ=255)                                                            tGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCg                                                                   <  2:74470/140‑1 (MQ=255)                                                              tATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCg                                                                   <  2:460867/138‑1 (MQ=255)                                                              tATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCg                                                                   <  2:466856/138‑1 (MQ=255)                                                              tATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCg                                                                   >  1:460867/1‑138 (MQ=255)                                                              tATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCg                                                                   >  1:466856/1‑138 (MQ=255)                                                                   cTGCTCAAATGTCAGTGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAAtt                                                            <  1:646280/140‑1 (MQ=255)                                                                   cTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAAtt                                                            >  1:820826/1‑140 (MQ=255)                                                                   cTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCAGCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAAtt                                                            >  2:756585/1‑140 (MQ=255)                                                                      cTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCAc                                                         >  2:79841/1‑140 (MQ=255)                                                                      cTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCAc                                                         >  1:7301/1‑140 (MQ=255)                                                                      cTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCAGGGATATTGGCGCTAAGGTCCGCTTGGCGTTGGATGTCGGTgctgtctc                                                         >  2:221273/1‑133 (MQ=255)                                                                       tCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCAcc                                                        <  1:745545/140‑1 (MQ=255)                                                                        cAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTCGCTTTCGATATCGGTGAATTCACCa                                                       >  1:340917/1‑140 (MQ=255)                                                                           aTGTCAGCGATGCTGTGAAGGGCATTTGTTGCTCAAGTTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCg                                                    <  1:273856/140‑1 (MQ=255)                                                                           aTGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCg                                                    <  1:716354/140‑1 (MQ=255)                                                                           aTGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGCCGGTGAATTCACCATCg                                                    <  2:496077/140‑1 (MQ=255)                                                                                aGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAat                                               >  1:633569/1‑140 (MQ=255)                                                                                     tGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGca                                          >  2:635812/1‑140 (MQ=255)                                                                                     tGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGca                                          >  1:773750/1‑140 (MQ=255)                                                                                     tGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGca                                          <  1:549825/140‑1 (MQ=255)                                                                                        tGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACaa                                       <  2:773750/140‑1 (MQ=255)                                                                                             aGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAg                                  <  2:716604/140‑1 (MQ=255)                                                                                                 gATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTaa                              >  1:485043/1‑140 (MQ=255)                                                                                                    cTGTTGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGc                           >  1:438783/1‑140 (MQ=255)                                                                                                        tGCTCAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCgaga                       <  1:756585/140‑1 (MQ=255)                                                                                                            cAAGCTGTAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGt                   >  1:540021/1‑140 (MQ=255)                                                                                                                   tAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCaa            <  2:633569/140‑1 (MQ=255)                                                                                                                   tAGTTCCGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCaa            <  1:79841/140‑1 (MQ=255)                                                                                                                         cGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAgcgc      >  2:562637/1‑140 (MQ=255)                                                                                                                         cGACATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAgcgc      >  2:655117/1‑140 (MQ=255)                                                                                                                           aCATCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGt    >  1:574101/1‑140 (MQ=255)                                                                                                                             aTCAGGAGCGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTgg  >  2:14591/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                     GAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGG  >  NC_000913/1016947‑1017210

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