Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	602,667	G→A	D237N (GAT→AAT)	pheP →	phenylalanine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	602,667	0	G	A	100.0%	201.6 / NA	52	D237N (GAT→AAT)	pheP	phenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/24);  total (28/24)

CGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCC  >  NC_000913/602525‑602811                                                                                                                                               |                                                                                                                                                 cgAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaa                                                                                                                                          <  2:50208/151‑1 (MQ=255)          cAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCAtt                                                                                                                                 <  2:670556/151‑1 (MQ=255)           aGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGGGCTGATTGGGGTTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCCtt                                                                                                                                 >  2:951711/1‑150 (MQ=255)               tCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCaa                                                                                                                             <  2:802272/150‑1 (MQ=255)                 gACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAg                                                                                                                          <  1:410825/151‑1 (MQ=255)                  aCAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAg                                                                                                                          >  1:995953/1‑150 (MQ=255)                  aCAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGc                                                                                                                         >  2:244032/1‑151 (MQ=255)                    aaCCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCgg                                                                                                                       >  1:359912/1‑151 (MQ=255)                    aaCCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCgg                                                                                                                       <  2:464129/151‑1 (MQ=255)                    aaCCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCgg                                                                                                                       >  1:603641/1‑151 (MQ=255)                      ccTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGt                                                                                                                      <  2:802444/150‑1 (MQ=255)                      ccTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGt                                                                                                                      >  1:339371/1‑150 (MQ=255)                                 aCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAggtgg                                                                                                           <  2:433005/150‑1 (MQ=255)                                  cGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGgtggtg                                                                                                         >  2:949482/1‑151 (MQ=255)                                    gtggtTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTa                                                                                                       <  1:572915/151‑1 (MQ=255)                                           cttcGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCAt                                                                                                 >  2:1052384/1‑150 (MQ=255)                                             tcGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCAt                                                                                                 <  1:455304/148‑1 (MQ=255)                                                ccACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCctgc                                                                                            >  2:967122/1‑150 (MQ=255)                                                ccACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCCGTAAATCAGGTGGCGTATCGCATCctg                                                                                             >  1:269734/1‑149 (MQ=255)                                                     ggCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGtttt                                                                                      <  1:151769/151‑1 (MQ=255)                                                      gCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTa                                                                                     <  2:804822/151‑1 (MQ=255)                                                        tGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTAc                                                                                    >  2:546163/1‑150 (MQ=255)                                                         ggAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACAt                                                                                  >  1:1036165/1‑151 (MQ=255)                                                         ggAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACAt                                                                                  >  1:998545/1‑151 (MQ=255)                                                           aaTGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCg                                                                                >  1:384534/1‑151 (MQ=255)                                                           aaTGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCg                                                                                <  2:931451/151‑1 (MQ=255)                                                            aTGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCg                                                                                <  2:542349/150‑1 (MQ=255)                                                              gggCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGt                                                                              <  1:217816/150‑1 (MQ=255)                                                                  tGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACt                                                                         >  2:1026124/1‑151 (MQ=255)                                                                   gATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtg                                                                        >  2:166665/1‑151 (MQ=255)                                                                     ttttGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtgg                                                                       >  1:1032921/1‑150 (MQ=255)                                                                       ttGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtggtgg                                                                    >  2:407058/1‑151 (MQ=255)                                                                       ttGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtggtg                                                                     <  2:395096/150‑1 (MQ=255)                                                                       ttGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtggtg                                                                     <  1:533528/150‑1 (MQ=255)                                                                          tCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGttt                                                                 >  2:880379/1‑151 (MQ=255)                                                                           cGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGtttt                                                                <  1:903988/151‑1 (MQ=255)                                                                           cGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGtttt                                                                >  2:459139/1‑151 (MQ=255)                                                                            gCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTa                                                               <  2:359901/151‑1 (MQ=255)                                                                               ggCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTg                                                            >  1:1015554/1‑151 (MQ=255)                                                                                gCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTgg                                                           <  1:63666/151‑1 (MQ=255)                                                                                   gTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGG‑AAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGct                                                       <  1:668351/151‑1 (MQ=255)                                                                                              ctcCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCgtgg                                             >  1:816942/1‑151 (MQ=255)                                                                                                ccTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCgtgg                                             <  2:157135/149‑1 (MQ=255)                                                                                                       cggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTg                                    <  2:558358/151‑1 (MQ=255)                                                                                                             ggAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATc                               <  2:40152/150‑1 (MQ=255)                                                                                                                gCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAAc                           >  2:713068/1‑151 (MQ=255)                                                                                                                    aTTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGt                        >  2:858246/1‑150 (MQ=255)                                                                                                                             aCTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTtg              >  1:263934/1‑151 (MQ=255)                                                                                                                                 ccGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGtgat          >  2:275521/1‑151 (MQ=255)                                                                                                                                 ccGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGAGGGCGGAAGTGAAATCCAACAGTAGCCCGTTTGtgat          >  1:585188/1‑151 (MQ=255)                                                                                                                                  cGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGtgatg         >  2:782702/1‑151 (MQ=255)                                                                                                                                   gCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGtgatga        <  2:413863/151‑1 (MQ=255)                                                                                                                                          cgcgcaATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTcc  >  2:744439/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                 CGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCC  >  NC_000913/602525‑602811

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