Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,042,861	A→G	R450R (CGT→CGC)	etk ←	protein‑tyrosine kinase Etk

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,042,861	0	A	G	100.0%	133.3 / NA	37	R450R (CGT→CGC)	etk	protein‑tyrosine kinase Etk
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/19);  total (18/19)

ATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACG  >  NC_000913/1042729‑1042993                                                                                                                                     |                                                                                                                                     aTGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACa                                                                                                                               >  2:689767/1‑140 (MQ=255)      gctgctGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCa                                                                                                                          <  1:969141/140‑1 (MQ=255)          ctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAGACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCAcgcg                                                                                                                      >  1:264310/1‑140 (MQ=255)            gATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCa                                                                                                                          >  2:581452/1‑134 (MQ=255)            gATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCa                                                                                                                          <  1:581452/134‑1 (MQ=255)               tAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGc                                                                                                                 <  2:474842/140‑1 (MQ=255)                    aataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCa                                                                                                                          <  1:238744/126‑1 (MQ=255)                    aataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCa                                                                                                                          >  2:238744/1‑126 (MQ=255)                         atTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCca                                                                                                       <  1:314365/140‑1 (MQ=255)                         atTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCca                                                                                                       <  1:531357/140‑1 (MQ=255)                                         cGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGcc                                                                                       <  2:794867/140‑1 (MQ=255)                                               cGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                 >  2:444590/1‑140 (MQ=255)                                                    aTCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAc                                                                            >  2:1004263/1‑140 (MQ=255)                                                      ccAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                          >  2:154221/1‑140 (MQ=255)                                                       cAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGTTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCaa                                                                         <  1:51871/140‑1 (MQ=255)                                                         gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                       >  2:78594/1‑140 (MQ=255)                                                          ccACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGca                                                                      >  1:632295/1‑140 (MQ=255)                                                              tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcaccac                                                                  >  1:581443/1‑140 (MQ=255)                                                                          aTAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttt                                                      >  1:571867/1‑140 (MQ=255)                                                                            aGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttttt                                                    >  2:112599/1‑140 (MQ=255)                                                                             gTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttttt                                                   <  1:39820/140‑1 (MQ=255)                                                                                gCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttt                                                     <  1:552836/135‑1 (MQ=255)                                                                                gCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttt                                                     >  2:552836/1‑135 (MQ=255)                                                                                gCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                >  1:731259/1‑140 (MQ=255)                                                                                gCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                >  1:840746/1‑140 (MQ=255)                                                                                gCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                >  1:337225/1‑140 (MQ=255)                                                                                 cATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGt                                               <  2:403740/140‑1 (MQ=255)                                                                                        gCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTg                                        >  2:551248/1‑140 (MQ=255)                                                                                          tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGc                                      >  1:779318/1‑140 (MQ=255)                                                                                            aTGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTg                                    <  1:689767/140‑1 (MQ=255)                                                                                             tGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGc                                   <  1:444590/140‑1 (MQ=255)                                                                                                 gTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCt                               <  1:1004263/140‑1 (MQ=255)                                                                                                 gTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCt                               <  2:779318/140‑1 (MQ=255)                                                                                                    ctctTCCAGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAg                            <  2:581443/140‑1 (MQ=255)                                                                                                           aGTTGTTCCGGGGCTTCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGc                     <  2:262347/140‑1 (MQ=255)                                                                                                                          tCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCg      <  1:629409/140‑1 (MQ=255)                                                                                                                          tCTACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCg      <  2:632295/140‑1 (MQ=255)                                                                                                                            tACACCACGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCa    <  2:786451/140‑1 (MQ=255)                                                                                                                              caccacGGCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGATCGATAATCCGCACg  <  2:150109/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     ATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACG  >  NC_000913/1042729‑1042993

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