Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	245,288	T→G	L247V (TTG→GTG)	betA2 →	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	245,288	0	T	G	100.0%	75.7 / NA	21	L247V (TTG→GTG)	betA2	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/7);  total (14/7)

GCGTGCCTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGTTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCA  >  AM260479/245158‑245403                                                                                                                                   |                                                                                                                    gcgTGCCTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTgg                                                                                                                  >  1:573901/1‑134 (MQ=255)                       aacaaGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAAct                                                                                            >  1:335683/1‑134 (MQ=255)                        acaaGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAActg                                                                                           >  2:432433/1‑134 (MQ=255)                                       ccTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGcgccgc                                                                            <  2:410233/134‑1 (MQ=255)                                                    gcgcgcgCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCgtgt                                                               >  2:57199/2‑134 (MQ=255)                                                    gcgcgcgCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCgtgt                                                               >  2:518045/2‑134 (MQ=255)                                                    ccGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCgtgt                                                               >  2:174442/1‑134 (MQ=255)                                                    ccGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCgtgt                                                               >  2:400788/1‑134 (MQ=255)                                                       cgcgcACCTGCGGGTGCTGACAGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTgcg                                                            <  2:1038410/134‑1 (MQ=255)                                                                  gggTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTc                                                 <  1:518045/134‑1 (MQ=255)                                                                  gggTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTc                                                 <  1:57199/134‑1 (MQ=255)                                                                       cTGACCGATACGCAGGCGCTGCTCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTc                                            <  2:573901/134‑1 (MQ=255)                                                                          aCCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTcgcc                                         >  2:1061347/1‑134 (MQ=255)                                                                                   cAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAgctgct                                >  2:691397/1‑134 (MQ=255)                                                                                        gcTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATgg                           >  2:423811/1‑134 (MQ=255)                                                                                                tAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGc                   >  1:999048/1‑134 (MQ=255)                                                                                                     ttCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTcgcc                                         >  2:151895/1‑107 (MQ=255)                                                                                                     ttCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTcgcc                                         <  1:151895/107‑1 (MQ=255)                                                                                                         aaGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGcccc          <  1:174442/134‑1 (MQ=255)                                                                                                                 cATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCa  >  2:543447/1‑134 (MQ=255)                                                                                                                 cATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCa  >  2:569975/1‑134 (MQ=255)                                                                                                                                   |                                                                                                                    GCGTGCCTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGTTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCA  >  AM260479/245158‑245403

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