Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,320,065	T→C	D200G (GAC→GGC)	phnI ←	ribophosphonate triphosphate synthase complex putative catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,320,065	0	T	C	100.0%	172.0 / NA	48	D200G (GAC→GGC)	phnI	ribophosphonate triphosphate synthase complex putative catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/17);  total (31/17)

CAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGA  >  NC_000913/4319917‑4320210                                                                                                                                                     |                                                                                                                                                  cagcagTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTtcgccg                                                                                                                                                 <  1:468234/151‑1 (MQ=255)     agTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTcgccgccgc                                                                                                                                             >  1:223934/1‑151 (MQ=255)        tCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCgcgc                                                                                                                                           >  2:570251/1‑150 (MQ=255)           ccGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCa                                                                                                                                       <  1:101606/151‑1 (MQ=255)             gACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAAc                                                                                                                                     <  1:401544/151‑1 (MQ=255)                 ttCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAActgct                                                                                                                                 >  1:59084/1‑151 (MQ=255)                   cACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACtgctgc                                                                                                                               >  1:457493/1‑151 (MQ=255)                     ccGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCaa                                                                                                                             <  1:303450/151‑1 (MQ=255)                     ccGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCa                                                                                                                              <  2:640589/150‑1 (MQ=255)                     ccGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCa                                                                                                                              >  1:34307/1‑150 (MQ=255)                       gCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAAc                                                                                                                           <  2:80662/151‑1 (MQ=255)                             cccAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGggagg                                                                                                                     >  1:654365/1‑151 (MQ=255)                                aGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGc                                                                                                                  <  2:609016/151‑1 (MQ=255)                                 gCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCg                                                                                                                 >  2:639008/1‑151 (MQ=255)                                 gCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCg                                                                                                                 >  1:432775/1‑151 (MQ=255)                                  ctctTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCg                                                                                                                 >  2:599771/1‑150 (MQ=255)                                    ctTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGt                                                                                                                <  2:560025/149‑1 (MQ=255)                                     ttCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTg                                                                                                               >  1:124928/1‑149 (MQ=255)                                       ccGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGc                                                                                                            <  1:67043/150‑1 (MQ=255)                                            acaATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGt                                                                                                      >  2:310705/1‑151 (MQ=255)                                                     acGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcggcg                                                                                             >  1:470481/1‑151 (MQ=255)                                                     acGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcggcg                                                                                             >  1:18440/1‑151 (MQ=255)                                                      cGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcggcgg                                                                                            >  1:397785/1‑151 (MQ=255)                                                      cGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcggcgg                                                                                            >  2:370923/1‑151 (MQ=255)                                                          gATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTg                                                                                        >  1:402179/1‑151 (MQ=255)                                                            tataGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCg                                                                                      >  2:179083/1‑151 (MQ=255)                                                             ataGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCgg                                                                                     <  1:35025/151‑1 (MQ=255)                                                               aGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGt                                                                                   <  1:144529/151‑1 (MQ=255)                                                                gCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTg                                                                                  >  1:323634/1‑151 (MQ=255)                                                                 ccGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTg                                                                                  >  2:450811/1‑150 (MQ=255)                                                                    cTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGt                                                                              >  1:230103/1‑151 (MQ=255)                                                                      gCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTc                                                                             >  2:46996/1‑150 (MQ=255)                                                                       cGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCAt                                                                           >  1:570420/1‑151 (MQ=255)                                                                                 cTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTgc                                                                 >  1:96085/1‑151 (MQ=255)                                                                                 cTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTgc                                                                 >  2:529395/1‑151 (MQ=255)                                                                                    cGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGcgcc                                                              >  2:312443/1‑151 (MQ=255)                                                                                     gAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGcgccg                                                             >  1:224795/1‑151 (MQ=255)                                                                                                   gcCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCaaa                                               <  2:159420/151‑1 (MQ=255)                                                                                                     ccGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACt                                             >  1:573785/1‑151 (MQ=255)                                                                                                            cgcgTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAg                                      <  1:70031/151‑1 (MQ=255)                                                                                                              cgTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGcc                                    <  2:661824/151‑1 (MQ=255)                                                                                                               gTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGccc                                   <  2:59082/151‑1 (MQ=255)                                                                                                                ttGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGcccc                                  >  1:92222/1‑151 (MQ=255)                                                                                                                      ggAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGAcgc                            <  1:55611/151‑1 (MQ=255)                                                                                                                       gAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGAcgcg                           >  1:519772/1‑151 (MQ=255)                                                                                                                            ggCCAGCGCCAGCAAATAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCcagc                      <  1:318104/151‑1 (MQ=255)                                                                                                                                         aaaTAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACg          >  2:392817/1‑150 (MQ=255)                                                                                                                                         aaaTAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGt         >  1:571600/1‑151 (MQ=255)                                                                                                                                          aaTAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTg        >  1:316438/1‑151 (MQ=255)                                                                                                                                          aaTAGCCTTCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTg        <  2:173802/151‑1 (MQ=255)                                                                                                                                                  tcgccgCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGa  <  1:154665/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGA  >  NC_000913/4319917‑4320210

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