Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,565,464	C→A	L129F (TTG→TTT)	yjiM ←	putative 2‑hydroxyglutaryl‑CoA dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,565,464	0	C	A	100.0%	166.2 / NA	52	L129F (TTG→TTT)	yjiM	putative 2‑hydroxyglutaryl‑CoA dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/28);  total (24/28)

TTTCAGCGCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAAAGTAGAAGTAGGGGCA  >  NC_000913/4565332‑4565574                                                                                                                                     |                                                                                                               tttCAGCGCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          >  2:709322/1‑139 (MQ=255)           aTGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCTTCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGaa                                                                                                <  2:577430/139‑1 (MQ=255)                    cgcAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTaa                                                                                       >  1:72101/1‑139 (MQ=255)                    cgcAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTaa                                                                                       >  1:1040678/1‑139 (MQ=255)                      cAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAAc                                                                                     <  2:72101/139‑1 (MQ=255)                       agagCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACt                                                                                    <  2:997380/139‑1 (MQ=255)                        gagCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGc                                                                                           <  1:1128593/131‑1 (MQ=255)                        gagCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGc                                                                                           >  2:1128593/1‑131 (MQ=255)                              cTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCAt                                                                             <  1:19566/139‑1 (MQ=255)                               ttCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATg                                                                            <  1:1270789/139‑1 (MQ=255)                                   cTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGaa                                                                                                <  2:816034/115‑1 (MQ=255)                                   cTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGaa                                                                                                >  1:816034/1‑115 (MQ=255)                                            tGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTaaa                                                                                      >  2:678101/1‑116 (MQ=255)                                            tGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTaaa                                                                                      <  1:678101/116‑1 (MQ=255)                                               ccAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCt                                                            <  1:709322/139‑1 (MQ=255)                                                cAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATcac                                                                                                     <  1:1329208/97‑1 (MQ=255)                                                cAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATcac                                                                                                     >  2:1329208/1‑97 (MQ=255)                                                 aaaaCGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGaa                                                                                                <  2:1119267/101‑1 (MQ=255)                                                 aaaaCGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGaa                                                                                                >  1:1119267/1‑101 (MQ=255)                                                        tcttctACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGt                                                   <  2:1040678/139‑1 (MQ=255)                                                          ttctACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCg                                                 <  1:824386/139‑1 (MQ=255)                                                               ccGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          <  2:9771/77‑1 (MQ=255)                                                               ccGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          >  1:9771/1‑77 (MQ=255)                                                               ccGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCa                                                                              <  1:1155927/105‑1 (MQ=255)                                                               ccGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCa                                                                              >  2:1155927/1‑105 (MQ=255)                                                                 gTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          <  2:146454/75‑1 (MQ=255)                                                                 gTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          >  1:146454/1‑75 (MQ=255)                                                                  ttttttGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCa                                                                              >  2:109846/1‑102 (MQ=255)                                                                  ttttttGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCa                                                                              <  1:109846/102‑1 (MQ=255)                                                                                                   acgcacgcGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCg                                                 >  1:838522/1‑98 (MQ=255)                                                                                                   acgcacgcGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCg                                                 <  2:838522/98‑1 (MQ=255)                                                                                                    cgcacgcGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGAt                        >  1:317912/1‑122 (MQ=255)                                                                                                    cgcacgcGAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGAt                        <  2:317912/122‑1 (MQ=255)                                                                                                         gcgAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          <  1:428791/35‑1 (MQ=255)                                                                                                         gcgAGGCATCGTCCTTAACGCTGTTGGGAAATTGc                                                                                                          >  2:428791/1‑35 (MQ=255)                                                                                                         gcgAGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAAAGTAGAAGTAGGGGCa  >  2:777002/1‑139 (MQ=255)                                                                                                            aGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAaagtagaagtag       >  2:897255/1‑131 (MQ=255)                                                                                                            aGGCATCGTCCTTAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAaagtagaagtag       <  1:897255/131‑1 (MQ=255)                                                                                                                       ttAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACt                                                                                    >  1:151421/1‑43 (MQ=255)                                                                                                                       ttAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACt                                                                                    <  2:151421/43‑1 (MQ=255)                                                                                                                        tAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACt                                                                                    >  1:125587/1‑42 (MQ=255)                                                                                                                        tAACGCTGTTGGGAAATTGCATCACATGAACAGGCTTAAACt                                                                                    <  2:125587/42‑1 (MQ=255)                                                                                                                             cTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAaagtagaagt         <  1:1214793/112‑1 (MQ=255)                                                                                                                             cTGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAaagtagaagt         >  2:1214793/1‑112 (MQ=255)                                                                                                                              tGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGttt                                       >  1:1331250/1‑81 (MQ=255)                                                                                                                              tGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGttt                                       <  2:1331250/81‑1 (MQ=255)                                                                                                                              tGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCg                                  >  1:1168265/1‑86 (MQ=255)                                                                                                                              tGTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCg                                  <  2:1168265/86‑1 (MQ=255)                                                                                                                               gTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAAAGTAGAAGTAGGGGc   <  2:979391/116‑1 (MQ=255)                                                                                                                               gTTGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGATTTCACCGACCACCAGATCCGAAAAGTAGAAGTAGGGGc   >  1:979391/1‑116 (MQ=255)                                                                                                                                ttGGGAAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCAt                                                                     <  2:813843/49‑1 (MQ=255)                                                                                                                                ttGGGAAATTGCATCACATGAACAGGCATAAAATCCGCCATGTATTCAt                                                                     >  1:813843/1‑49 (MQ=38)                                                                                                                                     |                                                                                                               TTTCAGCGCAATGGCATCGCGCAGAGCATCTTCGCTAATCTCGTGCCCAAAACGTTCTTCTACCGTTTTTTGCAAGCGCAGCATCTCGGCTTTCCATAACGCACGCGAGGCATCGTCCTTAACGCTGTTGGGCAATTGCATCACATGAACAGGCTTAAACTCCGCCATGTATTCATACATTTTCTTTTTGCCGTCGCAGGTGGTTTCACCGACCACCAGATCCGAAAAGTAGAAGTAGGGGCA  >  NC_000913/4565332‑4565574

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