Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	602,835	G→A	A293T (GCT→ACT)	pheP →	phenylalanine:H(+) symporter PheP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	602,835	0	G	A	100.0%	66.2 / NA	23	A293T (GCT→ACT)	pheP	phenylalanine:H(+) symporter PheP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/7);  total (16/7)

TTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATG  >  NC_000913/602593‑603071                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             ttttGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTgcg                                                                                                                                                                                                                                      >  1:158992/1‑251 (MQ=255)                                       cgggAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAATCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGATTTACTGTCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCTCTTTTGAGGATTTTCCAGATTCTCGACAGCAATGTGGTAACTTCTGCGCTGAACCTCGTCATTCTGGTAGCATCGCTGCCAgtgt                                                                                                                                                                                                 <  1:237703/248‑1 (MQ=255)                                       cTGGAGCGGCTTGGGCTTACTGCCGCTGAAGCGCGCGATCCGGAACACAGCATTCCAAACGCGGTAAATCCGGTGGTGTATCGAATCCTGCTGTTTTACATCGGTCCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGCAGTGAACTCCAACAGTAGCCAGTTTGTGATGCTTTTCCATAATCTCGTCAGCAACGTGGTAACTTCTGCGCTGAACTTCGCCATTCTGGTAGCATCGCTGTCAgtgt                                                                                                                                                                                                 >  2:177831/1‑250 (MQ=255)                                                   gggATTACTGCCGCTGAAGCGCGTGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTTGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCTACAGCAACGTGGTAACTTCTGCGCTTAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGt                                                                                                                                                                                    <  1:260149/251‑1 (MQ=255)                                                                         gcgATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTg                                                                                                                                                              <  1:65544/251‑1 (MQ=255)                                                                              ccGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTgg                                                                                                                                                         >  2:198741/1‑251 (MQ=255)                                                                                    aaaaGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTc                                                                                                                                                   <  1:238535/251‑1 (MQ=255)                                                                                                      gTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGcc                                                                                                                                 >  1:21715/1‑251 (MQ=255)                                                                                                                                 ctgTTTTACATCGGTTAACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACCGCAACGTGGTAACTTCTGCGCTGAACATCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCAAGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAGGTTTTTGACTAGCGTCAGCCGTcg                                                                                                      >  1:34978/1‑251 (MQ=255)                                                                                                                                         cATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGc                                                                                              >  1:176113/1‑251 (MQ=255)                                                                                                                                                                         aTCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCGGAACTTCGTCATTCGGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTTGCCTTGCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGGCAGCCGTCGCGGTGTTCCGAGTAACTCTCTGTGGCTTTCCGTAGCGAt                                                               <  2:172003/250‑1 (MQ=255)                                                                                                                                                                              tgggtggAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTc                                                         >  1:134030/1‑251 (MQ=255)                                                                                                                                                                                ggtggAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGc                                                       >  1:19700/1‑251 (MQ=255)                                                                                                                                                                                ggtggAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGc                                                       >  2:87046/1‑251 (MQ=255)                                                                                                                                                                                          acaTCCAACAGTAGCCCGTTTGTGATGATTTTCCAACCTCTCGCCAGCAACGCGGTAACTTCTGCGCTCAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATCACCGCGGGGTTTACTCTAACAGCCGCAAGCTGTCTAGCCTTTCTGTGCAGGGTACTGCGCCGACGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGCTGCTTTCAGGAGCGATCACTTCGCTGGTGGTGtt                                             >  2:137258/3‑251 (MQ=255)                                                                                                                                                                                                    gTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTTGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGGGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTAt                                   <  2:158992/251‑1 (MQ=255)                                                                                                                                                                                                      aGCCCGTTTGTGATGATTTTCCATAATCTCGACCGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGCCTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCACCTATct                                 >  2:247098/1‑251 (MQ=255)                                                                                                                                                                                                       gCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATctg                                >  1:158356/1‑251 (MQ=255)                                                                                                                                                                                                        cccGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATctgc                               >  1:49538/1‑251 (MQ=255)                                                                                                                                                                                                         ccGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATctgct                              <  1:198741/251‑1 (MQ=255)                                                                                                                                                                                                                tgatgaTTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCa                        >  2:65339/1‑250 (MQ=255)                                                                                                                                                                                                                              aaTCTCGACAGAAGAGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTATCAGTGTATCACAGCGGGGTTTACTCTAAAAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGAAATGTGCAGAAGTTTTTGACTCGCGTCAGCCGACGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGGGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTct         >  2:40887/1‑251 (MQ=255)                                                                                                                                                                                                                                     aCAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGGGCAGGGTAATGCGCCGCAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCACAAAGCGTTTGGTCTGCTGATg  >  2:11107/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             TTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATG  >  NC_000913/602593‑603071

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