Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	2,753,670	A→G	W20R (TGG→CGG)	ypfG ←	DUF1176 domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	2,753,670	0	A	G	100.0%	108.1 / NA	28	W20R (TGG→CGG)	ypfG	DUF1176 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/4);  total (24/4)

ACAGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATG  >  NC_017635/2753623‑2753716                                                |                                               aCAGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCcgca                                              >  1:1750882/1‑50 (MQ=255)  cAGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCgcac                                             >  2:1673293/1‑50 (MQ=255)   aGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCgcac                                             >  2:1033080/1‑49 (MQ=255)    ggTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCAcc                                            >  1:1548231/1‑49 (MQ=255)       gACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCa                                           >  1:491979/1‑47 (MQ=255)       gACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCa                                           >  1:601674/1‑47 (MQ=255)           tGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGt                                    >  1:1423768/1‑50 (MQ=255)           tGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGt                                    >  2:351322/1‑50 (MQ=255)           tGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGt                                    >  2:217418/1‑50 (MQ=255)           tGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGt                                    >  1:42776/1‑50 (MQ=255)            gCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTc                                   >  1:553616/1‑50 (MQ=255)              cAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCgg                                 >  1:328367/1‑50 (MQ=255)               aGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCgg                                 >  2:1624905/1‑49 (MQ=255)               aGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGc                                <  1:1879315/50‑1 (MQ=255)                 tCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGcaa                              >  2:793526/1‑50 (MQ=255)                  cGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGcaa                              >  2:767844/1‑49 (MQ=255)                         cGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGa                      <  1:940325/50‑1 (MQ=255)                          gCCCGTTGCGTTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGa                      >  2:396692/1‑49 (MQ=38)                          gCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGa                      >  2:1142214/1‑49 (MQ=255)                             cGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaa                   >  2:487038/1‑49 (MQ=255)                               ttGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaaaa                 >  2:1195751/1‑49 (MQ=255)                               ttGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaaaa                 >  2:1392524/1‑49 (MQ=255)                                    cTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAagag            >  2:1052466/1‑49 (MQ=255)                                     tGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAagag            >  1:868705/1‑48 (MQ=255)                                      ggTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAg         >  1:782820/1‑50 (MQ=255)                                       gTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAg         <  2:220549/49‑1 (MQ=255)                                       gTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAg         >  1:204520/1‑49 (MQ=255)                                       gTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAg         <  2:672223/49‑1 (MQ=255)                                         gccgccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGa       >  1:40007/1‑49 (MQ=255)                                            gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAAt   >  1:776935/1‑50 (MQ=255)                                             cccGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATg  >  1:1073741/1‑50 (MQ=255)                                                |                                               ACAGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATG  >  NC_017635/2753623‑2753716

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