Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,568,721	A→G	I145T (ATC→ACC)	mcrB ←	5‑methylcytosine‑specific restriction endonuclease subunit McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,568,721	0	A	G	100.0%	82.6 / NA	23	I145T (ATC→ACC)	mcrB	5‑methylcytosine‑specific restriction endonuclease subunit McrB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/13);  total (10/13)

GCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTC  >  NZ_CP009273/4568637‑4568802                                                                                     |                                                                                  gCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtgt                                                                             <  1:503313/91‑1 (MQ=255)  cATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtgtt                                                                            >  1:143166/1‑91 (MQ=255)  cATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtgtt                                                                            >  1:342585/1‑91 (MQ=255)          aGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCa                                                                    >  2:215982/1‑91 (MQ=255)               gTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATagag                                                               >  2:868903/1‑91 (MQ=255)               gTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATagag                                                               >  1:260621/1‑91 (MQ=255)                      tCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATc                                                        <  1:312487/91‑1 (MQ=255)                       cAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCg                                                       >  1:100327/1‑91 (MQ=255)                           ttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGgtg                                                   <  2:710333/91‑1 (MQ=255)                           ttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGgtg                                                   <  2:273670/91‑1 (MQ=255)                             ttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTa                                                 <  1:365887/91‑1 (MQ=255)                             ttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTa                                                 <  2:119746/91‑1 (MQ=255)                                                           aaTTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACa                   <  2:100327/91‑1 (MQ=255)                                                              tAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGc                <  1:599857/91‑1 (MQ=255)                                                              tAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGc                <  1:528353/91‑1 (MQ=255)                                                                aaTATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGt              <  1:528101/91‑1 (MQ=255)                                                                aaTATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGt              <  1:435947/91‑1 (MQ=255)                                                                    ttAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATa          >  1:34148/1‑91 (MQ=255)                                                                           ataGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATa          >  1:572748/1‑84 (MQ=255)                                                                           ataGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATa          <  2:572748/84‑1 (MQ=255)                                                                            taGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTc  >  2:197956/1‑91 (MQ=255)                                                                                cGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGc                >  1:672647/1‑73 (MQ=255)                                                                                cGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGc                <  2:672647/73‑1 (MQ=255)                                                                                     |                                                                                  GCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTC  >  NZ_CP009273/4568637‑4568802

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