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

TTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGG  >  NC_000913/3083488‑3083751                                                                                                                                     |                                                                                                                                    tttCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGc                                                                                                                               >  1:258638/1‑139 (MQ=255)                     tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                           >  1:123756/1‑139 (MQ=255)                     tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                           >  1:291575/1‑139 (MQ=255)                     tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                           >  1:964282/1‑139 (MQ=255)                         cAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       <  2:717462/139‑1 (MQ=255)                                 tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAg                                                                                               >  2:344584/1‑139 (MQ=255)                                 tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAg                                                                                               >  2:146033/1‑139 (MQ=255)                                        gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacca                                                                                                            <  1:666284/119‑1 (MQ=255)                                        gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacca                                                                                                            >  2:666284/1‑119 (MQ=255)                                        gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                        >  2:102599/1‑139 (MQ=255)                                           cgccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCat                                                                                     >  2:1001510/1‑139 (MQ=255)                                             ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCatat                                                                                   >  2:315581/1‑139 (MQ=255)                                             ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCatat                                                                                   >  2:741701/1‑139 (MQ=255)                                              cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATc                                                                                  >  2:575436/1‑139 (MQ=255)                                                    ccACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                                            <  2:720739/139‑1 (MQ=255)                                                      acgacgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCa                                                                          <  2:696100/139‑1 (MQ=255)                                                      acgacgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCa                                                                          <  1:344584/139‑1 (MQ=255)                                                          cgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                                      <  2:677210/139‑1 (MQ=255)                                                          cgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                                      <  1:146033/139‑1 (MQ=255)                                                             tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                                   >  2:327763/1‑139 (MQ=255)                                                                               cGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                                            <  2:845065/112‑1 (MQ=255)                                                                               cGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACgg                                                                            >  1:845065/1‑112 (MQ=255)                                                                                gAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       <  1:741404/84‑1 (MQ=255)                                                                                gAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       >  2:741404/1‑84 (MQ=255)                                                                                   aTTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCg                                                >  2:181562/1‑136 (MQ=255)                                                                                   aTTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCg                                                <  1:181562/136‑1 (MQ=255)                                                                                        aCGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAccc                                                           >  2:175743/1‑120 (MQ=255)                                                                                        aCGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAccc                                                           <  1:175743/120‑1 (MQ=255)                                                                                        aCGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                                        <  2:291575/139‑1 (MQ=255)                                                                                                          cccAGGCCAGATTCGTCGAAACCTGACAGATGGCTGCCGGa                                                                                                                        <  1:708160/41‑1 (MQ=255)                                                                                                          cccAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                        >  2:708160/1‑41 (MQ=255)                                                                                                           ccAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGc                     <  2:144160/139‑1 (MQ=255)                                                                                                                   gATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTaa             <  1:848545/139‑1 (MQ=255)                                                                                                                         tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAAcc       >  2:158990/1‑139 (MQ=255)                                                                                                                         tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAAcc       >  2:753822/1‑139 (MQ=255)                                                                                                                              aCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAgg  <  1:327763/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    TTTCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGG  >  NC_000913/3083488‑3083751

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