Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,022,002	C→T	E628K (GAG→AAG)	yccS ←	putative transporter YccS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,022,002	0	C	T	100.0%	92.0 / NA	27	E628K (GAG→AAG)	yccS	putative transporter YccS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/8);  total (19/8)

TGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTCCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTT  >  NC_000913/1021761‑1022236                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           tGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGa                                                                                                                                                                                                                                     <  1:274291/249‑1 (MQ=255)          tGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAggg                                                                                                                                                                                                                           <  1:127103/250‑1 (MQ=255)          tGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAggg                                                                                                                                                                                                                           >  2:127103/1‑250 (MQ=255)                      cgacCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCt                                                                                                                                                                                                              <  2:311730/251‑1 (MQ=255)                               tcAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGt                                                                                                                                                                                                     >  2:347687/1‑251 (MQ=255)                                      gTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGtt                                                                                                                                                                                               >  2:196860/1‑250 (MQ=255)                                                       ccAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGaa                                                                                                                                                                              >  1:84832/1‑250 (MQ=255)                                                                    cGACCAGGGGCTCTTTGCTGTCTGCCCGCGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTATTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAACCGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTcgcg                                                                                                                                                                 >  2:259274/1‑250 (MQ=255)                                                                                    gCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGt                                                                                                                                                >  2:32794/1‑251 (MQ=255)                                                                                                                            aGGCTAGCTAATGCCTCCTTGAAGCGTTCTTCATCAGACGGTTGATGAAGTAACGCGTCATCAGCATAGCAAACTGCGTCATCAAGAAAAGCAAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGAGTTAGCTGGTAAACGTATGGTTAAGGCACAGCAACAGAAACGCGGCTTCGCGAATTTGCGGGGTCACGATCGACTCGCAGGACatcgttgacaacaacacggccagccc                                                                                                        >  2:388296/1‑227 (MQ=255)                                                                                                                            aGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACGGCGTCATCAAGAAACGCCAGAATTTCAGTATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGTTCGCTGGACATATTTGAGACCTCCGACGCCAGCTc                                                                                                        <  2:380287/251‑1 (MQ=255)                                                                                                                               ctagctaATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACCCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGc                                                                                                     >  2:170364/1‑251 (MQ=255)                                                                                                                                 agctaATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCAt                                                                                                   <  2:378344/251‑1 (MQ=255)                                                                                                                                     aaTGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACg                                                                                               >  2:153211/1‑251 (MQ=255)                                                                                                                                                 aCGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATc                                                                                   <  1:180830/251‑1 (MQ=255)                                                                                                                                                     gTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCgcgg                                                                               <  2:319827/251‑1 (MQ=255)                                                                                                                                                             tcaGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAat                                                                       >  1:368487/1‑251 (MQ=255)                                                                                                                                                                               aaCGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGACTTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGCTGTGTGCATCCCGGCGGGCAATACGATACGccctacggct                                                     >  2:376/1‑243 (MQ=255)                                                                                                                                                                                 cgcgTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTAt                                                   >  1:411004/1‑251 (MQ=255)                                                                                                                                                                                                    aCTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTg                                >  1:136464/1‑251 (MQ=255)                                                                                                                                                                                                      tGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCt                              >  1:62224/1‑251 (MQ=255)                                                                                                                                                                                                         gTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAAAGCGGCTTCGCGAATTTGCGGGGCAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTccc                           >  1:26882/1‑250 (MQ=255)                                                                                                                                                                                                           catcaAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCGTGATGGTATTGCTCCAGt                         <  2:62224/251‑1 (MQ=255)                                                                                                                                                                                                               aaGAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTATCGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATgg                     >  2:289626/1‑251 (MQ=255)                                                                                                                                                                                                                 gAAACGCCAGAATTTCAGGATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGc                    >  1:260101/1‑250 (MQ=255)                                                                                                                                                                                                                              ttCAGGCTTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGAc      >  1:135478/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggATTAGTTAACTGCTTCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGtt  >  1:214028/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           TGAGTAATCTGGCGTTGCAGACGACCAATCTCAGGCAGTAATGCAATCAATAATCCAACTTGTTGTACGACCAGAGGTTCTTTGCTGTCTGCCCGTGGTTCAAGTTGCTGCATCCGCTGTTTCAGGCTAGCTAATGCCTCATTGACGCGTTCTTCATCAGCAGGTTGATGATGTAACGCGTCATCAACATAGCAAACTGCGTCATCAAGAAACGCCAGAATTTCAGGATTAGTTAACTGCTCCCGGTGAGCACCGAGGGCTGAGATATAGCTGGTAAACGTATGGTTAAGGCACAGCAACCGAAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTT  >  NC_000913/1021761‑1022236

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