Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,426,929	C→T	C165Y (TGC→TAC)	yjfZ ←	DUF2686 domain‑containing protein YjfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,426,929	0	C	T	100.0%	83.9 / NA	26	C165Y (TGC→TAC)	yjfZ	DUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

CATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTT  >  NC_000913/4426700‑4427151                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               aatacGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAAACCGAGTAGAGTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAAGCATGGGGGAAAGCAGAATACGCGTTTAGCAATGCAGCAATGATGAAATCGAGGGCAGTAGAAACGCCTGTAACAGTGCGAAGTTGGTTGATATAAAGATCGGTTGTTGGTAGTTCAGATGACCATTAAATAGGCATTTGGGCAGGaa                                                                                                                                                                                                           <  2:168530/250‑1 (MQ=255)     cGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGACACCCAAACCCGAGTAGAGTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGAAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAACGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGt                                                                                                                                                                                                       >  2:353787/1‑251 (MQ=255)       aGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTa                                                                                                                                                                                                     <  1:353787/251‑1 (MQ=255)             aCTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGt                                                                                                                                                                                               >  2:334659/1‑251 (MQ=255)                   gTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCt                                                                                                                                                                                          >  1:18713/1‑250 (MQ=255)                         cTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGttt                                                                                                                                                                                   >  2:222113/1‑251 (MQ=255)                                       gctcgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGACCGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGAAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACAAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCt                                                                                                                                                                     >  2:144964/1‑251 (MQ=255)                                                                  cGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAAcc                                                                                                                                          >  1:85763/1‑251 (MQ=255)                                                                    gCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCAc                                                                                                                                        <  1:334659/251‑1 (MQ=255)                                                                             gATAAAGAGCGCACACCACTTTTCAAATCTTGGGGTAAATCAGAATACGCGATTAGCAATGACGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGCTCGATTGTTGGACGTTAAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATACTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAACACACTTCgtgtat                                                                                                                               >  2:354179/1‑251 (MQ=255)                                                                                   gAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtattgtgt                                                                                                                          >  1:116428/1‑250 (MQ=255)                                                                                        cAGACCAATTTTCAAATCATGGGGTAAAACAGAATACGCGTTTAGCAATGCAGCAATGATGAAGCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTGGGGCGATATAAAGATCGATTGGTGGTAGTTCAGATGACCAGTAATTAGGCATTGGGGCAGGAAGGGTTACAACGTATCCTATCGTTGGGTCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAATCGACTTCGTGTAGTGTGTATTTAc                                                                                                                    <  2:145501/251‑1 (MQ=255)                                                                                                                      aGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGc                                                                                       >  1:76532/1‑250 (MQ=255)                                                                                                                      aGAATACGCGTTTAGAAATGCAGCAATGATGAATCCGAGGTCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGAACAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTAGAGTCTGGAGTACCATAACCAGGTCACAAAACCACTTCGTGTAAGGTGTATTTACTGTGGGGTTCGAAGGATGCGACATATTGCa                                                                                      <  2:283817/251‑1 (MQ=255)                                                                                                                                    gCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGt                                                                        <  1:96222/251‑1 (MQ=255)                                                                                                                                              aaTGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGTGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGAGTCGACATCTTGCAATAAATATTTCAGCTCTGCGGTTa                                                              <  2:87786/251‑1 (MQ=255)                                                                                                                                               atgatgAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaat                                                                                    >  1:318767/1‑228 (MQ=255)                                                                                                                                               atgatgAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaat                                                                                    <  2:318767/228‑1 (MQ=255)                                                                                                                                                 gatgaATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAAcc                                                           >  1:198595/1‑251 (MQ=255)                                                                                                                                                  atgaATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATAGAAAGATCGATTGTTGGTAGGTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTGACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTGCTGTTGGGTTCGAAGGATTCGACAGCTTGCAATAAATATTTCAGTTCTGCGGTTaaact                                                          <  2:93573/251‑4 (MQ=255)                                                                                                                                                      aTCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCTGTATAGTCTGGAGTACCATGACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCa                                                      <  1:352234/251‑1 (MQ=255)                                                                                                                                                              gCATTAGAAACGCCTGTAACAGGGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTa                                              <  2:186603/251‑1 (MQ=255)                                                                                                                                                                   aGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCtt                                         <  1:169773/251‑1 (MQ=255)                                                                                                                                                                             gTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgta                                                                                                                                >  2:221587/1‑154 (MQ=255)                                                                                                                                                                             gTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgta                                                                                                                                <  1:221587/154‑1 (MQ=255)                                                                                                                                                                                                          aTCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCGTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCtt  <  2:18713/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               CATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTT  >  NC_000913/4426700‑4427151

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