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%	141.1 / NA	34	coding (186/921 nt)	yedA	amino acid exporter for phenylalanine, threonine
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (17/17);  total (17/17)

GGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029597‑2029858                                                                                                                                |                                                                                                                                       ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:430324/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:321633/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:219485/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:991548/1‑139 (MQ=255) ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTAGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:769624/1‑139 (MQ=255)         aCCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCgctgct                                                                                                                    >  2:994290/1‑139 (MQ=255)                   tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                                                                       <  2:710660/111‑1 (MQ=255)                   tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                                                                       >  1:710660/1‑111 (MQ=255)                   tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgcg                                                                                                          >  1:224170/1‑139 (MQ=255)                     aTTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCt                                                                                                        >  2:1064976/1‑139 (MQ=255)                       tCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGa                                                                                                      <  1:867737/139‑1 (MQ=255)                                    gAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTAtt                                                                                         >  2:50176/1‑139 (MQ=255)                                     aaaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTg                                                                                        <  1:769624/139‑1 (MQ=255)                                         cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    <  1:850327/139‑1 (MQ=255)                                            gCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctggctg                                                                                 <  1:994290/139‑1 (MQ=255)                                             ccTCCGTTAATGATGGCGGGCGTTCGCTTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAATCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                <  2:426714/139‑1 (MQ=255)                                                 cGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGt                                                                            >  2:674154/1‑139 (MQ=255)                                                       tgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGc                                                                      <  2:224170/139‑1 (MQ=255)                                                           ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                                                                  >  2:629460/1‑139 (MQ=255)                                                                       aTTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCg                                                      <  1:629460/139‑1 (MQ=255)                                                                       aTTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCg                                                      <  1:1064976/139‑1 (MQ=255)                                                                        ttCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGa                                                     <  1:674154/139‑1 (MQ=255)                                                                                  ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGc                                                                                                 >  2:364868/1‑85 (MQ=255)                                                                                  ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGc                                                                                                 <  1:364868/85‑1 (MQ=255)                                                                                   cGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATg                                          <  1:219485/139‑1 (MQ=255)                                                                                                cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  <  1:1023659/70‑1 (MQ=255)                                                                                                cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  >  2:1023659/1‑70 (MQ=255)                                                                                                cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  1:369274/1‑139 (MQ=255)                                                                                                cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  2:832435/1‑139 (MQ=255)                                                                                                   gCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcgc                          <  2:369274/139‑1 (MQ=255)                                                                                                                  cgcgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGcc                                                       <  1:490249/95‑1 (MQ=255)                                                                                                                  cgcgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGcc                                                       >  2:490249/1‑95 (MQ=255)                                                                                                                    cgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCaa         <  2:611104/139‑1 (MQ=255)                                                                                                                    cgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCaa         >  1:611104/1‑139 (MQ=255)                                                                                                                          caAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGc   <  1:50176/139‑1 (MQ=255)                                                                                                                           aaaCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  2:991548/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 < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑