Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,267,564	T→C	Y151H (TAT→CAT)	tyrB →	tyrosine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,267,564	0	T	C	100.0%	58.1 / NA	18	Y151H (TAT→CAT)	tyrB	tyrosine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/10);  total (8/10)

GGCGCGGATTTCCTGAAACGCTACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGA  >  NC_000913/4267444‑4267691                                                                                                                         |                                                                                                                                ggCGCGGATTTCCTGAAACGCTACTTCCCGGAATCGGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACATACTCCTGGCATGACGAAGCGa                                                                                                                     >  1:225019/1‑133 (MQ=255) ggCGCGGATTTCCTGAAACGCTACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGa                                                                                                                     >  1:375699/1‑133 (MQ=255)   cgcgGATTTCCTGAAACGCTACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACt                                                                                                                   >  1:345478/1‑133 (MQ=255)           tCCTGAAATGCTACTTCCTGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGt                                                                                                           <  1:215989/133‑1 (MQ=255)            ccTGAAACGCTACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTg                                                                                                          <  1:512293/133‑1 (MQ=255)                    gCTACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATAAAGAAGCGACTAACGGCGTGCGCtttca                                                                                                  >  2:251504/1‑131 (MQ=255)                      tACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCTTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATg                                                                                                >  2:423493/1‑133 (MQ=255)                      tACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATAACGAAGCGACTGACTGCGTGCACTTTAATg                                                                                                >  1:241980/1‑133 (MQ=255)                             cGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGCCATGACGATGCGCCGACGGGCGTGCACTTTCGTGACCTggt                                                                                         >  2:251606/1‑131 (MQ=255)                                 aTCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCg                                                                                     <  2:345478/133‑1 (MQ=255)                                  tCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCTTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGa                                                                                    <  1:423493/133‑1 (MQ=255)                                           tGGGTCAGCCGTCCTACCTGGGAGGACCCCGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGaaaa                                                                           <  2:241980/133‑1 (MQ=255)                                                                                  ttCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACATGCCCGCGGTATTGTGATGCTACATCCATGtt                                    >  2:220775/1‑133 (MQ=255)                                                                                       cGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGGAGTATTGTGTTGCTACATCCATGTTGCcac                               >  1:229695/1‑133 (MQ=255)                                                                                       cGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCcac                               <  1:507901/133‑1 (MQ=255)                                                                                                 ttCGAAGTCAGTACTGACTCCTAGCATGACTAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACgg                     <  1:251606/133‑1 (MQ=255)                                                                                                         gAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAt             <  1:439042/133‑1 (MQ=255)                                                                                                         gAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAt             <  1:507106/133‑1 (MQ=255)                                                                                                         gAGTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCACCCATGTTGCCATAATCCAACGTGTGCCGAt             >  1:177030/1‑133 (MQ=255)                                                                                                           gTACTTACCCCTGGCATGACGAAGTGCCCAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtct           <  2:225019/133‑1 (MQ=255)                                                                                                           gTACTTACCCCTGGCATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGAtct           <  1:412206/133‑1 (MQ=255)                                                                                                                    ccTGGCATGAGGAAGGGACTGACGGCCTTCGCATTAATGACCTGTTGACGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGa  <  2:229695/133‑1 (MQ=255)                                                                                                                         |                                                                                                                                GGCGCGGATTTCCTGAAACGCTACTTCCCGGAATCAGGCGTCTGGGTCAGCGATCCTACCTGGGAAAACCACGTAGCAATATTCGCCGGGGCTGGATTCGAAGTGAGTACTTACCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGA  >  NC_000913/4267444‑4267691

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