Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,051,733	A→G	Y17C (TAC→TGC)	ymcF →	protein YmcF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,051,733	0	A	G	100.0%	176.7 / NA	50	Y17C (TAC→TGC)	ymcF	protein YmcF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/21);  total (29/21)

ACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAA  >  NC_000913/1051597‑1051869                                                                                                                                         |                                                                                                                                         aCGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCg                                                                                                                                       <  2:243432/140‑1 (MQ=255)               gTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCg                                                                                                                         <  2:45461/140‑1 (MQ=255)                      tttCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGAc                                                                                                                  <  1:684152/140‑1 (MQ=255)                        tCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACgg                                                                                                                >  2:221331/1‑140 (MQ=255)                                 cAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAAcc                                                                                                       <  1:221331/140‑1 (MQ=255)                                    gggCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCtt                                                                                                    >  1:896524/1‑140 (MQ=255)                                      gCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAAGATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTg                                                                                                  >  2:690413/1‑140 (MQ=255)                                       cAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTgg                                                                                                 <  1:734891/140‑1 (MQ=255)                                         aCGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGgc                                                                                               >  1:889609/1‑140 (MQ=255)                                         aCGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGgc                                                                                               <  2:457403/140‑1 (MQ=255)                                          cGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAAc                                                                                                        >  2:98292/1‑130 (MQ=255)                                          cGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAAc                                                                                                        <  1:98292/130‑1 (MQ=255)                                               ccccGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAg                                                                                         >  2:845364/1‑140 (MQ=255)                                                  cgcgGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAgtgt                                                                                      >  2:429446/1‑140 (MQ=255)                                                  cgcgGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAgtgt                                                                                      >  2:795448/1‑140 (MQ=255)                                                    cgGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTAt                                                                                    >  2:92143/1‑140 (MQ=255)                                                     ggCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTAtt                                                                                   >  1:774433/1‑140 (MQ=255)                                                               gttgttACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAAt                                                                         <  2:774433/140‑1 (MQ=255)                                                                   ttACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATcgac                                                                     <  1:795448/140‑1 (MQ=255)                                                                       gCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACgatg                                                                 >  1:578355/1‑140 (MQ=255)                                                                        ctctAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACgatga                                                                <  2:896524/140‑1 (MQ=255)                                                                             aGGTTGCCATTATTACTTAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCAcc                                                           >  2:441688/1‑140 (MQ=255)                                                                                  gCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGa                                                      >  1:19371/1‑140 (MQ=255)                                                                                   ccATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGAt                                                     >  2:947588/1‑140 (MQ=255)                                                                                     attattACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATaa                                                   >  1:586064/1‑140 (MQ=255)                                                                                       tattaCTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACg                                                 >  1:617955/1‑140 (MQ=255)                                                                                        attaCTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGt                                                <  2:19371/140‑1 (MQ=255)                                                                                               aaCATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTg                                         >  2:19557/1‑140 (MQ=255)                                                                                               aaCATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTg                                         <  2:617955/140‑1 (MQ=255)                                                                                               aaCATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTg                                         <  2:178593/140‑1 (MQ=255)                                                                                                 cATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCa                                       >  2:93343/1‑140 (MQ=255)                                                                                                 cATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCa                                       >  2:667989/1‑140 (MQ=255)                                                                                                    ctcCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAAt                                    <  2:262070/140‑1 (MQ=255)                                                                                                       cATTTCCGCTGTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACg                                 >  2:724382/1‑140 (MQ=255)                                                                                                                 gTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTg                                         >  2:591833/1‑122 (MQ=255)                                                                                                                 gTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTg                                         <  1:591833/122‑1 (MQ=255)                                                                                                                 gTCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCAt                       >  2:336678/1‑140 (MQ=255)                                                                                                                  tCCATGTTGTCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATg                      <  1:19557/140‑1 (MQ=255)                                                                                                                      tgttgtCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGcgcc                  >  2:359559/1‑140 (MQ=255)                                                                                                                        ttgtCATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGcgccgc                <  1:336678/140‑1 (MQ=255)                                                                                                                            cATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGa            >  1:684962/1‑140 (MQ=255)                                                                                                                            cATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGa            >  2:355702/1‑140 (MQ=255)                                                                                                                            cATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGa            >  2:14727/1‑140 (MQ=255)                                                                                                                            cATGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGa            >  1:134746/1‑140 (MQ=255)                                                                                                                             aTGGTTCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGAt           >  1:692613/1‑140 (MQ=255)                                                                                                                                 ttCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTca       <  2:853652/140‑1 (MQ=255)                                                                                                                                 ttCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTca       <  2:14552/140‑1 (MQ=255)                                                                                                                                  tCACAGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTcac      <  2:586064/140‑1 (MQ=255)                                                                                                                                     caGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACaaa   >  1:658503/1‑140 (MQ=255)                                                                                                                                     caGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACaaa   >  1:470802/1‑140 (MQ=255)                                                                                                                                      aGTGCCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACaaaa  <  1:359559/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         ACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAA  >  NC_000913/1051597‑1051869

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