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%	200.0 / NA	51	D237N (GAT→AAT)	pheP	phenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/21);  total (30/21)

CGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCC  >  NC_000913/602522‑602811                                                                                                                                                  |                                                                                                                                                 cggcgAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCg                                                                                                                                             <  1:739725/151‑1 (MQ=255)      agaAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaaaa                                                                                                                                        <  2:353417/151‑1 (MQ=255)      agaAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaaaa                                                                                                                                        <  1:666041/151‑1 (MQ=255)      agaAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTGCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaa                                                                                                                                          <  2:72912/149‑1 (MQ=255)       gaAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaaaa                                                                                                                                        >  1:291720/1‑150 (MQ=255)        aaaGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaaaaa                                                                                                                                       >  2:351760/1‑150 (MQ=255)         aaGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGaaaaa                                                                                                                                       >  2:311482/1‑149 (MQ=255)          aGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGc                                                                                                                                    <  1:107876/151‑1 (MQ=255)           gCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCCGAAGCGCGCAATCCGGAAAAAAGc                                                                                                                                    >  2:303206/1‑150 (MQ=255)               gTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTcc                                                                                                                               <  2:748103/151‑1 (MQ=255)                    gACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAg                                                                                                                          <  2:722327/151‑1 (MQ=255)                     aCAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCaaaa                                                                                                                           >  2:727501/1‑149 (MQ=255)                     aCAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGc                                                                                                                         <  1:150822/151‑1 (MQ=255)                        aCCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCgg                                                                                                                       >  2:737321/1‑150 (MQ=255)                           tctGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTaa                                                                                                                    >  1:426635/1‑150 (MQ=255)                                 gcTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAggt                                                                                                             <  1:167653/151‑1 (MQ=255)                                    aCGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGGAAATCAggtgg                                                                                                           >  2:75586/1‑150 (MQ=255)                                             tcttcGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCa                                                                                                  <  2:407410/150‑1 (MQ=255)                                                 cGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCctg                                                                                             >  1:241170/1‑151 (MQ=255)                                                        ggCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGtttt                                                                                      >  2:284956/1‑151 (MQ=255)                                                        ggCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGtt                                                                                        <  1:346696/149‑1 (MQ=255)                                                          cTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGGGGTGTATCGCATCCTGCTGTTTTa                                                                                     >  2:47860/1‑150 (MQ=255)                                                           tGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCctgct                                                                                           >  2:486815/1‑143 (MQ=255)                                                           tGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCctgct                                                                                           <  1:486838/143‑1 (MQ=255)                                                                tGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGt                                                                              >  1:512140/1‑151 (MQ=255)                                                                 gggCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGtt                                                                             <  2:484695/151‑1 (MQ=255)                                                                 gggCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGtt                                                                             <  1:169866/151‑1 (MQ=255)                                                                  ggCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTc                                                                            <  2:89321/151‑1 (MQ=255)                                                                  ggCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTc                                                                            >  1:243467/1‑151 (MQ=255)                                                                  ggCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTc                                                                            >  1:746304/1‑151 (MQ=255)                                                                        ttttGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACtgg                                                                       >  1:704148/1‑150 (MQ=255)                                                                                cTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTAc                                                              >  1:28984/1‑151 (MQ=255)                                                                                  ggCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTg                                                            <  2:109523/151‑1 (MQ=255)                                                                                     ggTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGgcg                                                         >  1:397736/1‑151 (MQ=255)                                                                                        aaTTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGctc                                                      >  1:47961/1‑151 (MQ=255)                                                                                         aTTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGctct                                                     >  1:698605/1‑151 (MQ=255)                                                                                                   ccTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCgtgggt                                           <  1:376228/151‑1 (MQ=255)                                                                                                       cggcggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGa                                        <  1:195384/150‑1 (MQ=255)                                                                                                        ggcggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAg                                      <  2:631000/151‑1 (MQ=255)                                                                                                           ggTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTg                                    >  2:428051/1‑150 (MQ=255)                                                                                                             tCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGaa                                  >  2:11864/1‑150 (MQ=255)                                                                                                              cTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAAt                                >  1:708737/1‑151 (MQ=255)                                                                                                                               tACTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGttt               >  1:411516/1‑151 (MQ=255)                                                                                                                                aCTGCCGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGttt               <  2:383227/150‑1 (MQ=255)                                                                                                                                     cGCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGtgatg         >  2:327112/1‑151 (MQ=255)                                                                                                                                      gCTGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGGGGAAGTGAAATCCAACAGTAGCCCGTTTGtgat          >  1:422406/1‑149 (MQ=255)                                                                                                                                        tGAAGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGAtt      >  2:218295/1‑151 (MQ=255)                                                                                                                                          aaGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGAtttt    >  2:173838/1‑151 (MQ=255)                                                                                                                                          aaGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGCGGTGTATCGCATCCTGCTGTTTTACATCGGATCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGAtttt    >  2:606787/1‑151 (MQ=255)                                                                                                                                           aGCGCGCAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTc   <  1:543862/151‑1 (MQ=255)                                                                                                                                            gcgcgcAATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTcc  >  2:701900/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGGCGAGAAAGCCAGTATCGACAACCTCTGGCGCTACGGTGGTTTCTTCGCCACCGGCTGGAATGGGCTGATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCC  >  NC_000913/602522‑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: ‑