Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,847,724	(C)6→7	coding (285/1044 nt)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,847,718	1	.	C	90.9%	32.2 / ‑2.5	11	coding (291/1044 nt)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (7/3);  total (8/3)
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

AGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAT  >  NC_000913/1847671‑1847762                                                  |                                              aGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑cc                                            >  1:534214/1‑50 (MQ=255)  gTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑ccc                                           <  1:1794657/50‑1 (MQ=255)   tttGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑cccc                                          >  1:585504/1‑50 (MQ=255)     tGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑cccccc                                        >  1:864199/1‑50 (MQ=255)     tGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑cccccc                                        >  1:513766/1‑50 (MQ=255)      gTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑cccccc                                        >  1:1133607/1‑49 (MQ=255)      gTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑cccccc                                        >  1:2267876/1‑49 (MQ=255)       ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑ccccccc                                       <  1:1321719/49‑2 (MQ=255)       ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑ccccccc                                       <  1:790500/49‑2 (MQ=255)                  gCCGAGGATCGCAATCGCCATAATCGG‑TTTGC‑CCCCCCATCGCGAAGata                            >  1:270511/1‑50 (MQ=38)                      aGGATCGCAATCGCCATAATCGG‑TTTGCCCCCCCCATCGCGAAGATATCg                         >  1:462472/1‑50 (MQ=37)                       ggATCGCAATCGCCATAATCGG‑TTTGC‑CCCCCCATCGCGAAGATATCGCt                       >  1:2239338/1‑50 (MQ=38)                        gATCGCAATCGCCATAATCGG‑TTTGC‑CCCCCCATCGCGAAGATATCGCt                       >  1:896725/1‑49 (MQ=38)                             cAATCGCCATAATCGG‑TTTGC‑CCCCCCATCGCGAAGATATCGCTGATGGc                 >  1:905804/1‑50 (MQ=39)                                  gCCATAATCGGTTTTC‑‑CCCCCCATCGCGAAGATATCGCTGATGATGTTAg            >  1:1256339/1‑50 (MQ=255)                                    cATAATCGG‑TTTGC‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTc          >  1:2056707/1‑50 (MQ=255)                                    cATAATCGG‑TTTGC‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTc          <  1:1186589/50‑1 (MQ=255)                                     aTAATCGG‑TTTGC‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcg         >  1:320934/1‑50 (MQ=255)                                           gg‑TTTGC‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCaa   <  1:210596/50‑1 (MQ=39)                                            g‑TTTGC‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAt  <  1:2043779/50‑1 (MQ=39)                                                  |                                              AGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGG‑TTTG‑‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAAT  >  NC_000913/1847671‑1847762

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