Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,488,579	(C)7→8	coding (1485/1503 nt)	yjgR ←	DUF853 family protein with NTPase fold

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,488,572	1	.	C	100.0%	143.6 / NA	40	R498G (AGA→GGA)	yjgR	DUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (21/19);  total (21/19)

TTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑CCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAAT  >  NC_000913/4488444‑4488701                                                                                                                                  |                                                                                                                                  ttAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑ccccccccag                                                                                                                         >  1:803607/1‑136 (MQ=255) ttAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑ccccccccag                                                                                                                         >  2:422457/1‑136 (MQ=255)   aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑ccccccccagca                                                                                                                       >  1:941976/1‑134 (MQ=255)   aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑ccccccccagca                                                                                                                       >  2:3005/1‑134 (MQ=255)   aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑ccccccccagca                                                                                                                       >  2:409615/1‑134 (MQ=255)   aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑ccccccccagca                                                                                                                       >  1:411734/1‑134 (MQ=255)    gCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCaa                                                                                                                       >  2:1165320/1‑139 (MQ=255)        aGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTc                                                                                                                   >  2:971629/1‑139 (MQ=255)         gtgtGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTcc                                                                                                                  <  2:411734/139‑1 (MQ=255)                             aaaaaTAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATc                                                                                              >  1:339439/1‑139 (MQ=255)                                                     tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACCCCCCAACATCCCGCGTACAATCTGATTCGTCACCTGGCGAacggcg                                                                      >  1:297010/1‑139 (MQ=39)                                                     tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑cccc                                                                                                                               <  1:844297/81‑1 (MQ=255)                                                     tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑cccc                                                                                                                               >  2:844297/1‑81 (MQ=255)                                                       cAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCt                                                                    <  1:971629/139‑1 (MQ=255)                                                          tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACt                                                                                                                    >  2:696167/1‑88 (MQ=255)                                                          tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACt                                                                                                                    <  1:696167/88‑1 (MQ=255)                                                          tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGt                                                                                      >  1:372149/1‑118 (MQ=255)                                                          tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGt                                                                                      <  2:372149/118‑1 (MQ=255)                                                            gggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTgg                                                               >  2:732201/1‑139 (MQ=255)                                                              gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgc                                                                          <  2:662938/126‑1 (MQ=255)                                                              gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAgc                                                                          >  1:662938/1‑126 (MQ=255)                                                              gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGcc                                                             <  2:297010/139‑1 (MQ=255)                                                                 ggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAAc                                                                                                             >  2:325742/1‑88 (MQ=255)                                                                 ggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAAc                                                                                                             <  1:325742/88‑1 (MQ=255)                                                                 ggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCAtt                                                          >  1:977654/1‑139 (MQ=255)                                                                    aTGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt                                                       <  2:339439/139‑1 (MQ=255)                                                                      gATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACt                                                                                                                    <  2:887337/76‑1 (MQ=255)                                                                      gATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACt                                                                                                                    >  1:887337/1‑76 (MQ=255)                                                                       aTAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGc                                                    >  2:395066/1‑139 (MQ=255)                                                                               cgcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCt                                                                    >  1:5656/1‑115 (MQ=255)                                                                               cgcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCt                                                                    <  2:5656/115‑1 (MQ=255)                                                                                gcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCaaa                                                                                                                      <  1:513513/64‑1 (MQ=255)                                                                                gcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCaaa                                                                                                                      >  2:513513/1‑64 (MQ=255)                                                                                 ccAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCAt                                          >  2:369592/1‑139 (MQ=255)                                                                                    gCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcaccacac                                             >  1:370672/1‑133 (MQ=255)                                                                                    gCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcaccacac                                             <  2:370672/133‑1 (MQ=255)                                                                                                gCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑cc                                                                                                                                 <  2:743115/36‑1 (MQ=255)                                                                                                gCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑cc                                                                                                                                 >  1:743115/1‑36 (MQ=255)                                                                                                 cAAAGACAAATCTTGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCcgcg                          <  2:572036/139‑1 (MQ=255)                                                                                                     gACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGt                      <  1:409615/139‑1 (MQ=255)                                                                                                          aTCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGt                 >  1:975386/1‑139 (MQ=255)                                                                                                          aTCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGt                 >  1:674700/1‑139 (MQ=255)                                                                                                          aTCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGt                 >  1:1101831/1‑139 (MQ=255)                                                                                                           tctcGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTa                <  1:422457/139‑1 (MQ=255)                                                                                                               gCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCcgcg                                                                                                     >  2:284230/1‑50 (MQ=255)                                                                                                               gCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCcgcg                                                                                                     <  1:284230/50‑1 (MQ=255)                                                                                                                ccTGACTTACCTTCTTCTCCCCCCCCAGCAAACTcc                                                                                                                  >  2:588518/1‑36 (MQ=255)                                                                                                                ccTGACTTACCTTCTTCTCCCCCCCCAGCAAACTcc                                                                                                                  <  1:588518/36‑1 (MQ=255)                                                                                                                       tACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACaaa    <  1:1165320/139‑1 (MQ=255)                                                                                                                         ccTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAat  <  1:395066/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                  TTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCT‑CCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACAAAAT  >  NC_000913/4488444‑4488701

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