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.5 / NA	49	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 (28/21);  total (28/21)

GCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTC  >  NC_000913/375171‑375441                                                                                                                                          |                                                                                                                                       gcgcGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAgtgt                                                                                                                                      <  1:1016111/139‑1 (MQ=37)        aaTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCATAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCa                                                                                                                               <  2:177435/139‑1 (MQ=255)                  ccTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGca                                                                                                                     >  2:801922/1‑139 (MQ=21)                     gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCacac                                                                                                                  >  2:1309861/1‑139 (MQ=25)                     gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCacac                                                                                                                  >  1:167702/1‑139 (MQ=25)                            cgcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCAt                                                                                                                          <  1:800296/124‑1 (MQ=18)                            cgcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCAt                                                                                                                          >  2:800296/1‑124 (MQ=17)                                ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                                                                       >  2:1151825/1‑139 (MQ=25)                                         cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAAc                                                                                              >  2:127485/1‑139 (MQ=255)                                                    tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGTATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                                                                   >  1:1100328/1‑139 (MQ=255)                                                      cgcACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                                                                 >  2:124087/1‑139 (MQ=255)                                                             cTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAttt                                                                          >  1:1530017/1‑139 (MQ=255)                                                                ccGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGcaca                                                                                                                   >  2:1192806/1‑95 (MQ=255)                                                                ccGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGcaca                                                                                                                   <  1:1192806/95‑1 (MQ=255)                                                                ccGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATATCGc                                                                       >  2:1208834/1‑139 (MQ=255)                                                                     tGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                                                                 >  2:917991/1‑108 (MQ=255)                                                                     tGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                                                                 <  1:917991/108‑1 (MQ=255)                                                                     tGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTa                                                                  <  1:801922/139‑1 (MQ=255)                                                                       cggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGc                                                                                                                                   >  1:434887/1‑72 (MQ=25)                                                                       cggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGc                                                                                                                                   <  2:434887/72‑1 (MQ=255)                                                                       cggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAgg                                                                >  1:1711352/1‑139 (MQ=255)                                                                        ggggttaaaggcTTTTCCGGGCTACGGGATGGGGCGGGATTTGGTAGG‑CTGGTTAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:358601/128‑1 (MQ=18)                                                                        ggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:1151825/139‑1 (MQ=255)                                                                        ggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:857824/139‑1 (MQ=255)                                                                        ggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:1309861/139‑1 (MQ=255)                                                                        ggcgTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  2:1530017/139‑1 (MQ=255)                                                                                    ttATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  1:1146192/1‑127 (MQ=255)                                                                                    ttATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  2:1146192/127‑1 (MQ=255)                                                                                      aTCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGt                                                 >  2:100594/1‑139 (MQ=255)                                                                                                c‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:629175/115‑1 (MQ=255)                                                                                                c‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  2:629175/1‑115 (MQ=255)                                                                                                      ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCt                                  >  1:706517/1‑139 (MQ=255)                                                                                                       gcgcGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTg                                 >  1:670924/1‑139 (MQ=255)                                                                                                          cgGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGTATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                                                     <  2:1767083/116‑1 (MQ=255)                                                                                                          cgGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:1432924/106‑1 (MQ=255)                                                                                                          cgGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  2:1432924/1‑106 (MQ=255)                                                                                                          cgGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTaaa                                                     >  1:1767083/1‑116 (MQ=255)                                                                                                                         cTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                                                                       <  2:199049/51‑1 (MQ=255)                                                                                                                         cTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                                                                                       >  1:199049/1‑51 (MQ=255)                                                                                                                            aTAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATg                                                                                                      >  2:1490618/1‑49 (MQ=255)                                                                                                                            aTAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATg                                                                                                      <  1:1490618/49‑1 (MQ=255)                                                                                                                                 aCGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCAtt                                                        <  2:711028/90‑1 (MQ=255)                                                                                                                                 aCGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCAtt                                                        >  1:711028/1‑90 (MQ=255)                                                                                                                                  cgcgTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGc      >  2:1020870/1‑139 (MQ=255)                                                                                                                                     gTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  2:1235712/1‑79 (MQ=255)                                                                                                                                     gTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:1235712/79‑1 (MQ=255)                                                                                                                                     gTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACt   <  2:167702/139‑1 (MQ=255)                                                                                                                                     gTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACt   <  2:706517/139‑1 (MQ=255)                                                                                                                                      ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTc  >  2:353208/1‑139 (MQ=255)                                                                                                                                      ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTc  >  2:373361/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                       GCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC‑GGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTC  >  NC_000913/375171‑375441

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