Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,215,832	(C)5→4	coding (645/732 nt)	yehW ←	glycine betaine ABC transporter membrane subunit YehW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,215,828	0	C	.	100.0%	117.5 / NA	28	coding (649/732 nt)	yehW	glycine betaine ABC transporter membrane subunit YehW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/15);  total (13/15)

TCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCA  >  NC_000913/2215691‑2215962                                                                                                                                          |                                                                                                                                       tcctgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGccc                                                                                                                                      <  1:147867/140‑1 (MQ=255) tcctgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGccc                                                                                                                                      <  2:129183/140‑1 (MQ=255)    tgttTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGcccctg                                                                                                                                   <  1:360293/140‑3 (MQ=255)          tGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGAt                                                                                                                               >  2:56262/1‑137 (MQ=255)          tGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGAt                                                                                                                               <  1:56262/137‑1 (MQ=255)          tGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGAt                                                                                                                               <  1:177380/137‑1 (MQ=255)          tGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGAt                                                                                                                               >  2:177380/1‑137 (MQ=255)               ttAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATAcg                                                                                                                       <  2:319225/140‑1 (MQ=255)               ttAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATAcg                                                                                                                       <  1:54547/140‑1 (MQ=255)                tAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATAcgc                                                                                                                      <  1:190555/140‑1 (MQ=255)                tAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATAcgc                                                                                                                      >  2:180163/1‑140 (MQ=255)                    cGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTa                                                                                                                  <  2:164779/140‑1 (MQ=255)                                 gTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTaa                                                                                                     >  2:288214/1‑140 (MQ=255)                                      gCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCg                                                                                                >  1:294839/1‑140 (MQ=255)                                             ggTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatgatga                                                                                         <  1:288214/140‑1 (MQ=255)                                             ggTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCgatgatga                                                                                         <  2:114568/140‑1 (MQ=255)                                                 ataCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATggg                                                                                     <  1:304319/140‑1 (MQ=255)                                                 ataCCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATggg                                                                                     >  2:304319/1‑140 (MQ=255)                                                                  gCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTAGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTgg                                                                    >  2:144445/1‑140 (MQ=255)                                                                    tGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGcc                                                                  >  2:48229/1‑140 (MQ=255)                                                                      gCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGcccc                                                                <  1:180163/140‑1 (MQ=255)                                                                                            tCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGcc                                          >  1:79618/1‑140 (MQ=255)                                                                                                  aGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCa                                    <  2:79618/140‑1 (MQ=255)                                                                                                     cGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATg                                 >  2:392851/1‑140 (MQ=255)                                                                                                        tCTGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTg                              <  1:127640/140‑1 (MQ=255)                                                                                                          tGCGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGAt                            >  2:183224/1‑140 (MQ=255)                                                                                                            cGATGATCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATaa                          <  2:232632/140‑1 (MQ=255)                                                                                                                 atCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCAcc                     >  2:29027/1‑140 (MQ=255)                                                                                                                  tCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTAcc                                     <  2:182264/123‑1 (MQ=255)                                                                                                                  tCGCCGCCAGTGCCACCAGTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTAcc                                     >  1:182264/1‑123 (MQ=255)                                                                                                                                    gTAACG‑CCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCa  <  2:294839/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       TCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCTTTATTTTGCGTGCTGGCTAAGCGCCTGCACCAGCCTTTCAAACAGGCGGTCTGCGATGATCGCCGCCAGTGCCACCAGTAACGCCCCCTGGATCACATACGCGGTATTAAATCCGCTAAGCCCGATGATGATGGGCGTACCCAGCGTGCTGGCCCCTACCGTTGAGGCGATCGTCGCCGTACCAATGTTGATAATCACCGAAGTTCGCACGCCCGCCA  >  NC_000913/2215691‑2215962

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