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%	258.3 / NA	54	coding (173/3513 nt)	hsdR	endonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (25/29);  total (25/29)

GGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGG  >  NC_000913/4586449‑4586735                                                                                                                                             |                                                                                                                                                   gggCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggggatg                                                                                                                                             <  2:800154/148‑4 (MQ=255)     aTCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGt                                                                                                                                      >  1:831899/1‑150 (MQ=255)     aTCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGt                                                                                                                                      >  1:636501/1‑150 (MQ=255)       cGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTaa                                                                                                                                    <  2:590087/150‑1 (MQ=255)          tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAg                                                                                                                                 <  2:26510/150‑1 (MQ=255)          tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAg                                                                                                                                 <  1:767785/150‑1 (MQ=255)          tGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAg                                                                                                                                 <  1:443008/150‑1 (MQ=255)            agaTCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGAcc                                                                                                                              >  1:219108/1‑151 (MQ=255)                cGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGa                                                                                                                          <  1:262132/151‑1 (MQ=255)                  ttATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATg                                                                                                                        >  1:121995/1‑151 (MQ=255)                  ttATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATg                                                                                                                        >  1:354567/1‑151 (MQ=255)                        tatTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                                  >  1:798676/1‑151 (MQ=255)                         gttCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGGGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGAAACAGACCAAGATGTATAGCt                                                                                                                 >  2:103024/2‑151 (MQ=255)                                cACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCtt                                                                                                          <  2:9331/151‑1 (MQ=255)                                 aCCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTc                                                                                                         <  1:1172/151‑1 (MQ=255)                                 aCCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTc                                                                                                         <  1:631579/151‑1 (MQ=255)                                   cGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGc                                                                                                       >  1:464709/1‑151 (MQ=255)                                     ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCa                                                                                                     <  1:204921/151‑1 (MQ=255)                                        ggTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaaa                                                                                                   <  2:563859/150‑1 (MQ=255)                                         gTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACa                                                                                                 <  2:630966/151‑1 (MQ=255)                                           tACCAATGCGGCGTAATTTTTGAAATACAGAGAGGATGTTGGCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaa                                                                                                    >  2:346266/1‑146 (MQ=255)                                           tACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACAt                                                                                                >  2:689644/1‑150 (MQ=255)                                             ccAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACg                                                                                             >  2:325069/1‑151 (MQ=255)                                                  gcggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAtt                                                                                         <  2:202115/150‑1 (MQ=255)                                                   cggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAtttt                                                                                       >  1:698993/1‑151 (MQ=255)                                                     gcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTaa                                                                                     <  2:15/151‑1 (MQ=255)                                                       gTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATc                                                                                   <  2:495150/151‑1 (MQ=255)                                                          aTTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGc                                                                                <  1:545748/151‑1 (MQ=255)                                                                 aaaTACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCgtgt                                                                         >  1:801084/1‑151 (MQ=255)                                                                          gagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatcat                                                                >  2:112081/1‑151 (MQ=255)                                                                           agGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatcatc                                                               >  2:763105/1‑151 (MQ=255)                                                                              aTGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCgg                                                            <  2:244103/151‑1 (MQ=255)                                                                                gttgtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGt                                                          >  1:242948/1‑151 (MQ=255)                                                                                  tgtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGt                                                          >  1:463246/1‑149 (MQ=255)                                                                                     catcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGtt                                                      >  1:436103/1‑150 (MQ=255)                                                                                     catcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTa                                                     <  1:420614/151‑1 (MQ=255)                                                                                        caacaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAttt                                                  <  1:804245/151‑1 (MQ=255)                                                                                        caacaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAtt                                                   <  1:132481/150‑1 (MQ=255)                                                                                                 cGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCcgccgc                                          >  2:871620/1‑150 (MQ=255)                                                                                                  gATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGcac                                        <  2:831832/151‑1 (MQ=255)                                                                                                   aTTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGcac                                        >  1:441154/1‑150 (MQ=255)                                                                                                         ccGAGTTCACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGAt                                 <  1:112103/151‑1 (MQ=255)                                                                                                                cACGCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATa                           >  2:907141/1‑150 (MQ=255)                                                                                                                   gCAGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTg                       >  1:284631/1‑151 (MQ=255)                                                                                                                     aGGAGATCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGa                      >  2:462935/1‑150 (MQ=255)                                                                                                                        agaTCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTc                  <  1:852297/151‑1 (MQ=255)                                                                                                                        agaTCGTGTTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTc                  >  1:446525/1‑151 (MQ=255)                                                                                                                             gtgtTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTg              <  2:108176/150‑1 (MQ=255)                                                                                                                             gtgtTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTg              >  1:706924/1‑150 (MQ=255)                                                                                                                             gtgtTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGa             <  2:706875/151‑1 (MQ=255)                                                                                                                             gtgtTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGa             <  2:909797/151‑1 (MQ=255)                                                                                                                               gtTGATTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGAc            <  1:876612/150‑1 (MQ=255)                                                                                                                                   aTTCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCt       <  1:849486/151‑1 (MQ=255)                                                                                                                                    ttCTCACAA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGTGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGCTGGCATAAGTGAAGTCGTTGACGCCCt       <  2:287262/150‑1 (MQ=255)                                                                                                                                        cacaA‑GGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAgg  <  1:457763/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                   GGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGG  >  NC_000913/4586449‑4586735

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