Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	803,662	C→A	L54I (CTT→ATT)	ybhJ →	putative hydratase YbhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,662	0	C	A	100.0%	96.6 / NA	27	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/12);  total (15/12)

CCTTCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTC  >  NC_000913/803418‑803889                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    ccTTCCTCCATCTCGAATTATGTAGAGTATATCCATTTATTATCGTTCTGCGCACTTCACGGGGCGCAGATAGCTGGAGCATTTGAGGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAATAAGCCAAAAAAGGCACTATTGCCTGGTATATTCTCTCTTCGCATAATTCGGCCGGAAATATGGATAAAATTaaaa                                                                                                                                                                                                                               <  2:285727/251‑1 (MQ=255)        tgaaCTCGAAATATTTATAGTATATCCATTTATTATCTTTCTGCGCACTGAACGGTTCGCAGAGATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCTTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAATAAGCCAAAAAAGGCACTATTGCCTGGTCTATGCTCTCTTCGCATAATACGTCCTTAAATATGGATAAAATTAAAATTAAg                                                                                                                                                                                                                          <  2:439597/247‑1 (MQ=255)                 attattTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                                                                                                                                                                                               <  2:37504/251‑1 (MQ=255)                                       attatCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCTTAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                                                                                                                                                                                          <  2:418904/250‑1 (MQ=255)                                                                                cATTTGATGATAAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACg                                                                                                                                                >  1:350926/1‑251 (MQ=255)                                                                                    tgatgaTCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTc                                                                                                                                            >  2:168421/1‑251 (MQ=255)                                                                                     gatgatCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTc                                                                                                                                            <  2:29487/250‑1 (MQ=255)                                                                                        gatCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCg                                                                                                                                         <  1:157523/250‑1 (MQ=255)                                                                                                     gAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCt                                                                                                                           <  2:6619/251‑1 (MQ=255)                                                                                                                                    aaaTAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGatat                                                                                                                                                                                                <  1:148995/151‑1 (MQ=255)                                                                                                                                    aaaTAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGatat                                                                                                                                                                                                >  2:148995/1‑151 (MQ=255)                                                                                                                                      ataatTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCGCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAGCTCACTCTTCGCCGTcgg                                                                                          <  2:412587/251‑1 (MQ=255)                                                                                                                                                          cACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTCTTAACGGtgatg                                                                      >  2:4183/1‑251 (MQ=255)                                                                                                                                                              ggCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGCTGAcc                                                                   >  1:83794/1‑250 (MQ=255)                                                                                                                                                                cGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATg                                                                <  1:168421/251‑1 (MQ=255)                                                                                                                                                                  aaaTTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTATATTCTCTCTTCGCATAATCCGTCCGGAAATATGGATAACATTAAACTTAAGTTTGATTCATTAGACTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGtt                                                              >  1:421093/1‑251 (MQ=255)                                                                                                                                                                               agaagCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATcgg                                                 >  1:78430/1‑251 (MQ=255)                                                                                                                                                                                             gCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTAt                                   >  1:265633/1‑251 (MQ=255)                                                                                                                                                                                               aCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCCGCGTTATgg                                 >  2:169731/1‑251 (MQ=255)                                                                                                                                                                                                     gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCgggat                           >  1:30329/1‑248 (MQ=255)                                                                                                                                                                                                       cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGAttt                         >  2:212326/1‑251 (MQ=255)                                                                                                                                                                                                              aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTAt                                                                                 <  2:248475/188‑1 (MQ=255)                                                                                                                                                                                                              aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTAt                                                                                 >  1:248475/1‑188 (MQ=255)                                                                                                                                                                                                                 ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTCACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTttttttttgtgtct               >  2:257349/1‑243 (MQ=255)                                                                                                                                                                                                                              ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCAt                                                                 >  2:366710/1‑188 (MQ=255)                                                                                                                                                                                                                              ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCAt                                                                 <  1:366710/188‑1 (MQ=255)                                                                                                                                                                                                                               taCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGGCTCCGCATAtttggggc  >  1:376638/1‑244 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                    CCTTCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTC  >  NC_000913/803418‑803889

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