Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	858,789	T→C	K91K (AAA→AAG)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,789	0	T	C	100.0%	129.5 / NA	35	K91K (AAA→AAG)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/20);  total (15/20)

AGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTG  >  NC_000913/858641‑858930                                                                                                                                                     |                                                                                                                                              aGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCctt                                                                                                                                             >  2:278156/1‑151 (MQ=255)     ggggTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATc                                                                                                                                         <  2:200607/151‑1 (MQ=255)      gggTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCa                                                                                                                                        <  2:83744/151‑1 (MQ=255)          aaaTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGc                                                                                                                                    <  1:396012/151‑1 (MQ=255)           aaTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGca                                                                                                                                   <  2:21194/151‑1 (MQ=255)                cGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACaa                                                                                                                                <  2:61901/149‑1 (MQ=255)                  ggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACcgc                                                                                                                            <  1:438123/151‑1 (MQ=255)                  ggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACcgc                                                                                                                            <  1:629134/151‑1 (MQ=255)                       tGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGt                                                                                                                       >  1:569653/1‑151 (MQ=255)                        gAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGt                                                                                                                       <  2:11159/150‑1 (MQ=255)                             gcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATc                                                                                                                 >  1:537775/1‑151 (MQ=255)                                gAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGt                                                                                                              >  1:390678/1‑151 (MQ=255)                                   cGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTcg                                                                                                           >  1:29970/1‑151 (MQ=255)                                       ccGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCgct                                                                                                       <  1:477818/151‑1 (MQ=255)                                         gaagaaATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCgctgc                                                                                                     >  2:487100/1‑151 (MQ=255)                                              aaTTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGTTTGTTCGCGGGGGAAATGGCCCTTTTTCACCGCACAACCGCGACGTTTAACCTGTTCGCGCTGctct                                                                                                  <  1:552867/149‑1 (MQ=255)                                               aTTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCa                                                                                               <  1:151735/151‑1 (MQ=255)                                                    cGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAAc                                                                                          <  1:583300/151‑1 (MQ=255)                                                    cGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAAc                                                                                          >  2:329636/1‑151 (MQ=255)                                                                       cTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCg                                                                       >  1:456818/1‑151 (MQ=255)                                                                       cTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCg                                                                       >  1:408609/1‑151 (MQ=255)                                                                            cAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATa                                                                  <  2:371051/151‑1 (MQ=255)                                                                                     tGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCa                                                          <  2:128359/150‑1 (MQ=255)                                                                                      gCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCAcc                                                        >  1:540974/1‑151 (MQ=255)                                                                                       cGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCg                                                       >  1:509150/1‑151 (MQ=255)                                                                                            cAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGTGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGcc                                                    >  2:298450/1‑149 (MQ=255)                                                                                            cAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAg                                                  <  2:295110/151‑1 (MQ=255)                                                                                                    ccAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCAt                                           <  1:134698/150‑1 (MQ=255)                                                                                                         aTCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTTTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCt                                     <  2:536362/151‑1 (MQ=255)                                                                                                         aTCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCt                                     >  2:264316/1‑151 (MQ=255)                                                                                                           cGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCtgt                                   <  2:269545/151‑1 (MQ=255)                                                                                                                 aTCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTg                             <  1:259053/151‑1 (MQ=255)                                                                                                                     cAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTg                         >  1:197361/1‑151 (MQ=255)                                                                                                                        ccTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGt                      >  1:153620/1‑151 (MQ=255)                                                                                                                        ccTAGCGCTTGTTCGCGGGGGAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGt                      <  2:167597/151‑1 (MQ=255)                                                                                                                                            gAAATGGCCCTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTg  >  1:305451/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                              AGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAATCTGTTCGCGCTGCTCTGCAGTAACATGACCTGCTTTGATATCGGCATAAGACAGCACCGGCCAGACGGCATCGCCCTGTGCTTTGAGTGCGTTGATTTCAGCCACTCGCGTG  >  NC_000913/858641‑858930

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