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%	80.8 / NA	27	C165Y (TGC→TAC)	yjfZ	DUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/15);  total (12/15)

ACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAA  >  NC_000913/4426699‑4427167                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               acattcGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGa                                                                                                                                                                                                                            <  1:122867/251‑1 (MQ=255)  catacGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGGGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGaa                                                                                                                                                                                                                           <  1:300727/251‑1 (MQ=255)   atacGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGa                                                                                                                                                                                                                            >  2:124924/1‑249 (MQ=255)             aaCTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACg                                                                                                                                                                                                                >  2:55458/1‑251 (MQ=255)                                           cgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATATGGACCGCACCCCCATTTTCAAATCATGGCGTACATCAGAATACGCGTTTCGCAATGCCGCACTGATGAACCCGAGGGCATTAGAAACGCCTGTAACCGTGCGAAGTTGGTCGATATACAGATCGATTGTTGGTAGTTCAGATGCCCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATCGTCTgg                                                                                                                                                                                   >  1:148533/1‑250 (MQ=255)                                                           gTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTc                                                                                                                                                                  <  2:16454/251‑1 (MQ=255)                                                             aGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTccc                                                                                                                                                                >  1:177225/1‑251 (MQ=255)                                                                   cGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAAcc                                                                                                                                                          <  2:303879/251‑1 (MQ=255)                                                                     gCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGta                                                                                                                                                                                           <  1:74299/216‑1 (MQ=255)                                                                     gCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGta                                                                                                                                                                                           >  2:74299/1‑216 (MQ=255)                                                                         aaCTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCg                                                                                                                                                    <  1:112629/251‑1 (MQ=255)                                                                              gATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtat                                                                                                                                               <  1:220679/251‑1 (MQ=255)                                                                                   ggAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtattgtgt                                                                                                                                          <  1:2238/251‑1 (MQ=255)                                                                                               aaTTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTgg                                                                                                                              <  1:52367/251‑1 (MQ=255)                                                                                                   ttCAAATCGTGGGGTAAATAAGAATACGCGTTTAGCAATACAGCAATGATGAATCCGAGAGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTAGTGAATTGTGAATTTACTGTTGGGTTc                                                                                                                          >  2:187653/1‑251 (MQ=255)                                                                                                      aaaTCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGa                                                                                                                        >  1:177148/1‑250 (MQ=255)                                                                                                           aTGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGAtt                                                                                                                  <  2:242444/251‑1 (MQ=255)                                                                                                           aTGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGCCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCTGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGAtt                                                                                                                  >  2:208878/1‑251 (MQ=255)                                                                                                                  aaaTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACAt                                                                                                            <  2:177225/250‑1 (MQ=255)                                                                                                                    aTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCtt                                                                                                         <  1:99504/251‑1 (MQ=255)                                                                                                                                       aaTGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGCTGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGttt                                                                                      >  2:74774/1‑250 (MQ=255)                                                                                                                                                         ccGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGGTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGGTCAGTAATTAGGCATGTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCGTAACCAGGTCCGAAAACCGCTTCGTGTATTGTGTATTTAGTTGTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCAc                                                                     <  2:20382/250‑1 (MQ=255)                                                                                                                                                                                         aaGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGaa                                    <  1:90561/251‑1 (MQ=255)                                                                                                                                                                                             tGGTCGATATAAAGCTCGATTGATGGTAGTTCAGATGACCAGTAATTAGGCACTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATCGTCTGGAGTACCATAACCAGGTCACAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTAGACATCTTGCAATAAATATTTCATTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGcatggaacc                                 >  2:290643/1‑242 (MQ=255)                                                                                                                                                                                                                  gttgGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGt            <  1:191250/250‑1 (MQ=255)                                                                                                                                                                                                                        aggtCAGATTACCAGTAATTAGGCATTTGGGCAGGGAGGGTTACAACGTAGCCTATCGTTTGGCCGGTATAGGCTGGAGTACCATAACCAGGTGCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTGCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTGAACCTCGCACCCGCGTACGCTTTATTGTATTTCTGGTCGTGAATATTATCTAATGTAGCTTCATGGTtccatcc     >  2:288643/4‑251 (MQ=255)                                                                                                                                                                                                                             aGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCaaa  >  2:136137/1‑249 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAA  >  NC_000913/4426699‑4427167

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