Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,669,362	T→C	S111P (TCT→CCT)	hcaE →	putative 3‑phenylpropionate/cinnamate dioxygenase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,669,362	0	T	C	100.0%	98.8 / NA	28	S111P (TCT→CCT)	hcaE	putative 3‑phenylpropionate/cinnamate dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/13);  total (15/13)

TGTCGTAGTGCGTCAGAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGTCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGTTATAAAGGG  >  NC_000913/2669235‑2669481                                                                                                                                |                                                                                                                        tGTCGTAGTGCGTCAGAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTAt                                                                                                                    <  2:325074/133‑1 (MQ=255)       aGTGCGTCAGAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATTCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCAtt                                                                                                              >  1:93052/1‑133 (MQ=255)         tGCGTCAGAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTaa                                                                                                            >  2:327853/1‑133 (MQ=255)             tCAGAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGc                                                                                                        >  1:448804/1‑133 (MQ=255)                gAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAg                                                                                                     >  1:439865/1‑133 (MQ=255)                     aCGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGAt                                                                                                <  1:498852/133‑1 (MQ=255)                           gcaTCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGt                                                                                          >  1:393196/1‑133 (MQ=255)                            caTCAAGGCGCTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTa                                                                                         >  1:338258/1‑133 (MQ=255)                                        tctcAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCt                                                                             <  2:601002/133‑1 (MQ=255)                                        tctcAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCt                                                                             <  1:327853/133‑1 (MQ=255)                                        tctcAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCt                                                                             <  2:580488/133‑1 (MQ=255)                                               cAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCt                                                                      <  2:108353/133‑1 (MQ=255)                                                 aTGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTAATCGATGCACCACTCGAACCTAGCGCCTAc                                                                    >  1:286570/1‑133 (MQ=255)                                                  tGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTAcc                                                                   >  2:343/1‑133 (MQ=255)                                                          ccGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAggg                                                           <  2:327319/133‑1 (MQ=255)                                                                     gtgtGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATc                                                <  2:439865/133‑1 (MQ=255)                                                                              aTGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTGCCCACAAGGGTTGTGTAAATCACACTgggg                                       >  1:252781/1‑133 (MQ=255)                                                                              aTGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTgggg                                       <  2:338258/133‑1 (MQ=255)                                                                                    aTTGAGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTaaa                                 >  1:445527/1‑133 (MQ=255)                                                                                         ggCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCTACTGGGGACTAAACGAAg                            >  1:265059/1‑133 (MQ=255)                                                                                                    cgcCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTtgtgtg                 <  2:448804/133‑1 (MQ=255)                                                                                                    cgcCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTtgtgtg                 <  2:467564/133‑1 (MQ=255)                                                                                                    cgcCTTTACCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTtgtgtg                 <  2:9643/133‑1 (MQ=255)                                                                                                           aCCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGtt          >  2:382468/1‑133 (MQ=255)                                                                                                           aCCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGtt          >  2:447815/1‑133 (MQ=255)                                                                                                           aCCTGCCCGTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGtt          >  1:142279/1‑133 (MQ=255)                                                                                                              tGCCCGTATCACGGCTGGCCTTATGGGATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGTtata       <  2:286570/133‑1 (MQ=255)                                                                                                                   gTATCACGGCTGGCCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGTTATAAAggg  >  2:344128/1‑133 (MQ=255)                                                                                                                                |                                                                                                                        TGTCGTAGTGCGTCAGAAAGACGGCAGCATCAAGGCGTTTCTCAACCAATGCCGCCACCGGGCCATGCGTGTGAGTTATGCAGATTGCGGCAACACTCGCGCCTTTACCTGCCCGTATCACGGCTGGTCTTATGGCATTAACGGCGAGTTGATCGATGTACCGCTGGAACCTCGCGCCTACCCACAAGGGTTGTGTAAATCCCACTGGGGACTAAACGAAGTTCCTTGTGTGGAGAGTTATAAAGGG  >  NC_000913/2669235‑2669481

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