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%	212.8 / NA	65	R60C (CGT→TGT)	speB	agmatinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/29);  total (36/29)

CTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGGCA  >  NC_000913/3083493‑3083753                                                                                                                                 |                                                                                                                                      cTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCg                                                                                                                            <  1:228767/139‑1 (MQ=255)           gggCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCg                                                                                                                  >  2:1266944/1‑139 (MQ=255)            ggCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGa                                                                                                                 >  1:1583986/1‑139 (MQ=255)              cATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    >  1:898139/1‑118 (MQ=255)              cATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                    <  2:898139/118‑1 (MQ=255)                tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGa                                                                                                                 <  1:1211858/135‑1 (MQ=255)                tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGa                                                                                                                 >  2:1211858/1‑135 (MQ=255)                    cAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                         <  2:404434/139‑1 (MQ=255)                       aGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACg                                                                                                      >  1:572907/1‑139 (MQ=255)                           atatACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCa                                                                                                  >  2:128188/1‑139 (MQ=255)                                 cAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGt                                                                                             >  2:1605293/1‑139 (MQ=255)                                   gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                          >  2:235650/1‑107 (MQ=255)                                   gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                          <  1:235650/107‑1 (MQ=255)                                   gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAg                                                                                              <  2:913860/135‑1 (MQ=255)                                   gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAg                                                                                              >  1:913860/1‑135 (MQ=255)                                   gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                          >  1:555192/1‑139 (MQ=255)                                    aTCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGcc                                                                                         >  1:1413812/1‑139 (MQ=255)                                    aTCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGcc                                                                                         >  1:311370/1‑139 (MQ=255)                                    aTCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGcc                                                                                         >  2:1292284/1‑139 (MQ=255)                                         cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                          <  2:392476/101‑1 (MQ=255)                                         cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                          >  1:392476/1‑101 (MQ=255)                                         cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATc                                                                                    >  1:167462/1‑139 (MQ=255)                                         cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATc                                                                                    >  2:1283298/1‑139 (MQ=255)                                         cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATc                                                                                    >  1:1450797/1‑139 (MQ=255)                                           cAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                         >  2:1488397/1‑116 (MQ=255)                                           cAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                         <  1:1488397/116‑1 (MQ=255)                                              tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACg                                                                               >  2:1035136/1‑139 (MQ=255)                                                    acgTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                                         <  2:366373/139‑1 (MQ=255)                                                    acgTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                                         <  1:965172/139‑1 (MQ=255)                                                    acgTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                                         <  2:555192/139‑1 (MQ=255)                                                        tCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                                     >  2:396929/1‑139 (MQ=255)                                                                 cACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCa                                                            <  1:1035136/139‑1 (MQ=255)                                                                 cACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCa                                                            <  1:1292284/139‑1 (MQ=255)                                                                        gTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTg                                                                                                                     >  1:1361633/1‑75 (MQ=255)                                                                        gTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTg                                                                                                                     <  2:1361633/75‑1 (MQ=255)                                                                        gTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCa                                                     <  1:735341/139‑1 (MQ=255)                                                                        gTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCa                                                     >  1:648320/1‑139 (MQ=255)                                                                        gTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCa                                                     <  1:467252/139‑1 (MQ=255)                                                                            aaaTTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGt                                                          <  1:133908/130‑1 (MQ=255)                                                                            aaaTTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGt                                                          >  2:133908/1‑130 (MQ=255)                                                                                   aC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGGCCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                                          <  1:1283298/139‑1 (MQ=255)                                                                                       ggAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTg                                                                                                                               >  1:1155665/1‑51 (MQ=255)                                                                                       ggAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTg                                                                                                                               <  2:1155665/51‑1 (MQ=255)                                                                                             ggTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGaa                                 <  2:572907/139‑1 (MQ=255)                                                                                              gTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                                          <  2:991678/129‑1 (MQ=255)                                                                                              gTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                                          >  1:991678/1‑129 (MQ=255)                                                                                                 gtgtTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCa                                                                >  2:1229146/1‑104 (MQ=255)                                                                                                 gtgtTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCa                                                                <  1:1229146/104‑1 (MQ=255)                                                                                                 gtgtTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTc                             >  2:205722/1‑139 (MQ=255)                                                                                                              aGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTa                >  1:1394230/1‑139 (MQ=255)                                                                                                              aGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTa                >  1:627453/1‑139 (MQ=255)                                                                                                                aTTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGc                                                                                                                   >  1:65240/1‑38 (MQ=38)                                                                                                                aTTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGc                                                                                                                   <  2:65240/38‑1 (MQ=38)                                                                                                                aTTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTaaa              >  1:799715/1‑139 (MQ=255)                                                                                                                aTTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTaaa              >  1:440126/1‑139 (MQ=255)                                                                                                                 ttCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAACAGAAGTGGCCATATAGGACGGCACGCCAGTAAGCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGGTCATCGGCAGGCGTaaaa             >  1:996794/1‑139 (MQ=255)                                                                                                                     tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGcc                                                                                         <  2:539562/59‑1 (MQ=255)                                                                                                                     tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGcc                                                                                         >  1:539562/1‑59 (MQ=255)                                                                                                                     tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                                                 >  1:1487639/1‑67 (MQ=255)                                                                                                                     tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                                                 <  2:1487639/67‑1 (MQ=255)                                                                                                                     tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAAcc         <  2:440126/139‑1 (MQ=255)                                                                                                                         aaCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAg     <  2:1007050/139‑1 (MQ=255)                                                                                                                         aaCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAg     >  1:1007050/1‑139 (MQ=255)                                                                                                                          aCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAgg    <  1:483287/139‑1 (MQ=255)                                                                                                                            cTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGGCa  <  1:1266944/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      CTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCAC‑GGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGGCA  >  NC_000913/3083493‑3083753

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