Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	92,951	G→C	D307E (GAC→GAG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	92,951	0	G	C	100.0%	30.8 / NA	13	D307E (GAC→GAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (12/1);  total (12/1)

TTCCGGCTGCAATCGGATAAACGGTGGTTAAAGTGCGGTCATAATTGCCCACGACCGCGCGTTTGAGCTTGATTTTGTTGGATTCAAGCGCACTAATCTCGACCATTTCTTTGTTTTTTCCTGTAATCAAAATGGCACGGCCACCGACGGCAAAATCATAGCCTGCGGTGTTAATAGCAAGCTCTGTGGCGCCCTGTTGCACATTTTGCAACAACTGGGCGCTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCG  >  1/92914‑93186                                      |                                                                                                                                                                                                                                            ccgccgaacgcagcgggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATc                                       >  1:839159/18‑236 (MQ=255)        ccgcagcGGATAAACGGTGTTTAAAGTGCGCTAATATTTCCCCACGCCCGCGCGTTTGAGTTTGTTTTTCCTCGGTTCAGGGGCGTTGTTTTCGACCTTTTTTTTGTTGCTTGCGGTCATAAAAACCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGc                 <  1:653339/245‑1 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATAcg                        >  1:187236/3‑236 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa         >  1:249910/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa         >  1:543314/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa         >  1:80681/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa         >  1:886681/3‑251 (MQ=255)                   aaaCAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATag      >  1:862101/1‑250 (MQ=255)                        gtggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAgcg  >  1:14596/1‑250 (MQ=255)                         tggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTaaa                                           >  1:462896/1‑208 (MQ=255)                          ggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAAc           >  1:28546/1‑239 (MQ=255)                          ggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAgc   >  1:491700/1‑247 (MQ=255)                                 gTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGAcc                             >  1:104724/1‑214 (MQ=255)                                      |                                                                                                                                                                                                                                            TTCCGGCTGCAATCGGATAAACGGTGGTTAAAGTGCGGTCATAATTGCCCACGACCGCGCGTTTGAGCTTGATTTTGTTGGATTCAAGCGCACTAATCTCGACCATTTCTTTGTTTTTTCCTGTAATCAAAATGGCACGGCCACCGACGGCAAAATCATAGCCTGCGGTGTTAATAGCAAGCTCTGTGGCGCCCTGTTGCACATTTTGCAACAACTGGGCGCTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCG  >  1/92914‑93186

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