Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,563,778	(G)6→5	coding (148/258 nt)	ytiA ←	uncharacterized protein YtiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,563,778	0	G	.	100.0%	138.1 / NA	31	coding (148/258 nt)	ytiA	uncharacterized protein YtiA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (16/15);  total (16/15)

TGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCA  >  NC_000913/4563538‑4564025                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        tGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGcc                                                                                                                                                                                                                                              >  2:7501/1‑251 (MQ=255)       aaTGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGc                                                                                                                                                                                                                                        <  2:96119/251‑1 (MQ=255)                                 agaTGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCa                                                                                                                                                                                                              <  2:139135/251‑1 (MQ=255)                                                 cGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCCAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGCATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCCTTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAACTGCATAACCCGCAGCCGCCGCGTACATCACGGGGAATCTCGCTg                                                                                                                                                                                               >  2:15187/1‑250 (MQ=255)                                                  gAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAGCCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACTGGGAATATCGCTGAc                                                                                                                                                                                             <  2:366496/251‑1 (MQ=255)                                                          ttttGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaaa                                                                                                                                                                                      <  1:7501/250‑1 (MQ=255)                                                                cACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATg                                                                                                                                                                                <  2:9931/250‑1 (MQ=255)                                                                       aTCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGCGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgc                                                                                                                                                                        >  2:239114/1‑251 (MQ=255)                                                                          ccTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCTTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcctg                                                                                                                                                                     <  2:177413/251‑1 (MQ=255)                                                                                          aTAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGtt                                                                                                                                                     >  2:282038/1‑251 (MQ=255)                                                                                                           ataCGTAGCGATCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGGCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCg                                                                                                                                     <  2:334482/250‑1 (MQ=255)                                                                                                                    gACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaa                                                                                                                            >  1:266447/1‑250 (MQ=255)                                                                                                                      ccTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaata                                                                                                                         <  1:282038/251‑1 (MQ=255)                                                                                                                                 ccTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAAc                                                                                                              <  2:305421/251‑1 (MQ=255)                                                                                                                                                aCCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCgtg                                                                                               >  1:141856/1‑251 (MQ=255)                                                                                                                                                            gCAATGCAGCGCCAGACGCCATCCTGCTGTCAATAAATTCACTCAGTCAGCCCAGTGATCACCCATTGTATTTCCTCGGCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTTCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAAGAAAAATAAAAACTCTTTCAGCGGCGTGTCCTCACTCGtt                                                                                   <  1:272353/251‑1 (MQ=255)                                                                                                                                                                       ccAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCCCTCGTTGACCAATTAcc                                                                        >  2:152635/1‑251 (MQ=255)                                                                                                                                                                        cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAccg                                                                       >  1:253734/1‑251 (MQ=255)                                                                                                                                                                         aGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAccgc                                                                      >  1:124584/1‑251 (MQ=255)                                                                                                                                                                             acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGCCCAATTACCGCagcg                                                                  >  2:79137/1‑251 (MQ=255)                                                                                                                                                                              ccATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGACCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGACAGCCACATGCATAACCCGCAGCCGCCGCGCAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGCTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACcgccgccc                                                                 >  2:217234/1‑249 (MQ=255)                                                                                                                                                                                     gctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAAtt                                                            <  2:5532/249‑1 (MQ=255)                                                                                                                                                                                      ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGa                                                                                                                                 <  2:111257/179‑1 (MQ=255)                                                                                                                                                                                      ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGa                                                                                                                                 >  1:111257/1‑179 (MQ=255)                                                                                                                                                                                                                cAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATgggg                               >  1:294726/1‑250 (MQ=255)                                                                                                                                                                                                                              tGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGCGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCCCTCGTTGCCCAATTACCGCCGCGCCAATTGCGCCCGCAAATTGCGCATCAGGATGGGTATTTACCGGccttt                 >  2:236095/1‑247 (MQ=255)                                                                                                                                                                                                                              tGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAAAGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAACGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTAACTCGTTGACCAATTACCGCCGCGCCCATTCAGCCAGCAAAATGCGCATCAGGATGGGTATTTAACGGCATTc                 >  1:319776/1‑251 (MQ=255)                                                                                                                                                                                                                              ggTATTGCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCGGACGCGAAAAGTCATGCCCGCGGCCTGCAACGCTTGGCGGGTTTGTATGACGCCGACCGTGGAGTGAAAGAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCAtt                  <  1:236095/249‑1 (MQ=255)                                                                                                                                                                                                                                tATTTCCGCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCCCGGGGAATATCGCTGCCGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGCCCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGTCGTGTCCTCACTCTTTGACCAATTACCGCCTCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCg                <  1:220380/250‑1 (MQ=255)                                                                                                                                                                                                                                    tCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTTCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAggg           >  2:36301/1‑250 (MQ=255)                                                                                                                                                                                                                                     ccTCGCCGGGGA‑ACAGGTCACCCAAATGCATCACCCGCAGCCGCCGCTTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGATTTGCGGTTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTTTCCTCCCTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGt           <  1:321566/250‑1 (MQ=255)                                                                                                                                                                                                                                              ggggACAGGTCAGCCAAATGCCTACCCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGATCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCa  <  1:251266/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        TGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCA  >  NC_000913/4563538‑4564025

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