Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,578,759	G→A	I193I (ATC→ATT)  I193I (ATC→ATT)	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,759	0	G	A	96.8%	89.3 / NA	32	I193I (ATC→ATT)  I193I (ATC→ATT)	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/1);  major base A (15/15);  minor base C (1/0);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATAT  >  NC_000913/4578518‑4579002                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     aCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAAt                                                                                                                                                                                                                                             >  1:11127/1‑251 (MQ=255)   tttttCACTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGCCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGCACTGAACCATCTTGACGCGTTGCGGAGACTTTTCTCCTGTCAGAAAGTCAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGTCGGCCCATGGAGACTAATATTTTTTTTAATGGTTAATCg                                                                                                                                                                                                                                           >  2:360918/1‑251 (MQ=255)          ctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAACATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGCCGATATCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCCACAAAGGTTTTTCCAACACCGGGCGGCCCATGGAGGATAATATTTTTTTTAATTGTTAATAGTTTGatt                                                                                                                                                                                                                                    >  2:255069/1‑249 (MQ=255)           tGGCTGCTCTTTAGCTCGCTGACAAAAATTGTCCACAATGCCGTCTTTACGTCGGCACCCGACGCCATTCGGGAGAACGACATGTATAACATCCTCATAGCTATAAGATTGCTGGAACTGAACTATATTGACGCGTTGCGGAGCATTTTCTCCTGTCAGCAAGTAAGACAGCCGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCAGCCCCTGGCGGATAATATTTTTTTTCATGGTTAGTCGTTTGAGt                                                                                                                                                                                                                                    >  2:448799/1‑250 (MQ=255)                 ctctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGtct                                                                                                                                                                                                                             >  1:420038/1‑251 (MQ=255)                                 cAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAgg                                                                                                                                                                                                             >  1:372692/1‑251 (MQ=255)                                    aaaTTGTAAAAAATGCCGACCGTACGTCGGAAGCCGTCGCCATTCGGACGATAGCCCTGTATAAAAACCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGATCCTTGTCGCCTGTCAGCAAGGAAGCCCGACGGCGGTCAACAAAGGTTTTTCCAACGCCGGGTGGCCCCGGGAGGATAATATTTTGTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGAt                                                                                                                                                                                                          <  2:368843/251‑1 (MQ=255)                                                 tGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCAttt                                                                                                                                                                                             >  1:180859/1‑251 (MQ=255)                                                  gccttCTTTACGTCTGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGATCCTGTTCTCCTGTCAGCAAGTAAGCCAGACTGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTa                                                                                                                                                                                            <  1:278842/247‑1 (MQ=255)                                                                     cagtcGCCATTCTGACGATATCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTGTTCTCCTGTCAGCAAGTAAGCCAGACGTCGTGCAACAATGGTTTTTCCAACGCCGGGCGTCCCCTGGAGTATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGTTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGActtta                                                                                                                                                                         <  2:100379/247‑5 (MQ=255)                                                                            cATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATAGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAATCAAATCATTTATCGCATCTTCCAGACAGTATGATTc                                                                                                                                                                   >  2:4744/1‑250 (MQ=255)                                                                                                ataAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGaa                                                                                                                                              <  1:74442/250‑1 (MQ=255)                                                                                                           aTAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGc                                                                                                                                  >  1:145274/1‑251 (MQ=255)                                                                                                            tAGCTATAAGATTGATGGAACTGATCCATATTGCCGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGcc                                                                                                                                 >  2:174756/1‑251 (MQ=255)                                                                                                             aGCTATAAGATTGATGGAACTGAACCATTTTGACGCGTTGCGGAGCTTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGcc                                                                                                                                 <  2:11127/250‑1 (MQ=255)                                                                                                                tataAGATTGATGGAACTGAACCATATTGACGCGTTGAGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGCTAATATTTTTTTTAATGTTTAATCGTTTGAGTATCGTCTCTATTTAGGTTTCAGTGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGCTTCA‑GTTTTTGACATAGGTGGAATAACACTATTGCCAGa                                                                                                                              >  2:258790/1‑250 (MQ=255)                                                                                                                 ataAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATAGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAAt                                                                                                                            <  1:4744/251‑1 (MQ=255)                                                                                                                                    aCCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTc                                                                                                         <  2:145274/251‑1 (MQ=255)                                                                                                                                            gACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGAtt                                                                                                  >  1:438841/1‑250 (MQ=255)                                                                                                                                                cgTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTAtgtt                                                                                             <  1:309224/251‑1 (MQ=255)                                                                                                                                                                       cAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCggg                                                                      >  1:425754/1‑251 (MQ=255)                                                                                                                                                                       cAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCgg                                                                       <  1:255069/250‑1 (MQ=255)                                                                                                                                                                                              caacaaAGGTTTTTCCAATGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCATGTTTTT‑ACATAGGTGGAATAACACTCTTGTCAGAATTAAAGATTAAGTTATAGTCGTTGATTATGTGGTCCAGCATATATGCAAATCGATTGTAATCAATTCCCTATGAGACt                                               <  2:172928/251‑1 (MQ=255)                                                                                                                                                                                                            ccAACGCCGGACGGCCCATGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCGCTATTGTGGTTTCAGGGATCAAAAACTCATTTACCGCATCTTCCAGACAGTATGATTCC‑GTTTTTGACATAGGTGGACTAACACTCTTGCCAGAATTAAATATTAATTTATAGTAGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGTGACTTTTTGGGAAAAGGc                                 >  2:34062/1‑251 (MQ=255)                                                                                                                                                                                                                    gggCGGCCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGCTTATGTTGTCCAGCATAGAGGCAACTCGGGTGTAATCCAAACCCTGTTATACTTTTTGGGAAAAGGCGAAAt                            >  2:380208/1‑248 (MQ=255)                                                                                                                                                                                                                       cggcCCCTGGAGGATAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATCCCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACt                      <  1:116320/251‑1 (MQ=255)                                                                                                                                                                                                                       cggcCCCTGGAGGATAAGATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCACTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAACCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACt                      <  1:431617/251‑1 (MQ=255)                                                                                                                                                                                                                               tgaggaTAATATTTTTTTTAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACAct                                                                                                                                       <  1:411409/129‑1 (MQ=255)                                                                                                                                                                                                                               ggagtaTAATATTGTTGTTAATGGTTAATTGTTTGAGTATCGTTTCTATGGTGGTTTCAGGGCTAACCCAATCATTTAACGCATCTTCCATACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCCGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGTGTGTAATCAATACCCTGTGATACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTAt              <  2:346857/251‑1 (MQ=255)                                                                                                                                                                                                                               ggaggaTACTATTTTTTTTAATGGTGAATCTTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAAACATTTAACGAATCTTCCAGACATTATTATTCA‑TTTTTTGACATAGGTGGAATAACTTTATTTCCAGAATTAAATATGAATTTATAGTCGATTATTATGTTGTCCAGCCTAGATGCAAATCGGGTGTAATCAATACCCTGTGAGTCTTTGTGGGATCAGGGGTAATAGGACTGTCCga                <  2:305625/249‑2 (MQ=255)                                                                                                                                                                                                                                 aggaTAAAATTGTTTTTAATGGTTAATCGTTTGAGGATCGTATCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTGATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCTTAATAGGTCGGTCCGTAttt            <  2:331925/251‑1 (MQ=255)                                                                                                                                                                                                                                           tttttttAATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGatat  >  1:93601/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGCCGGGCGGCCCCTGGAGGATAATATTTTTTTTGATGGTTAATCGTTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCA‑GTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATAT  >  NC_000913/4578518‑4579002

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