Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	2,303,559	A→C	S217R (AGT→AGG)	h16_A2116 ←	Putative DNA polymerase bacteriophage‑type,uracil‑DNA glycosylase superfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	2,303,559	0	A	C	100.0%	122.3 / NA	32	S217R (AGT→AGG)	h16_A2116	Putative DNA polymerase bacteriophage‑type,uracil‑DNA glycosylase superfamily
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/16);  total (16/16)

GGTAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACACTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGC  >  AM260479/2303441‑2303682                                                                                                                       |                                                                                                                            ggTAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAgc                                                                                                              >  1:870491/1‑134 (MQ=255)   tAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAgcgc                                                                                                            >  2:750171/1‑134 (MQ=255)             ccacGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTggcgg                                                                                                  >  2:674118/1‑134 (MQ=255)             ccacGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCACAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTggcgg                                                                                                  >  1:555061/1‑134 (MQ=255)                     aTACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCACAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGaa                                                                                          <  2:555061/134‑1 (MQ=255)                         cTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGc                                                                                      <  1:674118/134‑1 (MQ=255)                             aTAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCt                                                                                                                      >  1:366011/1‑98 (MQ=255)                             aTAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCt                                                                                                                      <  2:366011/98‑1 (MQ=255)                             aTAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCg                                                                                   <  2:326049/133‑1 (MQ=255)                                 cTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTg                                                                              >  2:136584/1‑134 (MQ=255)                                       aCGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTgc                                                                             >  2:745401/1‑129 (MQ=255)                                       aCGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTgc                                                                             <  1:745401/129‑1 (MQ=255)                                       aCGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCAc                                                                        >  1:293828/1‑134 (MQ=255)                                       aCGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCAc                                                                        >  1:181204/1‑134 (MQ=255)                                          gTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCAcctc                                                                     >  2:308840/1‑134 (MQ=255)                                               gcgcAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTAACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCgg                                                                >  1:816557/1‑134 (MQ=255)                                               gcgcAGCTTGCCGAGGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCgg                                                                >  2:413341/1‑134 (MQ=255)                                                    gCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTcc                                                           >  2:650335/1‑134 (MQ=255)                                                                 cGTGGTGGTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTAc                                              <  2:794362/134‑1 (MQ=255)                                                                     gtgCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGAt                                                                                                         <  2:357006/71‑1 (MQ=255)                                                                     gtgCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGAt                                                                                                         >  1:357006/1‑71 (MQ=255)                                                                     gtgCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcggg                                          <  2:193803/134‑1 (MQ=255)                                                                     gtgCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcggg                                          <  1:851573/134‑1 (MQ=255)                                                                     gtgCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTAACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcggg                                          <  2:816557/134‑1 (MQ=255)                                                                                    gACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAg                           >  1:328882/1‑134 (MQ=255)                                                                                     aCTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGc                          <  2:952839/134‑1 (MQ=255)                                                                                         cgccgcAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCAc                        <  1:308840/132‑1 (MQ=255)                                                                                         cgccgcAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGt                      <  2:593703/134‑1 (MQ=255)                                                                                                            accacGATCACCCTGGGCTTCCCCAGCGCGATCTGGCGGTGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTgcg   <  2:849693/134‑1 (MQ=255)                                                                                                             ccacGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTgcgc  <  1:650335/134‑1 (MQ=255)                                                                                                                 gATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGaa                                                                                          >  2:249135/1‑42 (MQ=255)                                                                                                                 gATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGaa                                                                                          <  1:249135/42‑1 (MQ=255)                                                                                                                       |                                                                                                                            GGTAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACACTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGC  >  AM260479/2303441‑2303682

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