Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,164,123	C→G	A162G (GCA→GGA)	btuB →	cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,164,123	0	C	G	100.0%	110.2 / NA	31	A162G (GCA→GGA)	btuB	cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (17/14);  total (17/14)

GTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGG  >  NC_000913/4163880‑4164357                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           gTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                                                                                                                                                                                     >  1:125903/1‑251 (MQ=255)     gCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATAGCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATAGCATCACGACGCGCGATGAGCCCGTAACGGAAATTTCAGGAGGGTGGGGaa                                                                                                                                                                                                                                 <  2:27695/251‑1 (MQ=255)           aTCTATTTTTATTCGCGGTACCAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGGAGCAAt                                                                                                                                                                                                                            >  1:91361/1‑250 (MQ=255)               aTTTTTAATCGCGGTACAAATGCCAGTCATGTGTTGGTGTTACTGGATGGCGTACGCCTGAATCTGGCGGGGGGGAGTGGTTCTGACGACCTTAGCCATTTCCCTATTGCGCTTGTCCAGCGTGTTGAATAGATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGAGATGAACCCGGAACGGAAATTTCAGGATGGTGGGGaatctata                                                                                                                                                                                                                           <  2:227863/247‑7 (MQ=255)                          cgGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                                                                                                                                                                            <  2:8447/251‑1 (MQ=255)                              aCAAAAGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGGAGGGTGGGGAAGCAATAGCTATCAGAAATATGATGt                                                                                                                                                                                                        >  2:122839/1‑251 (MQ=255)                                               gtTGGTGTTAATTGATGGCGTACGCATGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACTCTGCaacaac                                                                                                                                                                                       <  1:203102/251‑1 (MQ=255)                                                 tGGTGTTCATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                                                                                                                                                                     >  1:232974/1‑251 (MQ=255)                                                                gcttACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGATCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAGAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                                                                                                                                                                      <  1:24801/248‑1 (MQ=255)                                                                    aCGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                                                                                                                                                                   <  1:136335/250‑1 (MQ=255)                                                                          gAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt                                                                                                                                                            >  2:119166/1‑251 (MQ=255)                                                                                   gggggTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGCATATATCCGTGGGCCGAGCTCCGCTGTTTAAGGTTCCGATGCAATAGGCGGGGTGGTGAATATCCTCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGAGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt                                                                                                                                                   >  1:166938/1‑251 (MQ=255)                                                                                           gTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAt                                                                                                                                           >  2:141086/1‑251 (MQ=255)                                                                                              gTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACt                                                                                                                                        <  1:10921/251‑1 (MQ=255)                                                                                                   gCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCa                                                                                                                                      <  1:232664/248‑1 (MQ=255)                                                                                                         cTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatg                                                                                                                              >  2:196339/1‑250 (MQ=255)                                                                                                                   tCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGc                                                                                                                    <  2:71698/250‑1 (MQ=255)                                                                                                                             gcTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATAATCATGGTTATGATGTTGTTGCCTATGGGAAt                                                                                                          >  1:68475/1‑250 (MQ=255)                                                                                                                               ttGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATAcc                                                                                                       >  1:217431/1‑251 (MQ=255)                                                                                                                                        gtgtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCaa                                                                                              <  1:41124/251‑1 (MQ=255)                                                                                                                                                     ccGTGGCCCGCGCTCAGCTGTATATGGCTCAGATGCAATAGGCGGGGAGGAGAATCTCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAAGCAACAACAACTGGGGGATAAGACACGGGTAAAGCTGTTGGGCGATTATGCCCATACCCCTGGTTATGATGTTGTTGCCTATGGTAATTACGGCACGCAAGCGCAGACAGAt                                                                                  >  2:108893/1‑250 (MQ=255)                                                                                                                                                                 ctcagCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACTACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGCTTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGttgtt                                                                     >  2:6028/5‑248 (MQ=255)                                                                                                                                                                    cgcggTTTATGGTGCCGAGGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAAGTTCAGGAGGGTGCGGAAGCAATAGTTATCATATCTATGATGTCTCTACGCAGCAACAACTGGGCGATAAGACTCGTGTAGTGGAGTTGGGCGATTGTGTCCAGGCTCAGTGTTATGATGTGGTGGCCTATGGTAAGGCCCGATAGCAAGCCCAGACAGATAACTAGGTTTTTTTaa                                                                  <  2:152233/247‑1 (MQ=255)                                                                                                                                                                    cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTCGTGAATATCATCACGCCGCGAGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAAATGGGGGATAAGACACGGGTCACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGCTGTTGTTGCCTATGGAAATACCGGAACGCAAGCGCAGACAGATACAGATGGTTTTTTa                                                                   >  2:257017/1‑250 (MQ=255)                                                                                                                                                                                 tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAAGTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACTCttt                                                      <  2:154547/250‑1 (MQ=255)                                                                                                                                                                                  ccGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATGTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTAt                                                    <  2:206218/251‑1 (MQ=255)                                                                                                                                                                                   cGATGCAATAGGCGGGGTGGTGAATATCATCCCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGACGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATg                                                   >  2:136987/1‑251 (MQ=255)                                                                                                                                                                                         aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGgcgcgc                                             >  1:137392/1‑251 (MQ=255)                                                                                                                                                                                         aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATAAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATAACGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGgcgcgc                                             >  1:137384/1‑251 (MQ=255)                                                                                                                                                                                                                 caCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCt                     >  1:164281/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggaacggaaATTTCAGGAGGGAGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATTGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCTAGCGCAGACAGATAACTATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCgg  <  2:136/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           GTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGG  >  NC_000913/4163880‑4164357

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