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	93.3%	34.9 / ‑6.0	15	G11R (GGA→AGA)	ymfD	e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (3/11);  minor base A (0/1);  total (3/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAT  >  NC_000913/1197439‑1197541                                                                |                                        cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                     <  2:162165/68‑1 (MQ=255) cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                     <  2:183810/68‑1 (MQ=255) cGCAAAGATATCTAGAGGGCATCGTATGAGGTTATGCTGCATTCTCTGAGCGAATATTAACTCTGTGc                                     <  1:300800/68‑1 (MQ=255)   cAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCat                                   <  2:383885/68‑1 (MQ=255)      aGATTTCTAGAGTGCATCGTAGGAGGTTTTGCTGCATTCTCTGAGCGAATTTTAACTCTGTGCAtatt                                <  1:561003/68‑1 (MQ=255)           tCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGt                           <  2:731591/68‑1 (MQ=255)                 gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTACAGTTCAAAg                     <  1:575297/68‑1 (MQ=255)                   gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                              >  1:603328/1‑57 (MQ=255)                   gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCAtattat                              <  2:603328/57‑1 (MQ=255)                   gCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCa                   >  2:261008/1‑68 (MQ=255)                             aGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTc         >  2:459445/1‑68 (MQ=255)                              tgTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCTAAGCAAGTACCATTCa        <  2:309473/67‑1 (MQ=255)                              ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCa        <  2:573539/68‑1 (MQ=255)                               gTTTTGATGCATTCTCTGAGCGAATATTAACTTAGTGCATATTATAGTTCAAAGCAAGTACCATTCaa       <  1:532034/68‑1 (MQ=255)                                    gctgcATTCTCTGAGCGAATATTAACTCTGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAt  <  1:507492/68‑1 (MQ=255)                                                                |                                        CGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAAT  >  NC_000913/1197439‑1197541

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: ‑