Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,576,875	T→C	intergenic (+2944/+37)	yjiT → / ← mcrC	putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,576,875	0	T	C	100.0%	80.0 / NA	26	intergenic (+2944/+37)	yjiT/mcrC	putative uncharacterized protein YjiT/Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/16);  total (10/16)

CGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGCGGTATCTACGTGGGGATATATTAATAACCCCCCTATGTTTTCGCCATTTTCAGGCTT  >  NC_000913/4576628‑4577091                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                         cGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCTGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATgcggt                                                                                                                                                                                                                       <  1:273631/251‑1 (MQ=255)    aaaCTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTc                                                                                                                                                                                                                    >  1:122695/1‑251 (MQ=255)      aCTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCgg                                                                                                                                                                                                                  >  1:241140/1‑251 (MQ=255)               gacgtttGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTAt                                                                                                                                                                                                         <  2:304536/247‑1 (MQ=255)                  ttttGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCCGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTAGGGAGCGGTAtttt                                                                                                                                                                                                      >  1:35858/1‑251 (MQ=255)                    ttGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCt                                                                                                                                                                                                    >  1:288975/1‑251 (MQ=255)                            ggTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAAc                                                                                                                                                                                             <  2:203485/250‑1 (MQ=255)                                    aTTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGTGGTATTTTCGAGAAAACTTACCGCt                                                                                                                                                                                     <  2:19986/250‑1 (MQ=255)                                     ttGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCtt                                                                                                                                                                                    >  2:199519/1‑250 (MQ=255)                                        aaGATCCGCTTCGCTCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt                                                                                                                                                                                >  1:77643/1‑251 (MQ=255)                                        aaGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt                                                                                                                                                                                <  2:338961/251‑1 (MQ=255)                                                   cGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCAt                                                                                                                                                                     >  1:356555/1‑251 (MQ=255)                                                   cGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCATCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCAt                                                                                                                                                                     <  2:16353/251‑1 (MQ=255)                                                    gATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATc                                                                                                                                                                    <  1:353904/251‑1 (MQ=255)                                                                 gCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt                                                                                                                                                       <  1:302123/251‑1 (MQ=255)                                                                            aTCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg                                                                                                                                            <  2:233890/251‑1 (MQ=255)                                                                                             aTGGCATTACTTTTATTTCCGCGATGACGTTTCTTTCCAGCTAGTCAAAATTGCTCAGGCCTGCAATGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCCACGGATGTGCTTGTGGAGGTGAGCGCCTCCTGATGATGAGCCGCTCCCGATGCTGTGTCGGGAGCGCTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATTTCGCGTAAGTCTTGATGTATTCACGGCCATTcca                                                                                                                           <  2:175444/251‑2 (MQ=255)                                                                                                        ttGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAAtt                                                                                                                 <  1:85837/250‑1 (MQ=255)                                                                                                              tCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGt                                                                                                           <  2:241140/250‑1 (MQ=255)                                                                                                                                        gggaaaaTGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGTTGGATGTGATCGCCTCCTGATGATAAGCCGCTCTCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAAATTACCGCTTATTTGAAATATTCATCGAAAATGTCGAGTAATGCTAGATGTATACGCGACCATTCCTGACCTAAATTGACGGAACACAAGCCAATAGCGAAGCCATTAAt                                                                                <  2:268567/248‑1 (MQ=255)                                                                                                                                         tCAAAATGGCTCAGGACTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAAGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAAtt                                                                               >  2:179122/1‑251 (MQ=255)                                                                                                                                               tGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTtata                                                                          >  1:299069/1‑250 (MQ=255)                                                                                                                                                           gCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGTTGTGCTGGTGGAGGTGATCGCCTCCGGATGATGAGCCGTTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCa                                                               <  1:227844/249‑1 (MQ=255)                                                                                                                                                                                    cGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGCGGTATCTACGTGGGGATATAtt                                     >  1:319355/1‑250 (MQ=255)                                                                                                                                                                                               gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGCGGTATCTACGTGGGGATATATTAATAACCCCCCt                         <  1:271850/251‑1 (MQ=255)                                                                                                                                                                                                       gATGTGCTGGGGGAGGTGAGCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGTCCATTCCTGTCCTAAATTGACGGTATACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGCGGTATCTACGTGGGGATATATTAATAACCCCCCTATGTTTTc                 <  1:175703/251‑1 (MQ=255)                                                                                                                                                                                                                       tGATCGCCTCCTGATGATGAGCCGCTCCCGATGCGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGTCCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGCGGTATCTACGTGGGGATATATTAATAACCCCCCTATGTTTTCGCCATTTTCAGGCtt  <  1:144271/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                         CGTAAACTGACGGTGACTTTTGAGGATGGTCAGGTATTGAAGATCCGCTTCGATCAGGGGATGGGCTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAAGCCATTAATTTTATAACGATGTTTCACTGCGGTATCTACGTGGGGATATATTAATAACCCCCCTATGTTTTCGCCATTTTCAGGCTT  >  NC_000913/4576628‑4577091

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