Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,671,917	Δ1 bp	coding (693/813 nt)	hcaB →	2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,671,914	0	G	.	100.0%	92.7 / NA	23	coding (690/813 nt)	hcaB	2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/14);  total (9/14)

GGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTACGGGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAG  >  NC_000913/2671783‑2672042                                                                                                                                    |                                                                                                                                 ggTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTACgg                                                                                                                                 <  1:461170/133‑1 (MQ=255)  gTATGGCCAGGGACCTGCGCGGCCCACATGCGATCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTACggg                                                                                                                                <  2:231813/133‑1 (MQ=255)        ccAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTAt                                                                                                                         <  1:468007/133‑1 (MQ=255)            cGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGa                                                                                                                     >  2:311802/1‑133 (MQ=255)                        cccACAGGCGATCGGCCAAAGTGAGACCGCGACAATGCAGTGTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATcgc                                                                                                         <  2:173300/133‑1 (MQ=255)                                  cTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCaataa                                                                                               >  1:52373/1‑133 (MQ=255)                                    cGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATc                                                                                             <  2:590611/133‑1 (MQ=255)                                    cGGGCAAAGTGAAACCTCGATAATGCAGTCTCGGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATc                                                                                             <  2:625134/133‑1 (MQ=255)                                       gCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATcgcg                                                                                          >  2:89001/1‑133 (MQ=255)                                             tGAAACCTGGATAATGCAGTCTATGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTaa                                                                                    <  1:223872/133‑1 (MQ=255)                                             tGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTaa                                                                                    >  2:32221/1‑133 (MQ=255)                                              gAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAg                                                                                   <  1:89001/133‑1 (MQ=255)                                               aaaCCTCGATAATGCTGTCTCTGAAGCCTGAGAAAATTGCCACCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGc                                                                                  <  2:271662/133‑1 (MQ=255)                                                              aGTCTCTGACGCGGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAAc                                                                   <  2:469946/133‑1 (MQ=255)                                                                tctcTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGc                                                                 >  1:418296/1‑133 (MQ=255)                                                                 ctctGACGCCGTAGATAATTGCTGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCACAATAATCGCGCATTAAGCGGTGTGATGATCAACGCt                                                                <  1:277582/133‑1 (MQ=255)                                                                  tctGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCAGCACAATACTCGGGCATTAAGCGGTGGGATAATCAACGcgg                                                               >  2:277582/1‑131 (MQ=255)                                                                     gACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATg                                                            <  1:311802/133‑1 (MQ=255)                                                                     gACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATg                                                            <  1:132214/133‑1 (MQ=255)                                                                            gagaAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGttt                                                     >  2:490109/1‑133 (MQ=255)                                                                                         ccATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGc                                        >  1:110023/1‑133 (MQ=255)                                                                                                         aTTTTTCCCGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGCATTCGCCACGTAGCgg                        <  2:110023/133‑1 (MQ=255)                                                                                                               cccGCAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGCATTCGCCACGTAGCggctgggc                  >  2:452637/1‑133 (MQ=255)                                                                                                                   cAACCGGCGGATTTTAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGa              <  2:52373/133‑1 (MQ=255)                                                                                                                               tttAC‑GGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAg  <  1:490109/133‑1 (MQ=255)                                                                                                                                    |                                                                                                                                 GGTATGGCCAGCGACCTGCGCGGCCCACAGGCGCTCGGGCAAAGTGAAACCTCGATAATGCAGTCTCTGACGCCGGAGAAAATTGCCGCCATTTTACCGCTGCAATTTTTCCCGCAACCGGCGGATTTTACGGGGCCGTATGTGATGTTGACATCGCGGCGCAATAATCGCGCATTAAGCGGTGTGATGATCAACGCTGATGCGGGTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAG  >  NC_000913/2671783‑2672042

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: ‑