Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	964,882	G→A	R74Q (CGG→CAG)	ECW_RS04675 →	DUF2732 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	964,882	0	G	A	100.0%	65.1 / NA	20	R74Q (CGG→CAG)	ECW_RS04675	DUF2732 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/15);  total (5/15)

CTGCGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCGGGAGATCTACTAATGGCTGATGCAATGGATCTCGTACAGCA  >  NC_017635/964845‑964923                                      |                                          ctgcGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATctact                                <  2:1895662/49‑1 (MQ=255)    cGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAAt                             <  1:1402179/49‑1 (MQ=255)    cGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATg                            <  1:1110212/50‑1 (MQ=255)    cGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATg                            <  1:1298065/50‑1 (MQ=255)    cGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATg                            <  1:1963032/50‑1 (MQ=255)    cGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATg                            <  1:2054239/50‑1 (MQ=255)      tgtgACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATGGc                          <  1:195846/50‑1 (MQ=255)      tgtgACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATGGc                          >  2:1961865/1‑50 (MQ=255)       gtgACTGCTGAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATGGCt                         >  1:864391/1‑50 (MQ=255)             gctgAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGc                    >  2:511543/1‑49 (MQ=255)              ctgAAGCAATCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCaa                  <  1:179365/50‑1 (MQ=255)                   gCAATCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGAt             <  2:649656/50‑1 (MQ=255)                    cAATCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGAtc            >  1:659728/1‑50 (MQ=255)                     aaTCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGAtct           <  1:1843696/50‑1 (MQ=255)                     aaTCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGAtct           <  2:1566550/50‑1 (MQ=255)                      aTCCAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGAtctc          >  1:1229934/1‑50 (MQ=255)                         cAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGATCTCGt        <  1:866152/49‑1 (MQ=255)                         cAGAACGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGATCTCGt        <  1:857030/49‑1 (MQ=255)                            aaCGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGATCTCGTAcag    <  1:214947/50‑1 (MQ=255)                              cGAAGCGCAGGAGATCTACTAATGGCTGATGCAATGGATCTCGTAcagca  <  1:2021251/50‑1 (MQ=255)                                      |                                          CTGCGTGTGACTGCTGAAGCAATCCAGAACGAAGCGCGGGAGATCTACTAATGGCTGATGCAATGGATCTCGTACAGCA  >  NC_017635/964845‑964923

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