Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	93,029	2 bp→GC	coding (842‑843/1239 nt)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	93,029	0	T	G	100.0%	96.9 / NA	33	K281S (AAA→AGC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (32/1);  total (32/1)
*	1	93,030	0	T	C	100.0%	94.3 / NA	33	K281S (AAA→AGC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/1);  total (32/1)

TTCCGGCTGCAATCGGATAAACGGTGGTTAAAGTGCGGTCATAATTGCCCACGACCGCGCGTTTGAGCTTGATTTTGTTGGATTCAAGCGCACTAATCTCGACCATTTCTTTGTTTTTTCCTGTAATCAAAATGGCACGGCCACCGACGGCAAAATCATAGCCTGCGGTGTTAATAGCAAGCTCTGTGGCGCCCTGTTGCACATTTTGCAACAACTGGGCGCTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAATGTACGGCGCGGCGACAGGCGGCGGGCA  >  1/92914‑93266                                                                                                                    ||                                                                                                                                                                                                                                             ccgccgaacgcagcgggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATc                                                                                                                       >  1:839159/18‑236 (MQ=255)        ccgcagcGGATAAACGGTGTTTAAAGTGCGCTAATATTTCCCCACGCCCGCGCGTTTGAGTTTGTTTTTCCTCGGTTCAGGGGCGTTGTTTTCGACCTTTTTTTTGTTGCTTGCGGTCATAAAAACCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGc                                                                                                 <  1:653339/245‑1 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATAcg                                                                                                        >  1:187236/3‑236 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                                                                         >  1:886681/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                                                                         >  1:80681/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                                                                         >  1:543314/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                                                                         >  1:249910/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)                                     gctCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAgcg                                                                     >  1:752436/3‑250 (MQ=255)                                        cATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGa                                                                 >  1:1960/1‑251 (MQ=255)                                        cATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGa                                                                 >  1:555352/1‑251 (MQ=255)                                        cATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGa                                                                 >  1:788551/1‑251 (MQ=255)                                               gCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAAc                                                                                           >  1:518439/1‑218 (MQ=255)                                               gCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCaa                                                          >  1:13718/1‑251 (MQ=255)                                                                  aGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACttt                                       >  1:746417/1‑251 (MQ=255)                                                                             actcggttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCaaacg                            >  1:818211/23‑249 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCACCATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:762311/18‑251 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:83280/18‑251 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:97017/18‑251 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:762306/18‑251 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:151024/18‑251 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:123917/18‑251 (MQ=255)                                                                                  gttcaagggcgctgatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTAcgg                       >  1:112096/18‑251 (MQ=255)                                                                                               gatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACt                                  >  1:235910/5‑227 (MQ=255)                                                                                               gatttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACaag           >  1:33515/5‑248 (MQ=255)                                                                                                 tttcGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAagcgg        >  1:588933/3‑251 (MQ=255)                                                                                                       ccATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCa  >  1:51729/1‑251 (MQ=255)                                                                                                         aTTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAAcc                                                                                          >  1:658362/1‑161 (MQ=255)                                                                                                                    ||                                                                                                                                                                                                                                             TTCCGGCTGCAATCGGATAAACGGTGGTTAAAGTGCGGTCATAATTGCCCACGACCGCGCGTTTGAGCTTGATTTTGTTGGATTCAAGCGCACTAATCTCGACCATTTCTTTGTTTTTTCCTGTAATCAAAATGGCACGGCCACCGACGGCAAAATCATAGCCTGCGGTGTTAATAGCAAGCTCTGTGGCGCCCTGTTGCACATTTTGCAACAACTGGGCGCTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAATGTACGGCGCGGCGACAGGCGGCGGGCA  >  1/92914‑93266

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