Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,197,502	C→T	G11R (GGA→AGA)	ymfD ←	e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,197,502	0	C	T	100.0%	118.5 / NA	40	G11R (GGA→AGA)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/29);  total (11/29)

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTG  >  NC_000913/1197437‑1197554                                                                  |                                                     ctctCAATGATATCTAGTGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                    <  2:929383/68‑1 (MQ=255) cTCGCAATGTTTTCTTGAGTGCATAGTATGAGGTTTTGCTGCATTCTCTGAGCGAATTTTAACTCTGt                                                    <  2:205651/68‑1 (MQ=255) cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                    <  1:96139/68‑1 (MQ=255) cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                    <  1:1046104/68‑1 (MQ=255) cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                    <  2:284298/68‑1 (MQ=255)  tCGCAAAGATATCTAGAGTGCATCGTATGAGGTTTAGCTTCATTCTCTGAGCGAATATTAACTCtgtg                                                   <  1:754204/68‑1 (MQ=255)  acGCAAAGAAATCTAGAGGGAAAAGTATGAGGTTTTGATGAGTTCTCTGCGCGAAAATAAACTCtgtg                                                   >  1:205651/2‑68 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                  <  1:519098/68‑1 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                  <  2:29297/68‑1 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                  <  1:763305/68‑1 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGGTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                  <  1:435287/68‑1 (MQ=255)       aaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCatat                                              <  1:404524/68‑1 (MQ=255)        aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                             <  2:478860/68‑1 (MQ=255)        aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                             <  2:654123/68‑1 (MQ=255)        aGATATCTAGAGGGCATCGTATGAGGTCTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                             <  2:789804/68‑1 (MQ=255)        aGATATCTAGAGGGCATCGTATGAGGCTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                             <  2:674965/68‑1 (MQ=255)        aGACATCTCGAGCGCATCCTATGAGGCCCCGCTGCATTCTCTGAGCGAAAACTAACTCTGTGCataca                                             <  2:1071141/68‑3 (MQ=255)              cTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGtt                                       >  1:646426/1‑68 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                     >  1:937223/1‑68 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                     >  1:306105/1‑68 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                     >  2:315779/1‑68 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                     >  2:625733/1‑68 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCa                                     >  2:430509/1‑68 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCATAg                                  <  1:315779/68‑1 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                  <  2:822534/68‑1 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                  <  2:536870/68‑1 (MQ=255)                    ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGc                                 <  2:937220/68‑1 (MQ=255)                     gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCa                                >  1:917229/1‑68 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAgc                             <  2:991314/68‑2 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                             <  1:523609/68‑1 (MQ=255)                          gTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAc                           >  2:432121/1‑68 (MQ=255)                                ggTTTTGCTGCATTCTCTGATCGAATATTAACTCTTTTCATATTATAGTTCAAAGCAAGTACCATTCa                     <  1:996980/68‑1 (MQ=255)                                ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                     <  1:434906/68‑1 (MQ=255)                                ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                     <  1:625734/68‑1 (MQ=255)                                ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                     <  1:214761/68‑1 (MQ=255)                                ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                     <  1:894518/68‑1 (MQ=255)                                  ttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCTTTCAAt                   <  1:684177/68‑1 (MQ=255)                                      gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt               >  2:962331/1‑68 (MQ=255)                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt             <  2:646425/68‑1 (MQ=255)                                                ctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                    >  2:340640/1‑53 (MQ=255)                                                ctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                    <  1:340640/53‑1 (MQ=255)                                                   aGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTg  >  2:342398/1‑68 (MQ=255)                                                                  |                                                     CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTG  >  NC_000913/1197437‑1197554

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