Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,424,407	C→T	intergenic (+21/+109)	ulaF → / ← yjfY	L‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 domain‑containing protein YjfY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,424,407	0	C	T	100.0%	55.7 / NA	20	intergenic (+21/+109)	ulaF/yjfY	L‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 domain‑containing protein YjfY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/8);  total (12/8)

GCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTAC  >  NC_000913/4424170‑4424556                                                                                                                                                                                                                                              |                                                                                                                                                      gCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTAGCCTGGGGGTAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGCTGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCggg                                                                                                                                             >  1:238358/1‑248 (MQ=255)  naTACGCAGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCggg                                                                                                                                             >  1:84023/2‑247 (MQ=255)                          cAGCACGGGCCGTTCGCCTGGGGGAAAGCTGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAACATGGCGCGGATTGGCCGCGGCATTAACCCACGAATCAATCACATCGACAGCTTCCTGATGAATAAACCCTCCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCAGCAGCCGtt                                                                                                                 >  1:326576/1‑251 (MQ=255)                              aCGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTc                                                                                                             <  1:304595/251‑1 (MQ=255)                                 ggCCGTTCGCCTGGGGGAAAGATGCCCACGATGCGGGCCATAACGCGGTGGCGATGGAAGACGTGGCGAAACTGGCGTGGATTGCTCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGCATAACCAATTCATTCGTTAACACGGTCCTAACGCTTATTACGGGCAGCAGTAGAACACGCGCTGCGGAAACTTTCTTCCTAGGGAGACACATGGTCTACTTCCGCCGCCGTTTTTCaa                                                                                                           >  1:218159/1‑250 (MQ=255)                                         gCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGAGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACCGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGCTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGcc                                                                                                  >  1:223558/1‑251 (MQ=255)                                              ggggAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTg                                                                                             >  1:113175/1‑251 (MQ=255)                                                       tGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTAGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTGAAGCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCAAACGCTTATTACGGGCATAAGGAGAGCACGCGCGGCGGAAATGTTCTTCCTCGGGAGATAACTGGTCTAAGTCCGCCGCCGTCTTTCAAAAAAAAGCCCCCTGCGAAgggg                                                                                     <  1:145657/250‑1 (MQ=255)                                                                                     gATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCCCAACTCAATCACATCGCCAGCTTCCTGATGAATAAACACTTCATGCGTCAACACGGTCCTAACGCTTCTTATGGGCCGAAGTAGCACACGCGCCGCGGAAATTTTCTTCCTCGGGAGATCACTGGACTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                      >  1:307886/1‑251 (MQ=255)                                                                                                 ggCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGc                                          >  1:9246/1‑251 (MQ=255)                                                                                                   cGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAAATCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCt                                         >  1:112440/1‑250 (MQ=255)                                                                                                        aTGGCGTGAATTGCCCGCGGCAGTAACCCACAACTCAATCACATCGACAGCATCCTGATGAATAAACACTTTATGCCGAAACACGGTCCTAACGCTTATTACGGGCAGAACTAAAACCCGCGCTGCGGAAATTCTCTTCCGCGGGATAGAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGATgggg                                                                                     <  1:220933/201‑1 (MQ=255)                                                                                                                gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATCAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAg                                                                                                    >  1:134223/1‑178 (MQ=255)                                                                                                                    gCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCgcgc                       >  1:165448/1‑251 (MQ=255)                                                                                                                                    cacaACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGCCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAACTTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGc                                                                                                   >  1:273676/1‑159 (MQ=255)                                                                                                                                                                              tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAAGTTTCGTCCTCGGGATATAACTGGTCTAATTCCGCAGCCGTTTTGCAAAAAAAAGCCCCCGGCGAAGGGGGGAAAGCAAACTATGCCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGc               <  1:44941/201‑1 (MQ=255)                                                                                                                                                                               cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTAc  <  1:40965/213‑1 (MQ=255)                                                                                                                                                                               cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTTCTTCCGCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTAc  <  1:40974/213‑1 (MQ=255)                                                                                                                                                                                                                                  cTNCGGAANTTTTNTTCCTCGGGAGATAANTGGTCTAAGTCCGCAGCCGNNTNTCAAANAAAAg                                                                                                    <  2:134223/64‑1 (MQ=255)                                                                                                                                                                                                                                      ggAAATTTTCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCTGTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                                       <  1:52134/105‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                      GCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTAC  >  NC_000913/4424170‑4424556

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