Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,237	C→T	C188Y (TGC→TAC)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,237	0	C	T	100.0%	72.9 / NA	24	C188Y (TGC→TAC)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/10);  total (14/10)

TGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCAC  >  NC_000913/4307100‑4307359                                                                                                                                          |                                                                                                                           tGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTaa                                                                                                                          >  2:682012/1‑140 (MQ=255)      gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGt                                                                                                                     >  1:14490/1‑140 (MQ=255)      gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGt                                                                                                                     >  1:635979/1‑140 (MQ=255)                 ccGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGGTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTCCCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCAcc                                                                                                          <  2:174366/140‑1 (MQ=255)                   gTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCa                                                                                                        >  1:169652/1‑140 (MQ=255)                         cgTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTg                                                                                                  >  2:646935/1‑140 (MQ=255)                           ttttCAAGCAGACTCTGGACGAAAGGGGCGTTTT‑CGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCca                                                                                               <  1:455273/140‑1 (MQ=255)                                  gCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                         <  1:646935/140‑1 (MQ=255)                                  gCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                         <  1:682012/140‑1 (MQ=255)                                  gCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                         <  2:845121/140‑1 (MQ=255)                                  gCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                         <  2:635979/140‑1 (MQ=255)                                      aCTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAgt                                                                                     >  1:459409/1‑140 (MQ=255)                                        tctGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAgtgt                                                                                   <  2:671803/140‑1 (MQ=255)                                                aaaGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCat                                                                           >  1:22476/1‑140 (MQ=255)                                                  aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                         <  1:610800/124‑1 (MQ=255)                                                  aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACAcgcg                                                                                         >  2:610800/1‑124 (MQ=255)                                                      gCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTcc                                                                     >  1:406935/1‑140 (MQ=255)                                                         ttttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAg                                                                  <  2:22476/140‑1 (MQ=255)                                                                cATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGata                                                           <  2:707499/140‑1 (MQ=255)                                                                        aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAg                                                   >  1:315261/1‑140 (MQ=255)                                                                          aaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGtt                                                 >  2:867446/1‑140 (MQ=255)                                                                               cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCt                                            >  1:76932/1‑140 (MQ=255)                                                                                     aTGGGTATTTTGGGGGTTGTTTTTCGTCCCAGCGTCTTAGCGCCATTCCAGAGTATTTGGTTTCCAGGCACCCAGGTTTGCCACACGCGCGGTGTGGGGTCATATCCCCCGGGGGGATATGACCCAGTTCGCCTGCcaca                                      <  2:14490/140‑1 (MQ=255)                                                                                                  gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCg                         >  2:74186/1‑140 (MQ=255)                                                                                                                       tCTTAGCGCCATTCCAGAGTAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc    >  1:693635/1‑140 (MQ=255)                                                                                                                         ttAGCGCCATTCCAGAGTATTTGGTTCCCAGCCCCCCAGGATTCCCACCCGCGCAGTGTGGGGTAATACCCCCGGGGGGTATATGCCCCAGTCCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac  <  2:246684/140‑1 (MQ=255)                                                                                                                         ttAGCGCCATTCCAGAGTAATTGGTTCCCAGGCCCCCAGGTTTCCCCCCCGCGCGGTGTGGGGTCATATCCCCCGGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac  <  2:49847/140‑1 (MQ=255)                                                                                                                         ttACCGCCTTTCCAGAGTATTTGGTTTCCAGGCCCCCAGGATCGCCCCCCGCGCGGTGTGGGGTCATATCCCCGGGGGGGATATCACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac  <  2:406935/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                           TGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCAC  >  NC_000913/4307100‑4307359

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