Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,586,594	(G)6→5	coding (168/3513 nt)	hsdR ←	endonuclease R Type I restriction enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,586,589	0	G	.	100.0%	199.9 / NA	42	coding (173/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/22);  total (20/22)

GGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAG  >  NC_000913/4586440‑4586734                                                                                                                                                      |                                                                                                                                                  ggCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgg                                                                                                                                                  >  1:632506/1‑151 (MQ=255) ggCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAg                                                                                                                                                   <  2:396911/150‑1 (MQ=255)   cacaTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggg                                                                                                                                                <  2:108055/151‑1 (MQ=255)     caTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggg                                                                                                                                               >  2:291863/1‑150 (MQ=255)      aTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggggat                                                                                                                                             >  1:567861/1‑149 (MQ=255)         tGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGtt                                                                                                                                         >  2:56263/1‑151 (MQ=255)           ggCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGa                                                                                                                                       <  2:72015/151‑1 (MQ=255)                   tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAAc                                                                                                                                  >  2:222359/1‑148 (MQ=255)                   tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAAc                                                                                                                                  >  2:135521/1‑148 (MQ=255)                   tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAAc                                                                                                                                  >  1:119925/1‑148 (MQ=255)                         cGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAg                                                                                                                          <  1:844906/150‑1 (MQ=255)                               gATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGttgc                                                                                                                   <  1:279921/151‑3 (MQ=255)                                     cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCtgtg                                                                                                             >  2:221286/1‑151 (MQ=255)                                        gCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGc                                                                                                           <  2:135563/150‑1 (MQ=255)                                          aCCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCtt                                                                                                         >  2:194432/1‑150 (MQ=255)                                             gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGcc                                                                                                     >  1:463010/1‑151 (MQ=255)                                             gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGcc                                                                                                     >  2:554801/1‑151 (MQ=255)                                                 ggTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaa                                                                                                   <  1:520224/149‑1 (MQ=255)                                                    tACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACAt                                                                                               >  2:484819/1‑150 (MQ=255)                                                      ccAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACg                                                                                            <  1:606256/151‑1 (MQ=255)                                                        aaTGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCa                                                                                          <  2:632546/151‑1 (MQ=255)                                                              gcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTaa                                                                                    <  1:739900/151‑1 (MQ=255)                                                              gcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTaa                                                                                    <  1:316452/151‑1 (MQ=255)                                                                gTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATc                                                                                  <  1:777744/151‑1 (MQ=255)                                                                  aaTTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCa                                                                                 <  1:829295/150‑1 (MQ=255)                                                                       gtgAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtc                                                                            <  2:593612/150‑1 (MQ=255)                                                                            aTACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTg                                                                      >  1:656859/1‑151 (MQ=255)                                                                               cAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGa                                                                   <  2:732623/151‑1 (MQ=255)                                                                                 gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGat                                                                  >  2:104638/1‑150 (MQ=255)                                                                                  agagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGTAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatc                                                                 >  2:707820/1‑150 (MQ=255)                                                                                         gttgtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGt                                                         >  2:414176/1‑151 (MQ=255)                                                                                           tgtCATCAACAAAGGCGATTTTGCCGATTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGAGGTTGTGTAACATACCAAGATGTTTCGCTGTGGCTTGGCCAAACATACGCTTTTTAATCAGCGTCGTGTTGGGATCATCCGGGt                                                         >  1:421887/1‑149 (MQ=255)                                                                                             tcatcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGtt                                                     >  1:21115/1‑151 (MQ=255)                                                                                               atcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAt                                                   <  1:554763/151‑1 (MQ=255)                                                                                                   acaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTc                                               >  2:495125/1‑151 (MQ=255)                                                                                                    caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTcc                                              <  2:15114/151‑1 (MQ=255)                                                                                                      aaGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTccgc                                            >  1:163554/1‑151 (MQ=255)                                                                                                            aTTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGcac                                       >  1:444434/1‑150 (MQ=255)                                                                                                                  ccGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGAt                                <  2:187706/151‑1 (MQ=255)                                                                                                                  ccGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGAt                                <  2:724427/151‑1 (MQ=255)                                                                                                                   cGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATg                               <  1:709144/151‑1 (MQ=255)                                                                                                                               ggAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAg                   <  1:222335/151‑1 (MQ=255)                                                                                                                                 agaTCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTc                 <  2:599584/151‑1 (MQ=255)                                                                                                                                    tCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGtt              >  1:755525/1‑151 (MQ=255)                                                                                                                                     cGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGtt              <  1:28879/150‑1 (MQ=255)                                                                                                                                          tGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGcc        >  2:579286/1‑151 (MQ=255)                                                                                                                                                 cacaA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAg  <  1:527276/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAG  >  NC_000913/4586440‑4586734

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