Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	409,021	T→C	R25R (CGT→CGC)	ykiD →	protein YkiD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	409,021	0	T	C	100.0%	111.3 / NA	32	R25R (CGT→CGC)	ykiD	protein YkiD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/22);  total (10/22)

GAAAACCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGTAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAATCAGGTTTTGAATTAAC  >  NC_000913/408887‑409156                                                                                                                                       |                                                                                                                                        gAAAACCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGa                                                                                                                                    <  2:541741/140‑1 (MQ=255) gAAAACCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGa                                                                                                                                    <  1:405427/140‑1 (MQ=255)  aaaaCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaa                                                                                                                                   <  2:660943/140‑1 (MQ=255)  aaaaCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaa                                                                                                                                   <  1:402266/140‑1 (MQ=255)  aaaaCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaa                                                                                                                                   <  1:432734/140‑1 (MQ=255)   aaaCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaat                                                                                                                                  <  2:189785/140‑1 (MQ=255)    aaCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaata                                                                                                                                 <  2:9021/140‑1 (MQ=255)     aCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaa                                                                                                                                   <  1:524121/137‑1 (MQ=255)     aCCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGaa                                                                                                                                   >  2:524121/1‑137 (MQ=255)        ggggCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATta                                                                                                                             <  1:462554/140‑1 (MQ=255)              aCGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGcc                                                                                                                       <  2:731603/140‑1 (MQ=255)                             aGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCACTAAAAATTATTGCTCTCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTa                                                                                                        <  2:24293/140‑1 (MQ=255)                                           ccGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGAcg                                                                                          >  2:692719/1‑140 (MQ=255)                                              gCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCgg                                                                                       <  1:280443/140‑1 (MQ=255)                                                            aaTGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATcc                                                                         <  1:13715/140‑1 (MQ=255)                                                            aaTGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATcc                                                                         <  1:423262/140‑1 (MQ=255)                                                                    aaCGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTaa                                                                 >  2:588543/1‑140 (MQ=255)                                                                      cGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAg                                                               >  1:224687/1‑140 (MQ=255)                                                                            tGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAAt                                                         <  2:224687/140‑1 (MQ=255)                                                                            tGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAAt                                                         <  2:664248/140‑1 (MQ=255)                                                                                         aGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGt                                            >  2:613937/1‑140 (MQ=255)                                                                                               aaaCGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCt                                      <  1:588543/140‑1 (MQ=255)                                                                                                             gctgcGCTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCACATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATcc                        >  2:513587/1‑140 (MQ=255)                                                                                                                 cgcTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCg                                   <  2:678685/125‑1 (MQ=255)                                                                                                                 cgcTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCg                                   >  1:678685/1‑125 (MQ=255)                                                                                                                 cgcTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCa                    >  2:44105/1‑140 (MQ=255)                                                                                                                   cTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAAt                  >  2:633741/1‑140 (MQ=255)                                                                                                                   cTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAAt                  >  1:38536/1‑140 (MQ=255)                                                                                                                   cTAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAAt                  <  1:44105/140‑1 (MQ=255)                                                                                                                    tAAAAATTATTGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCACATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAATc                 <  1:513587/140‑1 (MQ=255)                                                                                                                          ttattGCTCGCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAATCAGGttt           <  1:692719/140‑1 (MQ=255)                                                                                                                                   gCCGCAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAATCAGGTTTTGAATTAAc  <  1:653097/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        GAAAACCGGGGCAACGCGGAAAATTCGCAGTGAAGCGCGTCACCGGCTTGAGGAAAAACAATGACACAACGTCCATGGAGCAAGTTACAGCGTAAAACGCATAATATTGCTGCGCTAAAAATTATTGCTCGCCGTAGCGAATAATTATATGCCTGGTGTGGCTTCGTACGCCGGATAAGACGCGGCAGGCGTCGCATCCGGCATTAAAGGAAAATCAGCAATTAACGTTGTGCTTCGCCACCTAATCCTTCAATCAGGTTTTGAATTAAC  >  NC_000913/408887‑409156

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