Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,029,729	(C)6→5	coding (191/921 nt)	yedA →	amino acid exporter for phenylalanine, threonine

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,029,724	0	C	.	100.0%	123.0 / NA	30	coding (186/921 nt)	yedA	amino acid exporter for phenylalanine, threonine
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (18/12);  total (18/12)

GGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029597‑2029858                                                                                                                                |                                                                                                                                       ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGTTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:439293/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:44850/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:490298/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:873749/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGACGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:535989/1‑139 (MQ=255)   ggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTcc                                                                                                                          >  1:226271/1‑139 (MQ=255)     cTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCgc                                                                                                                        >  2:519746/1‑139 (MQ=255)        aaCCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCgctgc                                                                                                                     >  2:197814/1‑139 (MQ=255)        aaCCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         <  1:671157/135‑1 (MQ=255)        aaCCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         >  2:671157/1‑135 (MQ=255)         aCCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCgctgct                                                                                                                    <  1:197814/139‑1 (MQ=255)                    cATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCgcgc                                                                                                         >  2:382693/1‑139 (MQ=255)                                      aaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgc                                                                                       <  2:406108/139‑1 (MQ=255)                                         cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    <  2:535989/139‑1 (MQ=255)                                                                                  ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAAt                                           >  2:788468/1‑139 (MQ=255)                                                                                                cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  1:845791/1‑139 (MQ=255)                                                                                                  ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTaa                                                                          <  2:781416/92‑1 (MQ=255)                                                                                                  ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTaa                                                                          >  1:781416/1‑92 (MQ=255)                                                                                                       ttttGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCgccgcc                      >  2:399600/1‑139 (MQ=255)                                                                                                          tGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTa                   <  1:382693/139‑1 (MQ=255)                                                                                                             tACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTg                >  2:257996/1‑139 (MQ=255)                                                                                                                 gcgcGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTg            >  1:661727/1‑139 (MQ=255)                                                                                                                  cgcgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                                                                             >  2:703560/1‑73 (MQ=255)                                                                                                                  cgcgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                                                                             <  1:703560/73‑1 (MQ=255)                                                                                                                    cgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCaa         >  2:549809/1‑139 (MQ=255)                                                                                                                         acaAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg    <  1:585958/139‑1 (MQ=255)                                                                                                                         acaAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg    <  2:226271/139‑1 (MQ=255)                                                                                                                         acaAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg    <  1:571123/139‑1 (MQ=255)                                                                                                                           aaaCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:490298/139‑1 (MQ=255)                                                                                                                           aaaCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  2:661727/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                       GGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029597‑2029858

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