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

GCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAG  >  NC_000913/3083485‑3083734                                                                                                                                         |                                                                                                                   gCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                 <  1:1453199/139‑1 (MQ=255)     ttCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCt                                                                                                             >  2:770633/1‑139 (MQ=255)         cTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCg                                                                                                         <  2:792979/139‑1 (MQ=255)              ctcACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCg                                                                                                    <  1:1214517/139‑1 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCg                                                                                                               <  1:54452/118‑1 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCg                                                                                                               >  2:54452/1‑118 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCg                                                                                                    >  1:549166/1‑129 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCg                                                                                                    <  2:549166/129‑1 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                          >  2:781722/1‑139 (MQ=255)                               aGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACg                                                                                   >  2:374711/1‑139 (MQ=255)                                 gCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                 <  1:770633/139‑1 (MQ=255)                                    tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATc                                                                                                                <  2:663371/105‑1 (MQ=255)                                    tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATc                                                                                                                >  1:663371/1‑105 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                 >  1:1532519/1‑98 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                 <  2:1532519/98‑1 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCTAGATTCGTCGAAACCTGACAGAt                                                                                                                 <  2:909060/97‑1 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCTAGATTCGTCGAAACCTGACAGAt                                                                                                                 >  1:909060/1‑97 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                       >  1:1122508/1‑139 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                       >  1:1104322/1‑139 (MQ=255)                                                ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGCAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCatat                                                                  >  2:1335493/1‑139 (MQ=255)                                                ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTgg                                                                        <  2:399311/133‑1 (MQ=255)                                                ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTgg                                                                        >  1:399311/1‑133 (MQ=255)                                                     gTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAAc                                                             <  1:687036/139‑1 (MQ=255)                                                      tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                       >  1:1115639/1‑96 (MQ=255)                                                      tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACg                                                            >  2:466161/1‑139 (MQ=255)                                                      tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATACGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                       <  2:1115639/96‑1 (MQ=255)                                                       ccACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                           <  2:1520630/139‑1 (MQ=255)                                                        cACGACGTTCAGACGCTCACGCATGTGGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGc                                                          >  2:174272/1‑139 (MQ=255)                                                           gacgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACg                                                       >  1:989665/1‑139 (MQ=255)                                                             cgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                     >  2:1146115/1‑139 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                             >  1:1642494/1‑96 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                             <  2:1642494/96‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                      >  1:762525/1‑103 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                      <  2:762525/103‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                  >  1:1717694/1‑139 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGAT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                      >  2:316352/1‑103 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGAT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                      <  1:316352/103‑1 (MQ=255)                                                                     cGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                      <  1:827074/130‑1 (MQ=255)                                                                     cGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                      >  2:827074/1‑130 (MQ=255)                                                                      gCTCACGCATGTCGAAAATCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCag                                            <  2:1104322/139‑2 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCa                                                                                     >  1:1451469/1‑97 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCa                                                                                     <  2:1451469/97‑1 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAcc                                           <  1:69803/139‑1 (MQ=255)                                                                       cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAcc                                           >  1:65422/1‑139 (MQ=255)                                                                                gTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCa                                  <  1:1377789/139‑1 (MQ=255)                                                                                 tCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCAt                                 <  1:1335493/139‑1 (MQ=255)                                                                                  cGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCg                                                                >  2:71957/1‑107 (MQ=255)                                                                                  cGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCg                                                                <  1:71957/107‑1 (MQ=255)                                                                                   gAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCg                               <  1:374711/139‑1 (MQ=255)                                                                                           aCGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                       <  1:781722/139‑1 (MQ=255)                                                                                           aCGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                       <  1:174272/139‑1 (MQ=255)                                                                                                  gCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                           <  1:705661/96‑1 (MQ=255)                                                                                                  gCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                           >  2:705661/1‑96 (MQ=255)                                                                                                         tgtT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCAt                                 >  2:711552/1‑115 (MQ=255)                                                                                                         tgtT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCAt                                 <  1:711552/115‑1 (MQ=255)                                                                                                         tgtT‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCa         >  2:218258/1‑139 (MQ=255)                                                                                                           tt‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                     >  2:1373266/1‑93 (MQ=255)                                                                                                           tt‑CCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                     <  1:1373266/93‑1 (MQ=255)                                                                                                              cccAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCgg     <  2:1717694/139‑1 (MQ=255)                                                                                                                 aGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAg  >  2:1798940/1‑139 (MQ=255)                                                                                                                   gCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                 <  2:704621/58‑1 (MQ=255)                                                                                                                   gCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                 >  1:704621/1‑58 (MQ=255)                                                                                                                      aGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                              >  1:221921/1‑74 (MQ=255)                                                                                                                      aGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                              <  2:221921/74‑1 (MQ=255)                                                                                                                             tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                      >  1:190266/1‑43 (MQ=255)                                                                                                                             tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                      <  2:190266/43‑1 (MQ=255)                                                                                                                                         |                                                                                                                   GCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTT‑CCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAG  >  NC_000913/3083485‑3083734

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