Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,101,728	(G)7→8	coding (39/984 nt)	kdgT →	2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,101,721	1	.	G	100.0%	88.0 / NA	26	P11R (CCG→CGG)	kdgT	2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (17/9);  total (17/9)

CTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCG  >  NC_000913/4101587‑4101845                                                                                                                                        |                                                                                                                              cTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggg                                                                                                                         <  1:162124/140‑1 (MQ=255)      aTTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggat                                                                                                                    <  2:243044/140‑4 (MQ=255)       ttGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑ggggggggatg                                                                                                                   >  1:479879/1‑136 (MQ=255)            tAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCt                                                                                                               >  2:413131/1‑140 (MQ=255)                                tttCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGgcg                                                                                           >  2:32242/1‑140 (MQ=255)                                  tCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCa                                                                                         <  2:55332/140‑1 (MQ=255)                                  tCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCa                                                                                         >  1:219934/1‑140 (MQ=255)                                             tGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCacac                                                                              >  1:461341/1‑140 (MQ=255)                                               aTGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCt                                                                            >  2:37227/1‑140 (MQ=255)                                                    ccATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCg                                                                       >  1:307878/1‑140 (MQ=255)                                                     cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCtt                                                                           >  1:50864/1‑135 (MQ=255)                                                     cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCtt                                                                           <  2:50864/135‑1 (MQ=255)                                                                  cGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGaa                                                         <  2:219934/140‑1 (MQ=255)                                                                       cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAtt                                                    >  2:411902/1‑140 (MQ=255)                                                                       cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAtt                                                    >  1:450779/1‑140 (MQ=255)                                                                       cAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATAtt                                                    >  2:194878/1‑140 (MQ=255)                                                                          cccTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTg                                                 >  1:384128/1‑140 (MQ=255)                                                                                      tGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCa                                     <  1:418572/140‑1 (MQ=255)                                                                                      tGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCa                                     >  2:354981/1‑140 (MQ=255)                                                                                           aCAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACggg                                >  2:12807/1‑140 (MQ=255)                                                                                                        aTGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTa                   <  2:384128/140‑1 (MQ=255)                                                                                                         tGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTAc                  >  2:104436/1‑140 (MQ=255)                                                                                                         tGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTAc                  >  2:505771/1‑140 (MQ=255)                                                                                                         tGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTAc                  >  1:315401/1‑140 (MQ=255)                                                                                                         tGCAGATAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTAc                  <  1:194878/140‑1 (MQ=255)                                                                                                               tAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGcc            >  2:180264/1‑140 (MQ=255)                                                                                                               tAAAACGCTCGATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGcc            <  1:293672/140‑1 (MQ=255)                                                                                                                         gATTGAGAAAATCCCG‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCg  <  1:413131/140‑1 (MQ=255)                                                                                                                          aTTGAGAAAATCCCGGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCg  <  1:404889/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                              CTTTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCC‑‑GGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCCCATTCTGGCG  >  NC_000913/4101587‑4101845

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: ‑