Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,581,621	G→A	A476A (GCC→GCT)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,621	0	G	A	100.0%	107.9 / NA	34	A476A (GCC→GCT)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/21);  total (13/21)

TCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTT  >  NC_000913/4581487‑4581751                                                                                                                                       |                                                                                                                                   tcCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCAt                                                                                                                                <  1:896228/139‑1 (MQ=255)      cAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                           <  2:146010/139‑1 (MQ=255)          tGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGt                                                                                                                       <  2:748571/139‑1 (MQ=255)            aCGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCt                                                                                                                     >  1:806347/1‑139 (MQ=255)                    aaTCGGCCTCATCGCTCGCCCCCATTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCttt                                                                                                             <  2:318045/139‑1 (MQ=255)                        ggCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                         <  2:798017/139‑1 (MQ=255)                         gCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                        <  2:6220/139‑1 (MQ=255)                         gCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                        <  1:849078/139‑1 (MQ=255)                                                 cGCATCAGCGCATCCAGTTCAGACAGCGCCTTTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATAc                                                                                                                          <  2:651809/97‑1 (MQ=255)                                                 cGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATAc                                                                                                                          >  1:651809/1‑97 (MQ=255)                                                                   tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                           <  2:981631/78‑1 (MQ=255)                                                                   tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                           >  1:981631/1‑78 (MQ=255)                                                                   tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACt                                                              >  2:179254/1‑139 (MQ=255)                                                                       aCAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACg                                                          >  1:732563/1‑139 (MQ=255)                                                                           cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTg                                                      >  2:972101/1‑139 (MQ=255)                                                                                      gTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCa                                           <  1:972101/139‑1 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                           >  1:758760/1‑54 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                           <  2:758760/54‑1 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGt                                                                                                                       >  1:264529/1‑58 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGt                                                                                                                       <  2:264529/58‑1 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGc                                      >  2:957912/1‑139 (MQ=255)                                                                                              aTCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTgg                                   <  2:806347/139‑1 (MQ=255)                                                                                                   ttCTGCCTCTAATGCATCCGGCTCCGGCAGGCTGTCAGGATTAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTTCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCa                              <  1:330827/139‑1 (MQ=255)                                                                                                      tGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAg                                >  2:479323/1‑134 (MQ=255)                                                                                                      tGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAg                                <  1:479323/134‑1 (MQ=255)                                                                                                                  aTCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGgtg               <  1:894641/139‑1 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                             >  2:322712/1‑73 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                             <  1:322712/73‑1 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGAt                                                                   <  2:762113/83‑1 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATAACCAGCGAATCGGATTTTGCGGTGCGGAt                                                                   >  1:762113/1‑83 (MQ=255)                                                                                                                       gctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGttttt          <  1:652471/139‑1 (MQ=255)                                                                                                                              cAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTG‑TTTTGTTCTCtt  <  1:957912/139‑1 (MQ=255)                                                                                                                                 gCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTtc      <  2:195601/133‑1 (MQ=255)                                                                                                                                 gCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTtc      >  1:195601/1‑133 (MQ=255)                                                                                                                                       |                                                                                                                                   TCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTT  >  NC_000913/4581487‑4581751

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