Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,080,981	G→A	Y395Y (TAC→TAT)	unknown ←	Dca

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,080,981	0	G	A	100.0%	34.3 / NA	14	Y395Y (TAC→TAT)	unknown	Dca
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/5);  total (9/5)

CACCAAAATCTAAATCTATTTCTTCTTTGCTTTCTTTATTTTTGTGCGATAATCCATGATCTGAAAAATAAATTAAAGAATAACTTTCATTCTGTGCTTTTAATAAATAGACGATCTCTTCAATTAATTTATCCGTTTTTAAAATACTATTGACATAACAGGATAAACTTTCATTTATAAATTCGTATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTGGCAAGAGTGTCATAT  >  1/1080797‑1081206                                                                                                                                                                                         |                                                                                                                                                                                                                                  cACCAAAATCTAAATCTATTCCTTCTTTGCTTTCTTTATTTTTGTGCGATAATCCATGATCTGAAAAATAAATTAAAGAATAACTTTCATTCTGTGCTTTTAATAAATAGACGAGCTCTTCAATTAATTTATCCGTTTTTAAAATACTATTGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATc                                                                                                                                                                             <  1:19167/239‑1 (MQ=255)                                                 aTAATCCATGATCTGAAAAATAAATTAAAGAATAACTTTCATTCTGTGCTTTTAATAAATAGACGATCTCTTCAATTAATTTATCCGTTTTTAAAATACTATTGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTAttt                                                                                                                                               <  1:900223/221‑1 (MQ=255)                                                                                                gCTTTTAATAAATAGACGATCTCTTCAATTAATTTATCCGTTTTTAAAATACTATTGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACtt                                                                                                                                      <  1:554634/183‑1 (MQ=255)                                                                                                                                                   cTATTGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTgg                >  1:42767/1‑250 (MQ=255)                                                                                                                                                   cTATTGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATCACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTgg                >  1:227393/1‑250 (MQ=255)                                                                                                                                                     aTTGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATc                                       >  1:757416/1‑225 (MQ=255)                                                                                                                                                       tGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGt                                                                                                     >  1:134386/1‑161 (MQ=255)                                                                                                                                                       tGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGt                                                                           >  1:160083/1‑187 (MQ=255)                                                                                                                                                       tGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATc                                       >  1:520878/1‑223 (MQ=255)                                                                                                                                                       tGACATAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTGGCAaga           >  1:58150/1‑251 (MQ=255)                                                                                                                                                            tAACAGGATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGAt                                                                                                                                                                                   <  1:226543/78‑1 (MQ=255)                                                                                                                                                                 ggATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGa                        >  1:517130/1‑228 (MQ=255)                                                                                                                                                                 ggATAAACTTTCATTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTGGCAAGAGTGTCatat  >  1:317422/1‑250 (MQ=255)                                                                                                                                                                             aTTTATAAATTCATATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTGGCAAGAgt         <  1:384265/231‑1 (MQ=255)                                                                                                                                                                                         |                                                                                                                                                                                                                                  CACCAAAATCTAAATCTATTTCTTCTTTGCTTTCTTTATTTTTGTGCGATAATCCATGATCTGAAAAATAAATTAAAGAATAACTTTCATTCTGTGCTTTTAATAAATAGACGATCTCTTCAATTAATTTATCCGTTTTTAAAATACTATTGACATAACAGGATAAACTTTCATTTATAAATTCGTATTTCTTTTCATCTTTTGTGATTCTCTGGCAAAAATTACGATGAGATCCCATCAAATGAATGACAAATAGTCGTGTATTATTTGCATAACTTTTCTCTTTAAATTTGACTTTAAGTTCATCCAGTAATTTGAAATCATCGGTATTATTGGTAATAAATCCACCTTTTTTCGTGAAGTAATGAAAATCTGCACTTGCACCTGATCTGCTGGCAAGAGTGTCATAT  >  1/1080797‑1081206

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: ‑