Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,229	A→G	S143S (AGT→AGC)	yijF ←	DUF1287 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,229	0	A	G	100.0%	222.4 / NA	63	S143S (AGT→AGC)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/30);  total (33/30)

CTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAATTC  >  NC_000913/4137093‑4137361                                                                                                                                         |                                                                                                                                     cttcttCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                                    >  2:1610167/1‑139 (MQ=255)         tGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                            >  2:1375561/1‑139 (MQ=255)                cTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                     >  1:1518424/1‑139 (MQ=255)                    gcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                 <  1:1623199/139‑1 (MQ=255)                     ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  1:1447289/139‑1 (MQ=255)                     ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  1:534721/139‑1 (MQ=255)                     ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  1:763221/139‑1 (MQ=255)                     ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  2:391967/139‑1 (MQ=255)                     ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  2:548274/139‑1 (MQ=255)                       aaTGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgt                                                                                                              <  1:147434/139‑1 (MQ=255)                       aaTGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgt                                                                                                              <  1:460943/139‑1 (MQ=255)                             aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                            >  2:476905/1‑119 (MQ=255)                             aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                            <  1:476905/119‑1 (MQ=255)                             aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                >  1:902793/1‑131 (MQ=255)                             aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  2:902793/131‑1 (MQ=255)                                    aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTggg                                                                                                                                 >  2:626326/1‑107 (MQ=255)                                    aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTggg                                                                                                                                 <  1:626326/107‑1 (MQ=255)                                    aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  2:1034128/124‑1 (MQ=255)                                    aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                >  1:1034128/1‑124 (MQ=255)                                     ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg                                                                                                >  2:13284/1‑139 (MQ=255)                                     ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg                                                                                                >  2:951175/1‑139 (MQ=255)                                        gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTc                                                                                                                             >  1:860174/1‑107 (MQ=255)                                        gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTc                                                                                                                             <  2:860174/107‑1 (MQ=255)                                            cgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc                                                                                         >  2:450486/1‑139 (MQ=255)                                             gTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc                                                                                        >  1:512654/1‑139 (MQ=255)                                                        ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATACAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGa                                                                             >  2:1174939/1‑139 (MQ=255)                                                               ccATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCa                                                                      <  1:1610167/139‑1 (MQ=255)                                                                 aTCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCAcc                                                                    <  1:1375561/139‑1 (MQ=255)                                                                       aaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGC‑‑GGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTg                                                            <  1:1235721/139‑1 (MQ=255)                                                                        accaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGt                                                             <  1:1618242/139‑1 (MQ=255)                                                                         ccaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTg                                                            <  1:1199215/139‑1 (MQ=255)                                                                         ccaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTg                                                            <  2:512654/139‑1 (MQ=255)                                                                            ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATc                                                         >  1:815929/1‑139 (MQ=255)                                                                            ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGACTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATc                                                         >  1:1564115/1‑139 (MQ=255)                                                                                     gtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCt                                                >  1:1027806/1‑139 (MQ=255)                                                                                       gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGt                                              >  1:180837/1‑139 (MQ=255)                                                                                       gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGt                                              >  2:181517/1‑139 (MQ=255)                                                                                       gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGt                                              >  2:816731/1‑139 (MQ=255)                                                                                       gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGt                                              >  1:1601313/1‑139 (MQ=255)                                                                                           gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                   <  1:1560716/50‑1 (MQ=255)                                                                                           gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                   >  2:1560716/1‑50 (MQ=255)                                                                                           gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCgg                                          >  1:15336/1‑139 (MQ=255)                                                                                           gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCgg                                          >  1:1638598/1‑139 (MQ=255)                                                                                            cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTggg                                                                                                                  <  2:1259504/66‑1 (MQ=255)                                                                                            cccGTTATCCAGTCGCCAGGAGACAATATCGCCAGCTTGATAGTCGCTGGGGTTCTTGCTGGTggg                                                                                                                  >  1:1259504/1‑66 (MQ=255)                                                                                            cccGTTATCCAGTCGCCAGGAGACAATATCGACCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTCAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCggg                                         >  1:755944/1‑139 (MQ=255)                                                                                                   tCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTa                                  >  1:1166839/1‑139 (MQ=255)                                                                                                      aGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACt                               >  2:1596620/1‑139 (MQ=255)                                                                                                       gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgtg                                                                                                             >  1:1011185/1‑60 (MQ=255)                                                                                                       gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgtg                                                                                                             <  2:1011185/60‑1 (MQ=255)                                                                                                            cAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtt                         <  2:1518424/139‑1 (MQ=255)                                                                                                                 gaCAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGc                    <  1:450486/139‑1 (MQ=255)                                                                                                                    aaTATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc                                                                                         <  2:1407298/67‑1 (MQ=255)                                                                                                                    aaTATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc                                                                                         >  1:1407298/1‑67 (MQ=255)                                                                                                                      tatCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTa               <  1:969559/139‑1 (MQ=255)                                                                                                                      tatCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTa               <  1:181517/139‑1 (MQ=255)                                                                                                                         cGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                >  2:1545151/1‑39 (MQ=255)                                                                                                                         cGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                <  1:1545151/39‑1 (MQ=255)                                                                                                                           cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAg          >  2:1142603/1‑139 (MQ=255)                                                                                                                             cGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCa        <  2:1601313/139‑1 (MQ=255)                                                                                                                              gCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCaa       >  2:1653184/1‑139 (MQ=255)                                                                                                                                ttGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTGTATCACGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTACGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCaaaa     >  1:1106721/1‑139 (MQ=255)                                                                                                                                   aTAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAAttc  <  2:1027806/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     CTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAATTC  >  NC_000913/4137093‑4137361

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