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%	97.8 / NA	33	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 (16/17);  total (16/17)

GCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTA  >  NC_000913/4578304‑4578766                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                           gCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATCTAAACATTCTCCGGGACATAGCATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCAATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATCGATTTCATCTATAAAAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaa                                                                                                                                                                                                                       >  1:300262/1‑250 (MQ=255)                 ccAGAGATCGATCGGCAGTATTCATTAAACCGATGATATAAAAATTCTCCGGGACATAGAATCGTTATTCATCGTTTTCGGAGTAGGTTAGGGGCACAGACCAGTTTTCACCTCGTTTATCATGTTCTATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTTTAAAATATGCCg                                                                                                                                                                                                      >  1:70736/1‑251 (MQ=255)                                  gtggtggccgggccGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCGCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGcc                                                                                                                                                                                     <  1:292190/239‑1 (MQ=255)                                       cATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCCGTTTTCACCTCGTTTATCATGTTCAATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCATGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGaa                                                                                                                                                                                        >  1:452155/1‑243 (MQ=255)                                                                               cGTTCTTCATCGTTTTCGGAGTCGGTTAGGGGAAAAGACCAGTTTTCACCTCGCTTCTCCTGTTCCATTAACATCCTCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTCGCTTGCTAACAAAACTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGAAGATAGCCCTGTATAAAATCCTCATAGCtat                                                                                                                                         >  1:382835/1‑250 (MQ=255)                                                                                    ttcATCGTTTTCGGAGTAGGTTAGGGGAACAGTCCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTGTCTCTGGGTGCTCTGTAGCTTGCTAACAAAAATTGTAAAATATGCCGGCTTTACGGCGGAAGCCGACGCCTTTCGGACGATAGCCCTGTATAAAAGCCTCATAGc                                                                                                                                            <  1:192677/242‑1 (MQ=255)                                                                                      cATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCTGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTg                                                                                                                                  <  1:297370/250‑1 (MQ=255)                                                                                         cGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTTTATAAAAGCCTCATAGCTATAAGATTGATg                                                                                                                               <  1:350866/250‑1 (MQ=255)                                                                                                ggAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCCCTTAGCAAAATACTTTACTGAGATTGGCACGATTGCTTTAATCTATACTATAAATATCCTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTAc                                                                                                                                                                                               >  1:254532/1‑179 (MQ=255)                                                                                                      ggTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCatat                                                                                                                 >  1:4354/1‑251 (MQ=255)                                                                                                       gTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATAATCACTTAGACAAATAAATTACTGAGATTGGCACGCTTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATAtt                                                                                                                >  1:307296/1‑251 (MQ=255)                                                                                                              gAACANNCCNNNTTNNACCTCGTTTATCATGTTCCATTANCATCATCACTTCGCCCAATACTTTACTGANANTNGCACNNTTGNTTTCATCTATNATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAANTNTAAAATATGCCGTCTTTACGTCGGAAGCCGANGCCATTCGGACGATAGCCCTGTATNAAATNNTCNtnnnnnnnn                                                                                                                                       >  2:333767/1‑213 (MQ=255)                                                                                                                 cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAGTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGc                                                                                                      <  1:100509/251‑1 (MQ=255)                                                                                                                 cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTgn                                                                                                      <  1:116109/251‑2 (MQ=255)                                                                                                                   gACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCAAAATACTTTACTGAGATTGGCCCGATTGATTTCATCTATAATAAACATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCgg                                                                                                    >  1:361453/1‑251 (MQ=255)                                                                                                                     ccAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGt                                                                                                                                                                                             <  1:176023/160‑1 (MQ=255)                                                                                                                        nttttCACCTCGTTTATCATGTTCCATTAACATCATCAATTACCCAAATACTTTCCTGAGATTGGCACGATTGACTTCATCTCTAATAAAAATATACTTTTTCTCAGCCTGCTCTTTAGCTTGCTAACAAACATTGTAAAATATGCCATCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATACTCATAGCTATAAGATTGATGGAACTGAAACATATTTACGAGTTGCGGAGCCt                                                                                               >  1:87973/2‑251 (MQ=255)                                                                                                                        gTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTAATGAGATTGGCACGATTGCTTTCATCTCTAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCt                                                                                               >  1:179079/1‑251 (MQ=255)                                                                                                                                                     aCATCATCACTGCGCCAAATACTTTACTGATATTGGCACGATTGATTTCATCTATAATAAAAATATCCGTTTTCTCTGGCTGCGCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTGACGTCGGAAGCCGACGCCAGTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGAGTGATGGAACTGAACCGTATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACgg                                                                  <  1:320130/251‑1 (MQ=255)                                                                                                                                                       annntnncgnnnncNNNNNNNNTANNNAGATTGNCNNGATTGNTTTNATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATANCCCNGTATNAANTCCTCATAGCTATAANATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTNNTNNNANCAAGTAAg                                                                          <  2:307296/228‑1 (MQ=255)                                                                                                                                                             aCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATCTTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGcaac                                                          >  1:31223/1‑251 (MQ=255)                                                                                                                                                              nttCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAACATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAAAAAAAATTg                                                                                                                                                                                                                  >  1:38634/2‑98 (MQ=255)                                                                                                                                                                         aCTTTACTGAGATTGTCACGATTGATTTCATCTATAATAAAAATATTCTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGTCGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTc                                               <  1:69057/250‑1 (MQ=255)                                                                                                                                                                              aCTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACg                                          <  1:203838/250‑1 (MQ=255)                                                                                                                                                                                 tatattttcACGATTGATTTCATCTATAATATAAATATACGTGTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTaa                                                                                                                                                                                                               <  1:394974/74‑1 (MQ=255)                                                                                                                                                                                                gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAAGATGCCGGCTTTACTTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTGGCTGAGCCTTTTCTCCTGTGAGCAATTAAGCGAGACGGCGTGGAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTggagga                       <  1:333767/251‑1 (MQ=255)                                                                                                                                                                                                    tCATCTATAATAAAAATATACTTTTTCGCTGGCTGCTCTTTAGCGTGCTAACAAAGATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGcn                                                                                                                                                                               <  1:92493/95‑2 (MQ=255)                                                                                                                                                                                                       tCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGcc                                                                                                                                                                                     <  2:282215/86‑1 (MQ=255)                                                                                                                                                                                                       ncTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGACGTCTTTACGTCGGAAGcc                                                                                                                                                                                     >  1:282215/2‑86 (MQ=255)                                                                                                                                                                                                           taataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGCCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGCTGGAAATGAACCATATTGACGCGTTGAGGAGCATTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAgg                                                     >  1:196718/1‑210 (MQ=255)                                                                                                                                                                                                                      taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATACGATTGATGGAACTGAACCATATTGACGCGTTGCCGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAAAAAAGGTTTTTGAAACGACGGTCGGCCCCTGGAGGATAATATTTTTTTTGATggata  >  1:8812/1‑247 (MQ=255)                                                                                                                                                                                                                                     tgctCTTTAGCTTGCGAACAAAAAGTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGACTGATGGAACTGAACCATAGTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTGTTTCCAACGCCGggcggc                                 <  1:320785/204‑1 (MQ=255)                                                                                                                                                                                                                                              gCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGc                                                                                 <  1:375722/147‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                           GCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTA  >  NC_000913/4578304‑4578766

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