Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,083,620	G→A	R60C (CGT→TGT)	speB ←	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,083,620	0	G	A	100.0%	181.7 / NA	57	R60C (CGT→TGT)	speB	agmatinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/22);  total (35/22)

GCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGGCA  >  NC_000913/3083485‑3083753                                                                                                                                        |                                                                                                                                      gCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    <  2:568704/139‑1 (MQ=255)    tttCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    <  2:828502/136‑1 (MQ=255)    tttCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    >  1:828502/1‑136 (MQ=255)    tttCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGc                                                                                                                                 >  1:1500638/1‑139 (MQ=255)     ttCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCt                                                                                                                                >  1:283423/1‑139 (MQ=255)      tCGCTCATCTCACGGGCATCGCCAAATGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTg                                                                                                                               <  1:79269/139‑1 (MQ=255)            aTCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGAc                                                                                                                         >  2:980931/1‑139 (MQ=255)            aTCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGAc                                                                                                                         >  1:1228131/1‑139 (MQ=255)            aTCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGAc                                                                                                                         >  2:1152467/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    <  1:912984/116‑1 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    >  2:912984/1‑116 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                             >  1:536037/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                             >  1:1463040/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                             >  2:1100492/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                             >  1:734358/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                             >  1:766605/1‑139 (MQ=255)                                 gCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                    <  1:1242075/139‑1 (MQ=255)                                 gCATATACCAGAGCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                    <  1:141204/139‑1 (MQ=255)                                    tataCCAGATCGCCGCAGTCCACGACGTTCAGACGTTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAg                                                                                                 >  1:136849/1‑139 (MQ=255)                                      taCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGaa                                                                                               >  2:447927/1‑139 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCg                                                                                                                                  >  2:867891/1‑99 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCg                                                                                                                                  <  1:867891/99‑1 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                    <  1:391012/129‑1 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                    >  2:391012/1‑129 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                          >  2:1218400/1‑139 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                          <  2:536037/139‑1 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                          >  2:367109/1‑139 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                          >  1:1274128/1‑139 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                          >  1:65865/1‑139 (MQ=255)                                              cgccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGTCCACCCGCACGACCAGAAGTGGCCat                                                                                       >  2:1053570/1‑139 (MQ=255)                                                ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGTCatat                                                                                     >  1:391937/1‑139 (MQ=255)                                                 cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATc                                                                                    >  1:1185961/1‑139 (MQ=255)                                                      tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACg                                                                               >  1:325351/1‑139 (MQ=255)                                                             cgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                                        <  1:1100492/139‑1 (MQ=255)                                                             cgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                                        <  2:766605/139‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTCCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                                     >  2:1203783/1‑139 (MQ=255)                                                                 cAGACGTTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTa                                                                    <  2:136849/139‑1 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                         <  1:579487/96‑1 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                         >  2:579487/1‑96 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAcc                                                              >  1:727407/1‑139 (MQ=255)                                                                         cACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCa                                                            >  2:497270/1‑139 (MQ=255)                                                                                 tCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCAt                                                    >  2:799350/1‑139 (MQ=255)                                                                                  cGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATc                                                   >  1:1172355/1‑139 (MQ=255)                                                                                   gAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGTCAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAAAAGTGGCCATCTCGAACGGCACGCCAGTAATCACCCAGTCTGCATCg                                                  <  1:1053570/139‑1 (MQ=255)                                                                                    aaaTTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGc                                                 <  2:734358/139‑1 (MQ=255)                                                                                    aaaTTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGc                                                 <  2:1463040/139‑1 (MQ=255)                                                                                        tCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTc                                             >  1:373699/1‑139 (MQ=255)                                                                                        tCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTc                                             <  2:64553/139‑1 (MQ=255)                                                                                          cACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCAt                                           <  1:497270/139‑1 (MQ=255)                                                                                          cACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCAt                                           >  1:526101/1‑139 (MQ=255)                                                                                          cACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCAt                                           <  2:283423/139‑1 (MQ=255)                                                                                           aCGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                                          >  2:483000/1‑139 (MQ=255)                                                                                                  gCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTgg                                   <  1:1218400/139‑1 (MQ=255)                                                                                                                aGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAg                     >  2:364985/1‑139 (MQ=255)                                                                                                                        ttCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTaaaa             >  2:95468/1‑139 (MQ=255)                                                                                                                               aaaCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCaaa      <  2:745094/139‑1 (MQ=255)                                                                                                                                   cTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGGCa  <  1:95468/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      GCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGGCA  >  NC_000913/3083485‑3083753

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