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

CATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029592‑2029858                                                                                                                                     |                                                                                                                                       cATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACATAAACTcccccgc                                                                                                                                  >  2:727360/1‑137 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:897173/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:967379/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:723241/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  1:971464/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:111746/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:1245334/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:27419/1‑139 (MQ=255)      ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGt                                                                                                                            >  2:1325167/1‑139 (MQ=255)           tCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCgct                                                                                                                       >  1:125019/1‑139 (MQ=255)                  aTTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGCTGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACCCAAACT‑CCCCCGCTACGTCCGCTGCTCAAt                                                                                                                >  1:1041375/1‑139 (MQ=255)                        tCATTCGGATTGGCGTGGAAAGCTGGTCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgcg                                                                                                          <  1:934591/139‑1 (MQ=255)                        tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgcg                                                                                                          <  1:111746/139‑1 (MQ=255)                        tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgcg                                                                                                          <  2:189085/139‑1 (MQ=255)                        tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCcgcg                                                                                                          >  1:1334312/1‑139 (MQ=255)                         cATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCgcgc                                                                                                         >  2:680114/1‑139 (MQ=255)                              ggATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGAtt                                                                                                    <  2:897173/139‑1 (MQ=255)                                aTTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  >  2:323684/1‑139 (MQ=255)                                 ttGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGc                                                                                                 <  2:971464/139‑1 (MQ=255)                                      gTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGtt                                                                                            <  1:1263602/139‑1 (MQ=255)                                            aGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgct                                                                                      <  2:108616/139‑1 (MQ=255)                                            aGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgct                                                                                      <  1:1325167/139‑1 (MQ=255)                                                           atgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTAc                                                                                                                              >  1:1207003/1‑84 (MQ=255)                                                           atgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTAc                                                                                                                              <  2:1207003/84‑1 (MQ=255)                                                            tgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGc                                                                      <  1:1245334/139‑1 (MQ=255)                                                                ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGAt                                                                                                     <  2:399043/104‑1 (MQ=255)                                                                ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGAt                                                                                                     >  1:399043/1‑104 (MQ=255)                                                                ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgc                                                                                       <  1:499765/118‑1 (MQ=255)                                                                ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgc                                                                                       >  2:499765/1‑118 (MQ=255)                                                                ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                                                                  >  1:523849/1‑139 (MQ=255)                                                                            aTTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCg                                                      <  1:323684/139‑1 (MQ=255)                                                                                   gCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATg                                                                                                               <  2:549799/75‑1 (MQ=255)                                                                                   gCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATg                                                                                                               >  1:549799/1‑75 (MQ=255)                                                                                       ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                                                                       >  2:497093/1‑48 (MQ=255)                                                                                       ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                                                                       <  1:497093/48‑1 (MQ=255)                                                                                       ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         >  1:343925/1‑61 (MQ=255)                                                                                       ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         <  2:343925/61‑1 (MQ=255)                                                                                                  ttGCTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCgg                                >  2:1106950/1‑139 (MQ=255)                                                                                                     cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGt                                                                                             >  2:873294/1‑75 (MQ=255)                                                                                                     cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGt                                                                                             <  1:873294/75‑1 (MQ=255)                                                                                                     cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTa                                                                           <  2:525220/93‑1 (MQ=255)                                                                                                     cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTa                                                                           >  1:525220/1‑93 (MQ=255)                                                                                                     cTGGCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  2:804375/1‑139 (MQ=255)                                                                                                       ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         <  2:583393/45‑1 (MQ=39)                                                                                                       ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCg                                                                                                                         >  1:583393/1‑45 (MQ=39)                                                                                                       ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcg                           >  2:1378695/1‑139 (MQ=255)                                                                                                       ggCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCT‑GCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcgc                          >  2:1317324/1‑139 (MQ=255)                                                                                                        gCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCTGAACATCAAAATGTTCCTTCCGGCATcgc                          >  2:995717/1‑139 (MQ=255)                                                                                                        gCATTTTTGCTACTGCGCGGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCACTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcgc                          >  1:883269/1‑139 (MQ=255)                                                                                                                         cgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCaa         <  2:883269/139‑1 (MQ=255)                                                                                                                         cgGACACAAACT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCaa         >  1:857084/1‑139 (MQ=255)                                                                                                                                aaaCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:680114/139‑1 (MQ=255)                                                                                                                                aaaCT‑CCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:727360/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029592‑2029858

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