Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	375,307	C→T	intergenic (+426/‑152)	mhpE → / → mhpT	4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	375,307	0	C	T	100.0%	152.4 / NA	50	intergenic (+426/‑152)	mhpE/mhpT	4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/25);  total (25/25)

GTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTC  >  NC_000913/375183‑375441                                                                                                                             |                                                                                                                                       gTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCAt                                                                                                                          <  1:1025803/139‑1 (MQ=25)     gCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGc                                                                                                                      >  2:1025803/1‑139 (MQ=21)         gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCacac                                                                                                                  <  1:933705/139‑1 (MQ=25)         gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCacac                                                                                                                  <  1:295994/139‑1 (MQ=37)          aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACg                                                                                                                 <  1:1733629/139‑1 (MQ=25)               aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGc                                                                                                            >  2:831602/1‑139 (MQ=21)               aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGc                                                                                                            >  1:122586/1‑139 (MQ=25)                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTATGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                                                                       >  2:308632/1‑139 (MQ=21)                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCAt                                                                                                                                 >  1:961849/1‑113 (MQ=17)                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCAt                                                                                                                                 <  2:961849/113‑1 (MQ=17)                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                                                                       >  2:1343998/1‑139 (MQ=25)                      aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGc                                                                                                     <  2:1290738/139‑1 (MQ=21)                             cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                                                                 <  1:1767855/136‑1 (MQ=37)                             cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                                                                 >  2:1767855/1‑136 (MQ=37)                                  ggCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCa                                                                                                                               >  1:1669645/1‑101 (MQ=18)                                  ggCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCa                                                                                                                               <  2:1669645/101‑1 (MQ=18)                                       cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                       <  2:1282632/136‑1 (MQ=255)                                       cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                       >  1:1282632/1‑136 (MQ=255)                                        tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                                                                   >  1:280537/1‑139 (MQ=255)                                          cgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                      <  2:598644/134‑1 (MQ=255)                                          cgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                      >  1:598644/1‑134 (MQ=255)                                          cgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                                                                 <  2:280537/139‑1 (MQ=255)                                          cgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                                                                 <  1:831602/139‑1 (MQ=255)                                          cgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                                                                 <  1:1343998/139‑1 (MQ=255)                                             aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTAc                                                                              >  2:357012/1‑139 (MQ=255)                                                 cTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                                                          >  1:722393/1‑139 (MQ=255)                                                    ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                                                          >  2:1703475/1‑136 (MQ=255)                                                    ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                                                          <  1:1703475/136‑1 (MQ=255)                                                           cggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                                                                >  2:160649/1‑139 (MQ=255)                                                            ggcgTAAACGCCTTATCCGGCCTATGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:308632/139‑1 (MQ=255)                                                            ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  2:122586/139‑1 (MQ=255)                                                              cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                       >  2:1189884/1‑113 (MQ=255)                                                              cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                                                                       <  1:1189884/113‑1 (MQ=255)                                                                                    cGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg                                       <  1:160649/139‑1 (MQ=255)                                                                                    cGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcg                                       <  1:306101/139‑1 (MQ=255)                                                                                     ggATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTAcgc                                      >  1:1336638/1‑139 (MQ=255)                                                                                       aTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCgg                                    <  2:722393/139‑1 (MQ=255)                                                                                        tGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGc                                   <  2:1114237/139‑1 (MQ=255)                                                                                          gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTg                                 >  2:562802/1‑139 (MQ=255)                                                                                             cgGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  2:1104422/106‑1 (MQ=255)                                                                                             cgGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  1:1104422/1‑106 (MQ=255)                                                                                                              gATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGAtt             >  1:90593/1‑139 (MQ=255)                                                                                                                    aCGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGt                          >  1:1589867/1‑120 (MQ=255)                                                                                                                    aCGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGt                          <  2:1589867/120‑1 (MQ=255)                                                                                                                    aCGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTg       >  1:966217/1‑139 (MQ=255)                                                                                                                      gcgTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  2:209680/81‑1 (MQ=255)                                                                                                                      gcgTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  1:209680/1‑81 (MQ=255)                                                                                                                         ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGcc                                                                      >  1:614373/1‑71 (MQ=255)                                                                                                                         ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGcc                                                                      <  2:614373/71‑1 (MQ=255)                                                                                                                         ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTc  >  1:1536094/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                       GTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTC  >  NC_000913/375183‑375441

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