Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,563,567	T→C	I13V (ATT→GTT)	yjiK ←	PF06977 family protein YjiK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,563,567	0	T	C	100.0%	56.3 / NA	20	I13V (ATT→GTT)	yjiK	PF06977 family protein YjiK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/9);  total (11/9)

AAATCACCATTGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCC  >  NC_000913/4563332‑4563808                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  aaaTCAGCATGGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGAAATTGCTAAAAAGTTTATTGCTTTGGGCCGACCAGGTCAGCGAAGAGAAATTATTCGTTATGCCTGCAATTTCTTTGACGTCAATGGTGTCGTGGGAATTCTGAAATGATGCTGCATGATTGCTTTTGCGGGCGCAAGATTGCAAGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGtttcct                                                                                                                                                                                                                                    >  2:206120/1‑248 (MQ=255)          ttGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACACTGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAAt                                                                                                                                                                                                                           >  2:71941/1‑251 (MQ=255)                        aacaaTGGCGGCAGGATTGATACTAGTGCAACAAAGCGTATTGCTTTGTGCCGACCAGGCCAGCGACGAGATATTATTCGTTATGACTGCACTTTCTTTGCCGCCAATGGTGGCGTGGTAATTCTTAAATGATGCTGCCTGATTGCTTTTACGGGCGCAAGATTGCACGCAGAAAGTACACACTGACACCATGGCAAAAAGAATCACAATCACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACgg                                                                                                                                                                                                             >  2:241741/1‑251 (MQ=255)                                 ggcAGGTTTGTTGATAGTGCTAAAAAGAGTCTTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTAGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCACGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCtt                                                                                                                                                                                                    >  2:116144/1‑251 (MQ=255)                                                      aaaaaGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGata                                                                                                                                                                                <  1:251321/250‑1 (MQ=255)                                                                       gtgCCGACCAGGTCAGCGAAGAGCTATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGAAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTACGGCAAATCAATAAGATATCGAGAAGTATACGTAGc                                                                                                                                                              >  2:39429/1‑251 (MQ=255)                                                                          cctaCCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCGCGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGAc                                                                                                                                                           <  1:22352/248‑1 (MQ=255)                                                                                     ggCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATc                                                                                                                                                <  1:35640/250‑1 (MQ=255)                                                                                                   attCGTTATGCCTGCAATTTCTTTGCCGTCTATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGTCGATGCGTTTACTTAACTAAATACTTTTTGTCACGTTTATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCCTAACGCGgaga                                                                                                                                   <  2:105995/250‑1 (MQ=255)                                                                                                       gTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGAGCAAGATTGCACGAAGAAAGTACCCACTGCCACACTGGCAAAAAGAATCAAAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCCCGGTAATCCCTTTAAGGCAAATCAATACGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAAccc                                                                                                                              >  1:155905/1‑251 (MQ=255)                                                                                                                             cGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGAAAGATTGCACGCAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACACCGTAATGTGCAATGCAGCGccc                                                                                                        >  2:174158/1‑250 (MQ=255)                                                                                                                                 aaTGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAACGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTAACGGAAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGac                                                                                                     >  1:238249/1‑250 (MQ=255)                                                                                                                                    ggtggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAAACAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTCAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGAAAcc                                                                                                  <  2:55317/250‑1 (MQ=255)                                                                                                                                          gTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCctgc                                                                                           >  2:18471/1‑251 (MQ=255)                                                                                                                                                     gAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAACAAAt                                                                                <  1:39429/251‑1 (MQ=255)                                                                                                                                                             gctgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATGGACtcag                                                                        <  1:192596/251‑1 (MQ=255)                                                                                                                                                                           tttACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCAACCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGAt                                                           <  2:41810/250‑1 (MQ=255)                                                                                                                                                                             tACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGAtaac                                                        >  1:64033/1‑248 (MQ=255)                                                                                                                                                                                    gcAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAACGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGt                                                 >  2:252154/1‑251 (MQ=255)                                                                                                                                                                                                                                   acacaATAACGGAGATGCGTTTACTGAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAAACCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAgcc  <  2:96200/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  AAATCACCATTGGTGGTCATTTCAACAATGGCGGCAGGTTTGTTGATAGTGCTAAAAAGAGTATTGCTTTGTGCCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCC  >  NC_000913/4563332‑4563808

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