Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,578,409	A→G	L310P (CTA→CCA)	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,409	0	A	G	100.0%	123.1 / NA	36	L310P (CTA→CCA)	mcrB	5‑methylcytosine‑specific restriction enzyme subunit McrB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/16);  total (20/16)

ATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTT  >  NC_000913/4578276‑4578538                                                                                                                                      |                                                                                                                                  atCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGTCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGaa                                                                                                                             <  1:488079/140‑1 (MQ=255)     aTGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCAATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGa                                                                                                                         >  2:564379/1‑140 (MQ=255)          agaaaATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGt                                                                                                                    <  1:151320/140‑1 (MQ=255)          agaaaATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGt                                                                                                                    >  2:765234/1‑140 (MQ=255)               aTCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCa                                                                                                               >  2:65385/1‑140 (MQ=255)                     tGCGTAGGGCATAGTCAACAACGGCCAGATAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGt                                                                                                         <  1:528472/140‑1 (MQ=255)                     tGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGt                                                                                                         <  1:541444/140‑1 (MQ=255)                            ggCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCAt                                                                                                  >  1:268846/1‑140 (MQ=255)                             gCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATg                                                                                                 <  2:178360/140‑1 (MQ=255)                               aTAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGtt                                                                                               <  2:553981/140‑1 (MQ=255)                                 aGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTcc                                                                                             >  1:153474/1‑140 (MQ=255)                                               agagagCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCAc                                                                               >  1:709532/1‑140 (MQ=255)                                                gagagCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACt                                                                              <  2:268846/140‑1 (MQ=255)                                                     cGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGcc                                                                         >  2:700869/1‑140 (MQ=255)                                                           gCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATAc                                                                   <  1:65385/140‑1 (MQ=255)                                                             aGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACtt                                                                 >  2:528608/1‑140 (MQ=255)                                                                  tCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTg                                                            >  2:202315/1‑140 (MQ=255)                                                                    aTTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTgag                                                          <  1:167115/140‑1 (MQ=255)                                                                         aCCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCGTTAACATCATCACTTCGCCAAATACTTTACTGAGATTgg                                                     >  1:769243/1‑140 (MQ=255)                                                                                    aaaCATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACgattgatt                                          <  1:564379/140‑1 (MQ=255)                                                                                       cATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCa                                       >  2:433539/1‑140 (MQ=255)                                                                                          tctcCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCAt                                      >  1:586093/1‑138 (MQ=255)                                                                                          tctcCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCAt                                      <  2:586093/138‑1 (MQ=255)                                                                                           cgccGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCta                                   >  1:734556/3‑140 (MQ=255)                                                                                              cGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa                                <  1:765234/140‑1 (MQ=255)                                                                                                ggACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataata                              >  1:318206/1‑140 (MQ=255)                                                                                                 gACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTAtaataa                             >  2:791578/1‑140 (MQ=255)                                                                                                      aGAATCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAata                        >  2:276364/1‑140 (MQ=255)                                                                                                        aaTCGTTCTTCATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAAtata                      >  2:650565/1‑140 (MQ=255)                                                                                                             ttcttcATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACtttt                 >  2:141824/1‑140 (MQ=255)                                                                                                               cttcATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTGACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTtc               >  2:712799/1‑140 (MQ=255)                                                                                                                ttcATCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTtct              <  1:498079/140‑1 (MQ=255)                                                                                                                    tCGTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTgg          <  1:791578/140‑1 (MQ=255)                                                                                                                      gTTTTCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCGTTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGgct        <  2:769243/140‑1 (MQ=255)                                                                                                                       ttttCGGAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGgctg       <  2:734556/140‑1 (MQ=255)                                                                                                                            ggAGTAGGTTGGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCtt  >  1:120836/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                  ATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCATTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTTCGGAGTAGGTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTT  >  NC_000913/4578276‑4578538

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