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%	199.1 / NA	62	A476A (GCC→GCT)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/32);  total (30/32)

TCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCT  >  NC_000913/4581487‑4581755                                                                                                                                       |                                                                                                                                       tcCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCGTCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCAt                                                                                                                                    <  1:852649/139‑1 (MQ=255)      cAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                               <  1:232498/139‑1 (MQ=255)       agcgtcagttGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATAc                                                                                                                              >  1:1377057/10‑139 (MQ=255)            aCGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCt                                                                                                                         >  2:230168/1‑139 (MQ=255)                         gCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                            <  1:518673/139‑1 (MQ=255)                         gCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                            <  1:453024/139‑1 (MQ=255)                          ccTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCa                                                                                                           >  1:1409445/1‑139 (MQ=255)                              aTCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGgag                                                                                                       >  1:963445/1‑139 (MQ=255)                                          cAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCttt                                                                                                                 >  2:380651/1‑117 (MQ=255)                                          cAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCttt                                                                                                                 <  1:380651/117‑1 (MQ=255)                                                    aTCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                                 <  1:230168/139‑1 (MQ=255)                                                       aGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCAt                                                                                                                                    >  1:1119561/1‑85 (MQ=255)                                                       aGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                              >  2:197636/1‑139 (MQ=255)                                                       aGCGCATCCAGTTCAGACAGCGCATGTACCAGTTCACCCATCGCTTCTGCTGCTAATACATCCGGCTCCGGCAGGCTGTCAGCAt                                                                                                                                    <  2:1119561/85‑1 (MQ=255)                                                         cgcATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                                 <  2:336344/134‑1 (MQ=255)                                                         cgcATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                                 >  1:336344/1‑134 (MQ=255)                                                           cATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGCATCGGATTTtgcggtgcg                                                                          <  2:963445/139‑1 (MQ=255)                                                            aTCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                                 <  2:649371/131‑1 (MQ=255)                                                            aTCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                                 >  1:649371/1‑131 (MQ=255)                                                              ccAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCtt                                                                                                                  >  2:544170/1‑96 (MQ=255)                                                              ccAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCtt                                                                                                                  <  1:544170/96‑1 (MQ=255)                                                              ccAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGAt                                                                       <  2:1058665/139‑1 (MQ=255)                                                              ccAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGAt                                                                       >  1:551864/1‑139 (MQ=255)                                                               cAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATc                                                                      <  1:197636/139‑1 (MQ=255)                                                                  ttCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                            <  2:862207/98‑1 (MQ=255)                                                                  ttCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                                            >  1:862207/1‑98 (MQ=255)                                                                           cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAg                                                                                                          >  2:1695509/1‑91 (MQ=255)                                                                           cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAg                                                                                                          <  1:1695509/91‑1 (MQ=255)                                                                           cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTg                                                          >  1:252377/1‑139 (MQ=255)                                                                            gcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGa                                                         >  1:76303/1‑139 (MQ=255)                                                                            gcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGa                                                         >  2:625132/1‑139 (MQ=255)                                                                                 tACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgc                                                                                <  1:601230/111‑1 (MQ=255)                                                                                 tACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgc                                                                                >  2:601230/1‑111 (MQ=255)                                                                                 tACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTg                                                    >  1:1629800/1‑139 (MQ=255)                                                                                      gTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcggtg                                                                            >  1:143256/1‑110 (MQ=255)                                                                                      gTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcggtg                                                                            <  2:143256/110‑1 (MQ=255)                                                                                        tCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTTTCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGAt                                                                       >  2:117444/1‑113 (MQ=255)                                                                                        tCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGAt                                                                       <  1:117444/113‑1 (MQ=255)                                                                                         cGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCg                                            <  2:1629800/139‑1 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATAc                                                                                                                              >  1:266954/1‑55 (MQ=255)                                                                                           cccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATAc                                                                                                                              <  2:266954/55‑1 (MQ=255)                                                                                            ccATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCt                                         <  2:1409445/139‑1 (MQ=255)                                                                                              aTCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGTGGCTgg                                       <  1:1488624/139‑1 (MQ=255)                                                                                                          gCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGc                                             >  2:772281/1‑121 (MQ=255)                                                                                                          gCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGc                                             <  1:772281/121‑1 (MQ=255)                                                                                                          gCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCt                           >  1:1326200/1‑139 (MQ=255)                                                                                                             aaTACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGt                        <  1:625132/139‑1 (MQ=255)                                                                                                                    ccggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCg                      >  1:941334/1‑134 (MQ=255)                                                                                                                    ccggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCg                      <  2:941334/134‑1 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAGTCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCg                      <  2:1125078/132‑1 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCg                      >  1:1125078/1‑132 (MQ=255)                                                                                                                      ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCAGTGAGGATCCACTC‑CGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGttttt              >  2:915620/1‑139 (MQ=255)                                                                                                                       gctccggcAGGCTGTCAGCATCAATACTGTCTTTATCtt                                                                                                                  >  2:483669/1‑39 (MQ=255)                                                                                                                       gctccggcAGGCTGTCAGCATCAATACTGTCTTTATCtt                                                                                                                  <  1:483669/39‑1 (MQ=255)                                                                                                                        ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                              <  2:401036/74‑1 (MQ=255)                                                                                                                        ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                              >  1:401036/1‑74 (MQ=255)                                                                                                                        ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGc                                          <  1:1205017/110‑1 (MQ=255)                                                                                                                        ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGc                                          >  2:1205017/1‑110 (MQ=255)                                                                                                                           cggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATg                             >  1:1643052/1‑120 (MQ=255)                                                                                                                           cggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATg                             <  2:1643052/120‑1 (MQ=255)                                                                                                                                 gCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCg    <  2:252377/139‑1 (MQ=255)                                                                                                                                   tGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCt  <  2:76303/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       TCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCT  >  NC_000913/4581487‑4581755

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: ‑