Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	647,162	A→G	V44V (GTA→GTG)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	647,162	0	A	G	100.0%	43.2 / NA	14	V44V (GTA→GTG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/6);  total (8/6)

GTTTTCGGGTCTATTTTTTTATCTGCAATTTTCCCTCCGTTCTGTTAGATGGAAATTAATATTCTTTACAATAACACCTTTCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTAACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAAAGAAGCGGGATATAATACCTTCTTACTTGATTCTAAAGA  >  1/646983‑647408                                                                                                                                                                                    |                                                                                                                                                                                                                                                       gTTTTCGGGTCTATTTTTTTATCTGCAATTTTCCCTCCGTTCTGTTAGATGGAAATTAATATTCTTTACAATAACACCTTTCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCca                                                                                                                                                                                                            >  1:850511/1‑224 (MQ=255)                                                        ttAATATTCTTTACAATAACACCTTTCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAAc                                                                                                                                                         <  1:82687/220‑1 (MQ=255)                                                                 tttACAATAACACCTTTCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCcagc                                                                                                                 <  1:888372/251‑1 (MQ=255)                                                                     cAATAACACCTTTCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTgaga                                                                                                                                                                                    >  1:640751/1‑180 (MQ=255)                                                                               tttCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGaa                                                                                                   >  1:191438/1‑251 (MQ=255)                                                                                                                                                                              atatGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTgag                                                                                                                                                                                     <  1:666753/74‑1 (MQ=255)                                                                                                                                                                              atatGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAaagaag                                    >  1:556150/1‑219 (MQ=255)                                                                                                                                                                              atatGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAAAGAAGCGGGATATAATACCTTCTTACTTGATTCTaaa    >  1:183307/1‑251 (MQ=255)                                                                                                                                                                              atatGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAAAGAAGCGGGATATAATACCTTCTTACTTGATTCTaaa    >  1:298709/1‑251 (MQ=255)                                                                                                                                                                              atatGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAAAGAAGCGGGATATAATACCTTCTTACTTGATTCTaaa    >  1:645261/1‑251 (MQ=255)                                                                                                                                                                                atGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGtat                                                                                                                                                                                                                 <  1:758393/44‑1 (MQ=255)                                                                                                                                                                                atGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCca                                                                                                                                                                                                            <  1:172378/49‑1 (MQ=255)                                                                                                                                                                                atGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAAc                                                                                    <  1:85548/169‑1 (MQ=255)                                                                                                                                                                                atGTGACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAAAGAAGCGGGATATAATACCTTCTTACTTGATTCTAAAGa  >  1:278707/1‑251 (MQ=255)                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GTTTTCGGGTCTATTTTTTTATCTGCAATTTTCCCTCCGTTCTGTTAGATGGAAATTAATATTCTTTACAATAACACCTTTCAAAACTCAATAATTTATCTAAAAATCACGCTATTTTATATAAAAAATTATTTATTTTATAAAAATAAACTAAATTTTAATAATTTTTTATTATATGTAACATTCCTCACAGTAATGCTAAAAATTTTTTATTAGTATAGCCACAACTCGATGATAGATGTTTGAGAAGCGAGTTACCAGCATGGACTTTAAACTCATCATAAACAGGAGACAAATTATGACATACTATCCAGCAGAGCCGTTCCGAATCAAAAGTGTTGAACCGGTTTCCATTTTACCGAAAGAAGAACGCGAAAAAGCCATGAAAGAAGCGGGATATAATACCTTCTTACTTGATTCTAAAGA  >  1/646983‑647408

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