Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,862,139	+C	coding (710/909 nt)	ygbJ →	putative L‑threonate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,862,135	1	.	C	100.0%	107.5 / NA	32	T236P (ACC→CCC)	ygbJ	putative L‑threonate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/18);  total (14/18)

TGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACA‑CCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGAATATGTTCACCAGCGCCAGTAACGCGGGT  >  NC_000913/2861996‑2862266                                                                                                                                             |                                                                                                                                    tgctgcCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTaca                                                                                                                                      <  1:726547/140‑1 (MQ=255)    tgcCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACA‑ccc                                                                                                                                  >  2:104081/1‑140 (MQ=255)      ccGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACA‑ccccc                                                                                                                                >  2:184942/1‑139 (MQ=255)       cGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACA‑cccccg                                                                                                                               >  1:797738/1‑138 (MQ=255)                         aCTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATAtt                                                                                                              <  2:797738/140‑1 (MQ=255)                                 cccGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTaa                                                                                                      <  2:433304/140‑1 (MQ=255)                                   cGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAgg                                                                                                    >  1:740170/1‑140 (MQ=255)                                            ttCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTc                                                                                           <  1:33487/139‑1 (MQ=255)                                                           atgtATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAg                                                                            >  1:134083/1‑140 (MQ=255)                                                                tGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCa                                                                         <  1:125955/138‑1 (MQ=255)                                                                tGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCa                                                                         >  2:125955/1‑138 (MQ=255)                                                                   cgtcgtGACCAATGCCGGCGGATATTCCTGGCTGTGCGACAACCGGATCCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGcc                                                                    <  2:215271/140‑1 (MQ=255)                                                                   cgtcgtGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGcc                                                                    <  2:740170/140‑1 (MQ=255)                                                                   cgtcgtGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGATAAGGATATTGATCGGGTTGCAGATACAGCCAAAGcc                                                                    >  1:215271/1‑140 (MQ=255)                                                                     tcgtGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCt                                                                  >  2:584980/1‑140 (MQ=255)                                                                      cgtGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTg                                                                 >  2:670841/1‑140 (MQ=255)                                                                      cgtGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTg                                                                 >  2:31549/1‑140 (MQ=255)                                                                        tGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCa                                                                         <  1:200655/130‑1 (MQ=255)                                                                        tGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCa                                                                         >  2:200655/1‑130 (MQ=255)                                                                        tGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATAGTTTTGTTAAGGACCTTTGTCTGGTTGCCGATACAGCCAAAGCCCTGCa                                                               >  2:267745/1‑140 (MQ=255)                                                                              aTGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTccc                                                         >  1:701712/1‑140 (MQ=255)                                                                                gccgccGGAAATTCCTGGATGTCCGAAAATCGGATGCATCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGc                                                       <  1:267745/140‑1 (MQ=255)                                                                                     cgGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCa                                                  <  2:599989/140‑1 (MQ=255)                                                                                     cgGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCa                                                  <  2:58211/140‑1 (MQ=255)                                                                                          aTTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGc                                             <  1:584980/140‑1 (MQ=255)                                                                                                    gTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCa                                   <  1:104081/140‑1 (MQ=255)                                                                                                    gTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCa                                   <  2:134083/140‑1 (MQ=255)                                                                                                       cGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAg                                     <  2:441512/135‑1 (MQ=255)                                                                                                       cGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAg                                     >  1:441512/1‑135 (MQ=255)                                                                                                       cGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTg                                >  2:66263/1‑140 (MQ=255)                                                                                                       cGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTg                                >  2:764547/1‑140 (MQ=255)                                                                                                        gAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGa                               <  1:184942/140‑1 (MQ=255)                                                                                                        gAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGa                               <  2:701712/140‑1 (MQ=255)                                                                                                          aaaCCGGATGCGTCATGTGGTGGATGGCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGAat                             <  1:66263/140‑1 (MQ=255)                                                                                                                                   ggCGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGAATATGTTCACCAGCGCCAGTAACGCgg    >  2:592599/1‑140 (MQ=255)                                                                                                                                     cGATTACACCCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGAATATGTTCACCAGCGCCAGTAACGCGGGt  <  1:764547/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                    TGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACA‑CCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTGCACTTCCCGCTGCCATTGGCCTCAACAGCATTGAATATGTTCACCAGCGCCAGTAACGCGGGT  >  NC_000913/2861996‑2862266

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