Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,578,548	G→A	Q264* (CAG→TAG)  Q264* (CAG→TAG)	mcrB ← mcrB ←	5‑methylcytosine‑specific restriction enzyme subunit McrB 5‑methylcytosine‑specific restriction enzyme subunit McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,578,548	0	G	A	100.0%	87.2 / NA	27	Q264* (CAG→TAG)  Q264* (CAG→TAG)	mcrB mcrB	5‑methylcytosine‑specific restriction enzyme subunit McrB 5‑methylcytosine‑specific restriction enzyme subunit McrB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/13);  total (14/13)

TAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTA  >  NC_000913/4578307‑4578787                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                gagacaacaaCGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGATTTTCATCGTTTTCGGATTAGGTTAGGGGAACAGACCAGTTTTCCCTTCGTTTATCATTTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCAGCTATAATAAAAATATACTTTTTCGCTGGCTGCTCTTTAGCTTGCTAACAAAAAtt                                                                                                                                                                                                                                        <  2:156478/247‑1 (MQ=255)               cAGAGAGCGATCGGCCGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGt                                                                                                                                                                                                                          >  2:35379/1‑251 (MQ=255)                agagagCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTc                                                                                                                                                                                                                         <  1:64691/251‑1 (MQ=255)                                  ttCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGAc                                                                                                                                                                                                       >  1:150410/1‑251 (MQ=255)                                  ctCATTAAACCGATGATATAAACATTCTCAGGGAAATAGAATAGCCTTTCATTGCTTTGGGTGTAGGTTAGGGGCACAGACCAGTACTCATCTCGTTTATCATGCTCCTTTAATATCATCACTTCGCTAAATACTTTAAAGTGATTGGAAAGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTAGACGCTCTTTAGCTTGCTAATAAAAATTGTAAAATATGCTGTATTTAGGTCGGAAGCCGAc                                                                                                                                                                                                       >  1:150413/2‑251 (MQ=255)                                        aaaCCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACAAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAACTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                                                                                                                                 >  1:403826/1‑251 (MQ=255)                                         aaCCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTATATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGAACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTc                                                                                                                                                                                                >  1:216736/1‑251 (MQ=255)                                                     aaaCATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGccc                                                                                                                                                                                    >  2:425207/1‑251 (MQ=255)                                                       aCATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGATATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTg                                                                                                                                                                                  <  2:3981/251‑1 (MQ=255)                                                           tatcCGTGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGGCCAGTTTTCACCTCGTTTATCATGTGCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAAGATTCCGTCTTTACGTCGGAATCCGACGCCATTCGGACGATAGCCCTGtat                                                                                                                                                                               <  1:189313/248‑1 (MQ=255)                                                               cGGGACATAGAATCGTTCTTCATCGTTTTCTGAGTAGGTTAGGGGAACAGACCAGATTTCACCTCTTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCTCGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAt                                                                                                                                                                          <  2:150413/251‑1 (MQ=255)                                                               cGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAt                                                                                                                                                                          <  2:150410/251‑1 (MQ=255)                                                                gggACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTGTAGCTTGCTAACAAAAATTGTAAAATATGTCGTATTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATc                                                                                                                                                                         <  1:210416/251‑1 (MQ=255)                                                                                       gTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATgg                                                                                                                                                   >  2:246307/1‑250 (MQ=255)                                                                                         tttCGGAGTAGGTTAGGGGCAAAGACCAGTTTTCACCTCGTTTATCCTGTTCCATTAACATCATCACTTCGCCCAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGCCGCCATTCGGAAGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAAc                                                                                                                                                >  2:403513/1‑251 (MQ=255)                                                                                            cGGAGTTGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGATCCATTAACATCATCACTTCGCCAAATACTATACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTCAAATATGCCGTCTTTACGTCGGCAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATACGATTGTTGGAACTg                                                                                                                                              >  1:65259/1‑250 (MQ=255)                                                                                                                     gTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACg                                                                                                                                                                                                                   <  2:17654/156‑1 (MQ=255)                                                                                                                     gTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACg                                                                                                                                                                                                                   >  1:17654/1‑156 (MQ=255)                                                                                                                               cGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGt                                                                                                          >  2:464709/1‑251 (MQ=255)                                                                                                                                                      catcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTg                                                                                   <  1:425207/251‑1 (MQ=255)                                                                                                                                                          acgtCGCCAAATACTTTACTGAGATTGTCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGTTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGTAACTGAACCATATTGACGCGTTTCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGcaac                                                                               <  1:185027/248‑1 (MQ=255)                                                                                                                                                               gCCAAATCCTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTATTTCGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAgg                                                                          >  2:58299/1‑251 (MQ=255)                                                                                                                                                                    aTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGttttt                                                                     <  2:237147/251‑1 (MQ=255)                                                                                                                                                                         ttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCa                                                                  <  1:100008/249‑1 (MQ=255)                                                                                                                                                                                       aCGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATAGTTACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTGTCCAACTCCGGGCGGCCCCt                                                  <  2:46245/251‑1 (MQ=255)                                                                                                                                                                                                                         ttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCACGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGttt                >  2:84634/1‑251 (MQ=255)                                                                                                                                                                                                                                       tttAGCTTGCTAACAACAATTGTAAAATATGCTGCCTTTACGTCGGAACCCGACGCCACTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGTCGGCGTGCACCAACGGTTTTTCCAACGCCGTGCGGCCCCTGGAGAATAATATTTTTTTTGATTGTTAGTCGTTTGAGAATCggcgctt  >  2:348708/1‑245 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                TAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTA  >  NC_000913/4578307‑4578787

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