Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,070,052	C→T	S117N (AGC→AAC)	yihR ←	putative sulfoquinovose mutarotase YihR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,070,052	0	C	T	100.0%	48.6 / NA	16	S117N (AGC→AAC)	yihR	putative sulfoquinovose mutarotase YihR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/11);  total (5/11)

TCGCGTGCTCATCGACGGCATAAACCTGGTTTGCCGGAAGCTGGAACAGATATTCATCAACAGATGTTAAGTTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTG  >  NC_000913/4069807‑4070269                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                          tCGCGTGCTCATCGACGGCATAAACCTGGTTTGCCGGAAGCTGGAACAGATATTCATCAACAGATGTTAAGTTGCAGGTTAAATAAGGATGAATCCCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTtggt                                                                                                                                                                                                                       <  2:67422/250‑1 (MQ=255)                            ggTTTGCCGGAAGCTGGAACAGATATTCATCAACAGATGTTAAGTTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGTGGTTAATTCACTAATTTGCCAGt                                                                                                                                                                                             <  1:190104/249‑1 (MQ=255)                             gTTTGCCGGAAGCAGGAACAGATATTCATCAACAGATGTTAAGTGGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGtctc                                                                                                                                                                                          <  2:161964/251‑1 (MQ=255)                                                  atatTCATCAACAGATGTTAAGTTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGaa                                                                                                                                                                     <  1:96195/251‑1 (MQ=255)                                                   tatTCATCAACAGATGTTAAGTTGCAGGTTAAATAAGGATGAATACCCAAGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGGTAATTCACCAATTTGCCAGGCTCGCCAGGCAAAGAGCCCGTGAAt                                                                                                                                                                    >  1:143395/1‑251 (MQ=255)                                                                                 aaaTAAGGATGAATACCCAAGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAACCAGCCCGTGAATAGCGGCCTTTGAACTATGTTCATTAATggg                                                                                                                                      >  1:195172/1‑251 (MQ=255)                                                                                                     gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCTGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTGACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCTTTAATGGGTAACTGATATTCCTGccccg                                                                                                                  <  2:157430/251‑2 (MQ=255)                                                                                                                  gCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTACCTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTaaa                                                                                                      >  1:257565/1‑249 (MQ=255)                                                                                                                                                    cGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGTTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGc                                                                   <  2:182234/251‑1 (MQ=255)                                                                                                                                                     ggAAAGTCCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCa                                                                  <  1:83715/251‑1 (MQ=255)                                                                                                                                                                                            tgAAGCCAGAATAAACGGAAAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGGGAATAACGGCTTTTGAAATTTTTTACTTAATGGGTAACTGATATTCCTTCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCAAGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTACTCGGGCTtgt                           >  2:41836/1‑251 (MQ=255)                                                                                                                                                                                            ggAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTtgt                           <  1:131053/250‑1 (MQ=255)                                                                                                                                                                                                  cAGCATAAACGGATAGCCGTAACTGGGTGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCCCTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGa                     <  2:160310/251‑1 (MQ=255)                                                                                                                                                                                                             gATAGCCTGAACTGGGCGGCTGAAAAGCTTTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCATCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATCCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGTCCACGGAATAAGCACTTTACCCAGATGCGCCAGCTGCATCTCGTCTGGTTTGTGCGGGATGACTAAATGa          <  2:223310/251‑1 (MQ=255)                                                                                                                                                                                                              aTAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGAc         <  1:173673/251‑1 (MQ=255)                                                                                                                                                                                                                     tAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGTTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACCGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCCGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTg  >  2:23826/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                          TCGCGTGCTCATCGACGGCATAAACCTGGTTTGCCGGAAGCTGGAACAGATATTCATCAACAGATGTTAAGTTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTGGCGGTTAATTCACTAATTTGCCAGTCTCGCCAGGCAAGCAGCCCGTGAATAGCGGCTTTTGAACTATGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTG  >  NC_000913/4069807‑4070269

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