Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,615,859	T→G	100%	T403P (ACC→CCC)	yjjI ←	DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,615,859	0	T	G	100.0%	275.5 / NA	80	T403P (ACC→CCC)	yjjI	DUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (41/39);  total (41/39)

CTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTG  >  NC_000913/4615764‑4615958                                                                                                |                                                                                                    cTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATc                                                                                                 <  1:1363212/100‑1 (MQ=255)  tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGCTGCCGGAATAATACTAAGCATGATGCAGGGCGACAGCTTGCAGATGGGGGAt                                                                                                  >  2:104301/1‑98 (MQ=255)  tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                                <  1:1982762/100‑1 (MQ=255)  tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                                <  2:16346/100‑1 (MQ=255)   gCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCgg                                                                                               >  1:1916995/1‑100 (MQ=255)   gCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATATCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCgg                                                                                               >  2:1363562/1‑100 (MQ=255)    cGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGa                                                                                              <  2:2199927/100‑1 (MQ=255)     gggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGAt                                                                                             >  1:897776/1‑100 (MQ=255)     gggTTACGCTTGCTGGTTACAACGAGCGTCTGAATGTCGCTGATGCCGGCATAATAATAAGCATACTGCAGGGCGAAAGTTTGCCGATGGGGTATCgaa                                                                                              >  2:884182/1‑97 (MQ=255)      ggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGggatcggatc                                                                                            >  1:532990/1‑100 (MQ=255)       gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                           >  2:1356054/1‑100 (MQ=255)        ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                          >  1:765813/1‑100 (MQ=255)          aCGTTTGATGGTTTCGTCGAGCGTCAGTATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                        <  1:1486569/100‑1 (MQ=255)          aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                        <  1:1636833/100‑1 (MQ=255)          aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                        <  1:371653/100‑1 (MQ=255)           cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                       >  1:158636/1‑100 (MQ=255)           cGTTTGATGGTTTCGTCGAGCGTCAGAAGGTCGCTGATGCCGGAATCATAATAAGCATGATGCGGGGCGACAGTCTGCAGATGGCGGATCGCATCTGGcc                                                                                       <  2:2139346/100‑2 (MQ=255)            gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATTATGCGGGGCGACAGTTTGCAGATGGGGGCTCGGATCTGGCTc                                                                                      <  1:1731825/100‑1 (MQ=255)                    gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACTGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                              <  2:1366551/100‑1 (MQ=255)                    gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                              >  2:2465978/1‑100 (MQ=255)                    gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                              >  1:2470019/1‑100 (MQ=255)                    gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                              >  1:1933821/1‑100 (MQ=255)                    gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                              >  2:1455118/1‑100 (MQ=255)                        cgtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACgg                                                                          >  2:2471467/1‑100 (MQ=255)                         gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                         <  1:1099788/100‑1 (MQ=255)                          tcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCa                                                                        <  1:1346600/100‑1 (MQ=255)                           cgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaa                                                                       <  2:848225/100‑1 (MQ=255)                             aGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAc                                                                     >  1:2496656/1‑100 (MQ=255)                                gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                                  >  2:1520442/1‑100 (MQ=255)                                 tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                                  >  2:934279/1‑99 (MQ=255)                                   aGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                                               >  2:770707/1‑100 (MQ=255)                                      aTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCg                                                            >  1:1057498/1‑100 (MQ=255)                                      aTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCg                                                            >  2:877885/1‑100 (MQ=255)                                           gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                       <  1:1639682/100‑1 (MQ=255)                                           gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                       >  1:1314576/1‑100 (MQ=255)                                           gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                       >  1:1251715/1‑100 (MQ=255)                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                                  <  2:1856856/100‑1 (MQ=255)                                                tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                                  >  1:995837/1‑100 (MQ=255)                                                tGCCAGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                                   >  1:1867425/1‑99 (MQ=255)                                                 gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCTGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                                 <  2:758961/100‑1 (MQ=255)                                                 gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                                 <  2:2461883/100‑1 (MQ=255)                                                  ccGGAATAATAATACGCATGATGCGGGGCGACAGTTGGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                                <  2:1039119/100‑1 (MQ=255)                                                  ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                                >  2:2072720/1‑100 (MQ=255)                                                    ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATCCGGCAAACGCGCGCCCGGCGTGGTGCCGata                                              >  2:2214338/1‑100 (MQ=255)                                                    ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGata                                              <  1:1479714/100‑1 (MQ=255)                                                        taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                          >  2:609542/1‑100 (MQ=255)                                                          ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                        >  1:1432752/1‑100 (MQ=255)                                                           taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                       <  1:1069771/100‑1 (MQ=255)                                                           taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                       <  1:422866/100‑1 (MQ=255)                                                            aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGCAc                                       >  1:1172931/1‑99 (MQ=255)                                                             ataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTg                                     >  1:2050243/1‑100 (MQ=255)                                                              taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                                    <  1:357971/100‑1 (MQ=255)                                                              taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                                    <  2:1646612/100‑1 (MQ=255)                                                              taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                                    <  2:1647635/100‑1 (MQ=255)                                                              taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCACCCGGCGTGGTGCCGATATCGCAACTGa                                    >  2:487923/1‑100 (MQ=255)                                                              taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGAGCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                                    >  2:1395750/1‑100 (MQ=255)                                                                aGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATc                                  >  1:1755973/1‑100 (MQ=255)                                                                 gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                                 <  1:1433471/100‑1 (MQ=255)                                                                  cATGTTGCGGCGCGACAGTTTGCAGATAGGGGCTCGGATATGGCCCCTCGGCATACCTGAAACGCGCGCCCGGCGTGGGGCCGatcttcggac                                       >  2:1780884/1‑85 (MQ=255)                                                                  cATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATccc                                <  2:1592018/100‑1 (MQ=255)                                                                   atgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc                               >  1:624800/1‑100 (MQ=255)                                                                    tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                              <  1:1615705/100‑1 (MQ=255)                                                                      aggCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGAc                            <  1:2146614/98‑1 (MQ=255)                                                                       tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACt                           >  2:1104077/1‑100 (MQ=255)                                                                        gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                          <  1:1881658/100‑1 (MQ=255)                                                                         cGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgt                         >  1:512306/1‑100 (MQ=255)                                                                         cGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgt                         <  2:1129791/100‑1 (MQ=255)                                                                          tctgCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtg                        <  2:122726/97‑1 (MQ=255)                                                                          ggggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtg                        <  1:2471481/100‑1 (MQ=255)                                                                          ggggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtg                        <  1:1880561/100‑1 (MQ=255)                                                                            ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                      >  1:1454736/1‑100 (MQ=255)                                                                             gCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtgaga                     >  1:122728/1‑97 (MQ=255)                                                                              cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                    <  1:298836/100‑1 (MQ=255)                                                                              cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                    <  1:2203377/100‑1 (MQ=255)                                                                               gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt                   <  1:1678109/100‑1 (MQ=255)                                                                                aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa                  >  1:776451/1‑100 (MQ=255)                                                                                aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa                  <  1:173622/100‑1 (MQ=255)                                                                                    tttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt              >  2:1320945/1‑100 (MQ=255)                                                                                      tGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg            <  2:336842/100‑1 (MQ=255)                                                                                          gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg        >  2:1769036/1‑100 (MQ=255)                                                                                           aTGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGt       <  2:438757/100‑1 (MQ=255)                                                                                             ggggGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGttt     >  2:1056707/1‑100 (MQ=255)                                                                                                tGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCCGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTg  <  1:412475/100‑1 (MQ=255)                                                                                                |                                                                                                    CTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTG  >  NC_000913/4615764‑4615958

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