Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	729,200	G→A	intergenic (‑248/‑145)	unknown ← / → unknown	hypothetical protein/Mobile element protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	729,200	0	G	A	100.0%	36.3 / NA	21	intergenic (‑248/‑145)	unknown/unknown	hypothetical protein/Mobile element protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/14);  total (7/14)

CTAATGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAG  >  1/728989‑729290                                                                                                                                                                                                                    |                                                                                             cTAATGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAag                                                       <  1:490904/251‑1 (MQ=17)    aTGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAgg                                                    >  1:237593/1‑251 (MQ=17)           ttttttCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAAt                                             <  1:624309/251‑1 (MQ=17)             ttttCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATgg                                           <  1:9734/251‑1 (MQ=17)             ttttCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATgg                                           <  1:798772/251‑1 (MQ=17)                           tCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                             >  1:155840/1‑251 (MQ=14)                           tCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                             >  1:853646/1‑251 (MQ=14)                           tCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                             >  1:616262/1‑251 (MQ=14)                            ccTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAagag                                                     <  1:183913/226‑1 (MQ=14)                            ccTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTta                            >  1:323201/1‑251 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:607018/251‑1 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:691095/251‑1 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:898492/251‑1 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:110887/251‑1 (MQ=14)                                    tATGCAAATTTTACCGGATTTTTTGATCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc                    <  1:548875/251‑1 (MQ=14)                                     aTGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                   <  1:481439/251‑1 (MQ=14)                                                     aTTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt    <  1:35198/250‑1 (MQ=11)                                                     aTTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt    <  1:170623/250‑1 (MQ=11)                                                      ttttttGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt    >  1:25008/1‑249 (MQ=11)                                                      ttttttGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg  <  1:729281/251‑1 (MQ=11)                                                      ttttttGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg  >  1:840211/1‑251 (MQ=11)                                                                                                                                                                                                                    |                                                                                             CTAATGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAG  >  1/728989‑729290

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