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

CTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTT  >  NC_000913/1197437‑1197565                                                                  |                                                                cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  1:131583/68‑1 (MQ=255) cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGt                                                               <  2:143150/68‑1 (MQ=255)  tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCtgtg                                                              >  2:405014/1‑68 (MQ=255)   cGCATAGATATCTAGAGGGCTTCTTATGATTTTTTGCTTCATTCTCTGAGCTAATATTAACTCTGttc                                                             <  1:661447/68‑3 (MQ=255)   cGCAAAGATATCTAGAGTGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  2:526431/68‑1 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  2:332314/68‑1 (MQ=255)   cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                                             <  2:407410/68‑1 (MQ=255)     cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                                           <  2:471048/68‑1 (MQ=255)        aGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtttt                                                        >  2:673892/1‑68 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                        >  1:401636/1‑60 (MQ=255)                agagGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                        <  2:401636/60‑1 (MQ=255)                  agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                        >  1:893340/1‑58 (MQ=255)                  agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtatt                                                        <  2:893337/58‑1 (MQ=255)                  agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCaaa                                              >  1:925355/1‑68 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             >  1:916195/1‑68 (MQ=255)                   gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             >  2:792310/1‑68 (MQ=255)                   gggCATCGTATGAGGTTTAGCTGCAATCTATGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAg                                             >  1:254406/1‑68 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                        <  2:549396/68‑1 (MQ=255)                        tCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGt                                        <  2:925352/68‑1 (MQ=255)                           tATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAcc                                     <  2:220411/68‑1 (MQ=255)                            aTGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCa                                    <  2:450077/68‑1 (MQ=255)                                gtttttGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                                <  1:599490/66‑1 (MQ=255)                                ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa                                <  1:288718/68‑1 (MQ=255)                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                               >  1:718005/1‑68 (MQ=255)                                 gTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCaa                               <  1:4094/68‑1 (MQ=255)                                      gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                          <  2:916192/68‑1 (MQ=255)                                      gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt                          <  1:108675/68‑1 (MQ=255)                                        tgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGt                        <  1:39798/68‑1 (MQ=255)                                             tctcTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAcc                                     <  1:12631/50‑1 (MQ=255)                                             tctcTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTAcc                                     >  2:12631/1‑50 (MQ=255)                                                    gCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGa            >  1:50764/1‑68 (MQ=255)                                                     cGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGaa           >  1:942409/1‑68 (MQ=255)                                                              aCTCTGTGCATATTATAGTTCAAAGCAAGTCCCATTCAATTAATGTTATTTTTAGTGAAAAATTCttt  >  1:762457/1‑68 (MQ=255)                                                              aCTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCttt  >  1:89574/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑