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%	126.3 / NA	27	coding (148/258 nt)	ytiA	uncharacterized protein YtiA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (13/14);  total (13/14)

ACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAG  >  NC_000913/4563542‑4564014                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             acaATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAAt                                                                                                                                                                                                                               <  1:41197/251‑1 (MQ=255)         aaaaaGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAAcc                                                                                                                                                                                                                       >  1:119238/1‑251 (MQ=255)                                tGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgg                                                                                                                                                                                                >  1:55227/1‑251 (MQ=255)                                              gAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGGATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                                                                                                                                                                  <  2:128964/251‑1 (MQ=255)                                               aaaTACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTACATCACGGGGAATATCGCTGAcg                                                                                                                                                                                 >  2:202708/1‑251 (MQ=255)                                                           tCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATg                                                                                                                                                                     >  2:199790/1‑251 (MQ=255)                                                                   aTCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgc                                                                                                                                                             >  1:248599/1‑251 (MQ=255)                                                                                  aTCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCg                                                                                                                                              >  1:150828/1‑251 (MQ=255)                                                                                               cGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGAc                                                                                                                                  <  1:37528/250‑1 (MQ=255)                                                                                                gagaAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGc                                                                                                                                >  1:189062/1‑251 (MQ=255)                                                                                                 agaAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGcc                                                                                                                               <  1:194552/251‑1 (MQ=255)                                                                                                   aaGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCg                                                                                                                              >  2:259258/1‑250 (MQ=255)                                                                                                         aCGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTg                                                                                                                        <  2:83626/250‑1 (MQ=255)                                                                                                                gACCTTAATGAATCCTTCACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaa                                                                                                                 >  2:257759/1‑250 (MQ=255)                                                                                                                    ttAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCAt                                                                                                                                                                                                                           <  2:141399/140‑1 (MQ=255)                                                                                                                    ttAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCAt                                                                                                                                                                                                                           >  1:141399/1‑140 (MQ=255)                                                                                                                                                        gCAATGCAGCGCCAGACACCATCCTTCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGTTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCTGCGTGTCCTCACTCGtt                                                                        <  2:129920/251‑1 (MQ=255)                                                                                                                                                        gCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTTACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGG‑‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCAGGCCCGCTGCCTGCAACGCGTTGCGGGTTTGTATGACGCGGACCTTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTg                                                                       <  2:133686/251‑1 (MQ=255)                                                                                                                                                        gCAAGGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGtt                                                                        <  1:66485/251‑1 (MQ=255)                                                                                                                                                                      gACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAccgcc                                                          <  1:221749/251‑1 (MQ=255)                                                                                                                                                                            aTCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTCACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCCCTCGTTGACCCATTAccgccccgcc                                                     >  1:176926/1‑250 (MQ=255)                                                                                                                                                                                       aaTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTgg                                                                                                                       <  1:263980/173‑1 (MQ=255)                                                                                                                                                                                       aaTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTgg                                                                                                                       >  2:263980/1‑173 (MQ=255)                                                                                                                                                                                                        agCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGAt                        <  2:119238/251‑1 (MQ=255)                                                                                                                                                                                                          cccAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATAACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTTCGCCCGCAAAATTCGCCTCAGGCTgg                      >  2:261110/1‑251 (MQ=255)                                                                                                                                                                                                              gggATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATACCCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTAt                  <  1:214683/251‑1 (MQ=255)                                                                                                                                                                                                                               ttCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATACCCCGAAGCCCCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAATTCATGCCCGCTGCATGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAACCTCTTTCATCGCCGTGGCCTCCTTCGTTTACCATGTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAg  <  2:96230/250‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             ACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAG  >  NC_000913/4563542‑4564014

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