Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,309,016	C→T	*418* (TAG→TAA)	kch ←	voltage‑gated potassium channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,309,016	0	C	T	100.0%	133.9 / NA	42	*418* (TAG→TAA)	kch	voltage‑gated potassium channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/15);  total (27/15)

GCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTCTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTAAAAC  >  NC_000913/1308909‑1309144                                                                                                            |                                                                                                                                 gCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCtt                                                                                                   >  2:641474/1‑139 (MQ=255)       tAAACAAACAGCGTAACTCGTTATTGTTTGTGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGttt                                                                                             >  1:305938/1‑139 (MQ=255)          acaaacaGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCtt                                                                                                               <  1:628283/118‑1 (MQ=255)          acaaacaGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCtt                                                                                                               >  2:628283/1‑118 (MQ=255)                 gCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATc                                                                                   >  2:492901/1‑139 (MQ=255)                 gCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATc                                                                                   >  1:736750/1‑139 (MQ=255)                 gCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATc                                                                                   >  2:743110/1‑139 (MQ=255)                   gTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                     <  1:166627/103‑1 (MQ=255)                   gTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                     >  2:166627/1‑103 (MQ=255)                       cTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATc                                                                             <  1:107521/139‑1 (MQ=255)                          gTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATc                                                                             <  2:742285/136‑1 (MQ=255)                          gTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATc                                                                             >  1:742285/1‑136 (MQ=255)                              ttGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGaaa                                                                      <  1:743110/139‑1 (MQ=255)                                                  cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTAtttt                                                                                                                            <  2:551101/65‑1 (MQ=255)                                                  cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTAtttt                                                                                                                            >  1:551101/1‑65 (MQ=255)                                                  cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGAc                                                            >  1:126914/1‑129 (MQ=255)                                                  cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGAc                                                            <  2:126914/129‑1 (MQ=255)                                                  cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTg                                                         <  1:330827/132‑1 (MQ=255)                                                  cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTg                                                         >  2:330827/1‑132 (MQ=255)                                                    aaaTGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaaca                                                >  2:676052/1‑139 (MQ=255)                                                    aaaTGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaaca                                                >  1:715293/1‑139 (MQ=255)                                                         tACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACAtt                                           >  1:138116/1‑139 (MQ=255)                                                         tACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACAtt                                           >  2:554098/1‑139 (MQ=255)                                                           cTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGa                                         >  1:223530/1‑139 (MQ=255)                                                             gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                     >  1:162475/1‑61 (MQ=255)                                                             gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                     <  2:162475/61‑1 (MQ=255)                                                             gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaaa                                       >  1:210531/1‑139 (MQ=255)                                                               gCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACg                                     >  2:308758/1‑139 (MQ=255)                                                                   cATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCa                                 <  2:736750/139‑1 (MQ=255)                                                                      tGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCatat                              <  2:715293/139‑1 (MQ=255)                                                                             cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCattatt                       >  2:479895/1‑139 (MQ=255)                                                                             cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCattatt                       >  2:677465/1‑139 (MQ=255)                                                                             cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCattatt                       >  1:354355/1‑139 (MQ=255)                                                                                       ttACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCa             >  1:210099/1‑139 (MQ=255)                                                                                       ttACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCa             <  2:223530/139‑1 (MQ=255)                                                                                       ttACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCa             <  2:210531/139‑1 (MQ=255)                                                                                       ttACTCTTCCGTCTTTAAACTTTATTATTGCGCCGATCCTGTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGCAAATACCATGACCTGAGTTCAACAACATTGAATCGAGCATATCATTATTAATCTCTTCa             <  2:378042/139‑1 (MQ=255)                                                                                          ctctTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCa          >  1:560061/1‑139 (MQ=255)                                                                                               ccGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTaa     >  2:659542/1‑139 (MQ=255)                                                                                                cGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCAt                                                               <  1:289595/80‑1 (MQ=255)                                                                                                cGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCAt                                                               >  2:289595/1‑80 (MQ=255)                                                                                                  tCTTTAAACGTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTAAAAc  >  1:516840/1‑139 (MQ=255)                                                                                                            |                                                                                                                                 GCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTCTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTAAAAC  >  NC_000913/1308909‑1309144

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