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

GCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAA  >  NC_000913/3083485‑3083749                                                                                                                                        |                                                                                                                                  gCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                <  2:1266881/139‑1 (MQ=255)                      cATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacc                                                                                                           >  2:53582/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                         >  2:1301921/1‑139 (MQ=255)                        tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                         >  1:282209/1‑139 (MQ=255)                            cAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                     <  1:1106087/139‑1 (MQ=255)                                 gCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                <  2:1156242/139‑1 (MQ=255)                                    tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                <  2:703945/104‑1 (MQ=255)                                    tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                >  1:703945/1‑104 (MQ=255)                                          aGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTgg                                                                                       <  1:141268/139‑1 (MQ=255)                                           gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                      >  1:667339/1‑139 (MQ=255)                                             tCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCa                                                                                    >  2:102878/1‑139 (MQ=255)                                              cgccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCAGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCat                                                                                   >  2:271595/1‑139 (MQ=255)                                                ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                >  1:377247/1‑92 (MQ=255)                                                ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                <  2:377247/92‑1 (MQ=255)                                                    aGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                                             <  2:667339/139‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                                            <  1:251740/96‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                                            >  2:251740/1‑96 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                                            >  2:469971/1‑96 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                                            <  1:469971/96‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                                          >  1:791720/1‑130 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                                          <  2:791720/130‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGc                                                                         <  2:1017755/131‑1 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGc                                                                         >  1:1017755/1‑131 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                                     >  2:974170/1‑135 (MQ=255)                                                                tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                                     <  1:974170/135‑1 (MQ=255)                                                                         cACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCgg                                                                                                     <  1:1417509/94‑2 (MQ=255)                                                                         cACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                     >  2:1417509/1‑94 (MQ=255)                                                                                            cGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATAc                                     >  2:339606/1‑139 (MQ=255)                                                                                                  gCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTgg                               <  2:282209/139‑1 (MQ=255)                                                                                                     gTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCa                                                                        >  2:559186/1‑95 (MQ=255)                                                                                                     gTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCa                                                                        <  1:559186/95‑1 (MQ=255)                                                                                                           ttCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATc                      >  2:667525/1‑139 (MQ=255)                                                                                                                aGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAg                 >  2:1295848/1‑139 (MQ=255)                                                                                                                              gAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCaa   <  1:1295848/139‑1 (MQ=255)                                                                                                                               aaaCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCaaa  <  1:53582/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GCTTTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAA  >  NC_000913/3083485‑3083749

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