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	100.0%	92.0 / NA	25	G97G (GGC→GGG)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/9);  total (16/9)

AAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAA  >  NC_000913/1847670‑1847761                                                  |                                            aaGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTg                                              >  1:850628/1‑49 (MQ=255)    tttGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑ccc                                          >  1:441775/1‑49 (MQ=255)    tttGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑ccc                                          >  1:1878608/1‑49 (MQ=255)       gTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑cccccc                                       >  1:527405/1‑49 (MQ=255)        ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑cccccc                                       >  1:555314/1‑48 (MQ=255)        ttAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑cccccc                                       >  1:1432998/1‑48 (MQ=255)          aaTCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑ccccccca                                     <  1:1020891/48‑3 (MQ=255)                 cAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATcgcg                                 >  1:950237/1‑46 (MQ=38)                 cAGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAg                              <  1:2070546/49‑1 (MQ=38)                  aGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATcgcg                                 >  1:83398/1‑45 (MQ=38)                  aGCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATcgcg                                 >  1:2093312/1‑45 (MQ=38)                   gCCGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGa                             >  1:127356/1‑48 (MQ=38)                    ccGAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGata                           <  1:394533/49‑1 (MQ=38)                      gAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATc                         <  1:1369834/49‑1 (MQ=38)                      gAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATc                         >  1:1228763/1‑49 (MQ=38)                      gAGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATc                         >  1:1825176/1‑49 (MQ=38)                       aGGATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCg                        >  1:7323/1‑49 (MQ=38)                         gATCGCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCt                      <  1:1017397/49‑1 (MQ=38)                             gCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATg                  >  1:2077633/1‑49 (MQ=38)                             gCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATg                  <  1:80064/49‑1 (MQ=38)                             gCAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATg                  >  1:1569998/1‑49 (MQ=38)                              cAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATgg                 >  1:1124785/1‑49 (MQ=38)                              cAATCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATgg                 <  1:1082825/49‑1 (MQ=38)                               aaTCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGc                >  1:1619108/1‑49 (MQ=39)                                aTCGCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCg               >  1:694052/1‑49 (MQ=255)                                   gCCATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGt              >  1:1678806/1‑47 (MQ=255)                                    ccATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAg           <  1:1626962/49‑1 (MQ=39)                                    ccATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAg           >  1:638422/1‑49 (MQ=39)                                    ccATAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAg           <  1:714782/49‑1 (MQ=39)                                       tAATCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcg        <  1:1526180/49‑1 (MQ=39)                                         aTCGGTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgcc      >  1:1274036/1‑49 (MQ=38)                                             gTTTGCCCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCaa  >  1:2109344/1‑49 (MQ=38)                                                  |                                            AAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑CCCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAA  >  NC_000913/1847670‑1847761

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