Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,627,625	T→C	L37L (TTG→CTG)	ydfG →	NADP‑dependent 3‑hydroxy acid dehydrogenase, malonic semialdehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,627,625	0	T	C	100.0%	149.2 / NA	41	L37L (TTG→CTG)	ydfG	NADP‑dependent 3‑hydroxy acid dehydrogenase, malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/22);  total (19/22)

ACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAA  >  NC_000913/1627489‑1627754                                                                                                                                         |                                                                                                                                  acacacAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTg                                                                                                                                 <  2:111593/139‑1 (MQ=255)   acacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCa                                                                                                                               >  2:757905/1‑139 (MQ=255)   acacaAAAGCGTAGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCa                                                                                                                               >  2:420606/1‑139 (MQ=255)    cacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAg                                                                                                                              <  1:288519/139‑1 (MQ=255)              ttGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAg                                                                                                                    <  2:15524/139‑1 (MQ=255)                 aGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACg                                                                                                                 >  1:666452/1‑139 (MQ=255)                         tGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGga                                                                                                         >  2:731198/1‑139 (MQ=255)                          gagaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                        <  2:532434/139‑1 (MQ=255)                          gagaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                        <  2:118179/139‑1 (MQ=255)                            gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      <  2:512364/139‑1 (MQ=255)                            gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      <  1:731198/139‑1 (MQ=255)                                atCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                             <  2:93732/128‑1 (MQ=255)                                atCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                             >  1:93732/1‑128 (MQ=255)                                 tCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      <  2:879249/134‑1 (MQ=255)                                 tCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      >  1:879249/1‑134 (MQ=255)                                 tCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCt                                                                                                 >  2:425941/1‑139 (MQ=255)                                                                    gAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCATGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              <  2:666452/139‑1 (MQ=255)                                                                    gAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              <  2:701596/139‑1 (MQ=255)                                                                    gAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              <  1:344436/139‑1 (MQ=255)                                                                                   cgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACg                                                                                                                 <  1:316222/73‑1 (MQ=255)                                                                                   cgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACg                                                                                                                 >  2:316222/1‑73 (MQ=255)                                                                                   cgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              >  1:142012/1‑124 (MQ=255)                                                                                   cgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              <  2:142012/124‑1 (MQ=255)                                                                                               caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              >  1:400642/1‑112 (MQ=255)                                                                                               caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              <  2:400642/112‑1 (MQ=255)                                                                                                         aaGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                         >  2:178273/1‑139 (MQ=255)                                                                                                             tATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                             <  1:57053/51‑1 (MQ=255)                                                                                                             tATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                             >  2:57053/1‑51 (MQ=255)                                                                                                             tATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGc                     >  1:931492/1‑139 (MQ=255)                                                                                                              aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAAt                                                                                                   >  2:257725/1‑60 (MQ=255)                                                                                                              aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAAt                                                                                                   <  1:257725/60‑1 (MQ=255)                                                                                                                    aCTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              <  2:931979/91‑1 (MQ=255)                                                                                                                    aCTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                              >  1:931979/1‑91 (MQ=255)                                                                                                                     cTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                        <  1:434471/48‑1 (MQ=255)                                                                                                                     cTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                        >  2:434471/1‑48 (MQ=255)                                                                                                                      tGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      >  1:822946/1‑49 (MQ=255)                                                                                                                      tGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      <  2:822946/49‑1 (MQ=255)                                                                                                                             cGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      <  2:625736/42‑1 (MQ=255)                                                                                                                             cGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                      >  1:625736/1‑42 (MQ=255)                                                                                                                                cAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTaa  >  2:154014/1‑139 (MQ=255)                                                                                                                                  ggAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGc                     >  1:198061/1‑118 (MQ=255)                                                                                                                                  ggAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGc                     <  2:198061/118‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  ACACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAA  >  NC_000913/1627489‑1627754

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