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

ATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029588‑2029858                                                                                                                                         |                                                                                                                                       atatCATTTGAGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                                                                      >  1:354781/1‑139 (MQ=255)          ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:27938/1‑139 (MQ=255)          ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:571359/1‑139 (MQ=255)          ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:49959/1‑139 (MQ=255)             gCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         <  1:229499/139‑1 (MQ=255)                      aTTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAAt                                                                                                                >  1:664415/1‑139 (MQ=255)                                     ttGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTTGATTCCTGGCAGCCGGTATTTTATTGCTGGTATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGc                                                                                                 <  1:579368/139‑1 (MQ=255)                                               aaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgc                                                                                       <  1:27938/139‑1 (MQ=255)                                                      ccTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                <  1:757007/139‑1 (MQ=255)                                                                tgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGc                                                                      <  2:49959/139‑1 (MQ=255)                                                                        ggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                <  1:300709/121‑1 (MQ=255)                                                                        ggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                >  2:300709/1‑121 (MQ=255)                                                                            ttCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           >  2:32787/1‑90 (MQ=255)                                                                            ttCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           <  1:32787/90‑1 (MQ=255)                                                                                  tCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGaa                                                    <  1:654805/139‑1 (MQ=255)                                                                                    cTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACa                                                  <  1:322695/139‑1 (MQ=255)                                                                                    cTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACa                                                  >  2:568252/1‑139 (MQ=255)                                                                                            cGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATg                                          <  2:664415/139‑1 (MQ=255)                                                                                                       tGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGc                               >  1:587838/1‑139 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACTcccccgct                                                                                                                                 >  1:179536/1‑37 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACTcccccgct                                                                                                                                 <  2:179536/40‑4 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                <  2:255699/88‑1 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                >  1:255699/1‑88 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                                                                             <  2:133234/91‑1 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                                                                             >  1:133234/1‑91 (MQ=255)                                                                                                         cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             <  1:584743/139‑1 (MQ=255)                                                                                                           ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  <  1:428880/68‑1 (MQ=255)                                                                                                           ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  >  2:428880/1‑68 (MQ=255)                                                                                                           ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                                                                             <  2:225071/89‑1 (MQ=255)                                                                                                           ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                                                                             >  1:225071/1‑89 (MQ=255)                                                                                                           ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcg                           >  2:138008/1‑139 (MQ=255)                                                                                                           ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcg                           >  1:765941/1‑139 (MQ=255)                                                                                                            gCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcgc                          <  2:765941/139‑1 (MQ=255)                                                                                                                          gcgcGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    <  1:339026/67‑1 (MQ=255)                                                                                                                          gcgcGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    >  2:339026/1‑67 (MQ=255)                                                                                                                                    aaaCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCGGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  2:314140/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       ATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029588‑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: ‑