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%	79.5 / NA	26	G11R (GGA→AGA)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/15);  total (11/15)

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTT  >  NC_000913/1197437‑1197565                                                                  |                                                                cTCGCAAAGTTATCTATTTGGCATCGTATGTGGTTTTGCTTCATTCTCTTAGCGAATATTAACTCTGt                                                               <  2:100907/68‑1 (MQ=255)  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               >  1:141473/1‑67 (MQ=255)  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  2:141474/67‑1 (MQ=255)   cGCAAAGATATCTTGATTTCATCTTATGAGGCTTTGCTGCATTCTCTGAGCGAATTTTAACTCTGTGc                                                             <  2:812630/68‑1 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  1:180062/68‑1 (MQ=255)     cAAAGATATCTAGAGTGCATCGTATTAGGTTTTGCTGCATTCTCTGAGCGAATATTATCTCTGTGCat                                                           <  2:559693/68‑1 (MQ=255)     cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                           <  2:431722/68‑1 (MQ=255)        aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATAAACTCTGTGCAtatt                                                        <  2:670124/68‑1 (MQ=255)          atatCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                                                      >  1:272180/1‑68 (MQ=255)             tCTAGAGTGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                   <  2:704598/68‑1 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                   <  1:430242/65‑1 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                                                   >  2:430243/1‑65 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             <  2:242731/68‑1 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             >  2:612317/1‑68 (MQ=255)                       aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  1:616232/48‑1 (MQ=255)                       aTCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             >  2:616233/1‑48 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                        <  2:170025/68‑1 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                        <  1:612316/68‑1 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                        <  1:552168/68‑1 (MQ=255)                           tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAcc                                     >  2:421779/1‑68 (MQ=255)                                      gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                          <  2:684474/68‑1 (MQ=255)                                                ctGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTa                >  1:537618/1‑68 (MQ=255)                                                     cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaa           >  2:404600/1‑68 (MQ=255)                                                     cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaa           >  1:250188/1‑68 (MQ=255)                                                             aaCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCtt   >  1:233946/1‑68 (MQ=255)                                                              aCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttt  >  2:623498/1‑68 (MQ=255)                                                                  |                                                                CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTT  >  NC_000913/1197437‑1197565

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