Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,210,492	A→G	E49E (GAA→GAG)	mcrA →	e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,210,492	0	A	G	100.0%	137.8 / NA	39	E49E (GAA→GAG)	mcrA	e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/21);  total (18/21)

TTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGCCCCCG  >  NC_000913/1210356‑1210626                                                                                                                                         |                                                                                                                                       ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCgg                                                                                                                                     >  2:405001/1‑140 (MQ=255) ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCgg                                                                                                                                     <  2:147227/140‑1 (MQ=255)    ataataATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCggttg                                                                                                                                  <  2:333701/140‑1 (MQ=255)      aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCggttggt                                                                                                                                <  2:477230/140‑1 (MQ=255)        taatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGgttggttg                                                                                                                              >  2:138194/1‑140 (MQ=255)           tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGAtt                                                                                                                           >  1:152680/1‑140 (MQ=255)           tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGAtt                                                                                                                           >  1:323369/1‑140 (MQ=255)           tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGAtt                                                                                                                           <  1:342626/140‑1 (MQ=255)           tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGAtt                                                                                                                           <  1:455174/140‑1 (MQ=255)                ttGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACGATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGgt                                                                                                                      >  2:388656/1‑140 (MQ=255)                   aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTa                                                                                                                   >  2:181802/1‑140 (MQ=255)                    aCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTAt                                                                                                                  >  1:225980/1‑140 (MQ=255)                     cTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATc                                                                                                                 >  2:528476/1‑140 (MQ=255)                        aaaGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCa                                                                                                              <  1:138194/140‑1 (MQ=255)                          aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAg                                                                                                            <  2:323369/140‑1 (MQ=255)                           gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGt                                                                                                           <  1:438362/140‑1 (MQ=255)                            cTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTc                                                                                                          <  2:152680/140‑1 (MQ=255)                                 tgtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCgtagt                                                                                                     <  1:405001/140‑1 (MQ=255)                                   ttCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGtat                                                                                                   >  2:495426/1‑140 (MQ=255)                                      gATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATc                                                                                                >  2:208751/1‑140 (MQ=255)                                          ggTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGc                                                                                            <  2:253150/140‑1 (MQ=255)                                             gAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGtc                                                                                         >  2:164245/1‑140 (MQ=255)                                                    aTGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGtcatcatca                                                                                  <  1:495426/140‑1 (MQ=255)                                                                 tCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACata                                                                     >  1:499641/1‑140 (MQ=255)                                                                     aTCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCa                                                                 <  1:407818/140‑1 (MQ=255)                                                                                gTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATg                                                      >  1:377616/1‑140 (MQ=255)                                                                                      gCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATgaaagaa                                                >  1:339805/1‑140 (MQ=255)                                                                                        cGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGaaaa                                              >  1:447834/1‑140 (MQ=255)                                                                                               cAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAatccatc                                       <  2:377616/140‑1 (MQ=255)                                                                                                    aCAAAGATTACAATATCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGt                                  <  1:260848/140‑1 (MQ=255)                                                                                                                 tatCGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACtt                     <  1:388656/140‑1 (MQ=255)                                                                                                                    cGCTCTTGATTATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGa                  <  1:181802/140‑1 (MQ=255)                                                                                                                             ttATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGAt                 <  2:207896/132‑1 (MQ=255)                                                                                                                             ttATATCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGAt                 >  1:207896/1‑132 (MQ=255)                                                                                                                                 atCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGccc     >  2:210188/1‑140 (MQ=255)                                                                                                                                 atCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGccc     >  1:529862/1‑140 (MQ=255)                                                                                                                                  tCATTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGcccc    <  2:499641/140‑1 (MQ=255)                                                                                                                                    aTTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGCCCCcg  <  1:210188/140‑1 (MQ=255)                                                                                                                                    aTTGAGCGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGCCCCcg  <  1:164245/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       TTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAATAGCCCCCG  >  NC_000913/1210356‑1210626

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