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%	214.8 / NA	66	*418* (TAG→TAA)	kch	voltage‑gated potassium channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/29);  total (37/29)

TCCACAGCATCAGGAAGTGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTCTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTAAAACTCGCG  >  NC_000913/1308891‑1309149                                                                                                                              |                                                                                                                                      tCCACAGCATCAGGAAGTGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                          >  2:753512/1‑139 (MQ=255)          tCAGGAAGTGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTAc                                                                                                                 >  1:11610/1‑139 (MQ=255)             ggAAGTGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTAAAACTTTATTTTTGCGCCGATTCTTTACTGt                                                                                                              >  1:774513/1‑139 (MQ=255)                  tGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCt                                                                                                         <  1:513817/139‑1 (MQ=255)                  tGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCt                                                                                                         >  2:620074/1‑139 (MQ=255)                     aTTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTc                                                                                                      <  1:620074/139‑1 (MQ=255)                            acaaacaGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTa                                                                                               >  1:35139/1‑139 (MQ=255)                                  aGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCAt                                                                                         <  1:753512/139‑1 (MQ=255)                                   gCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATc                                                                                        >  1:412032/1‑139 (MQ=255)                                                atGTTTTCCATCCACAATATCCAAATTTACTGGGCACCATTTCCGGTCCCAATAAGGTTTCCTTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGaaa                                                                           <  2:774513/138‑1 (MQ=255)                                                          tCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGAc                                                                 <  1:99114/139‑1 (MQ=255)                                                             aCAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTg                                                              >  2:345920/1‑139 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGt                                                                                             >  1:585367/1‑101 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGt                                                                                             <  2:585367/101‑1 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       >  2:523537/1‑139 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       >  1:677403/1‑139 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       >  2:718215/1‑139 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       >  1:9640/1‑139 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       <  2:11610/139‑1 (MQ=255)                                                                    cTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       >  2:78485/1‑139 (MQ=255)                                                                     tAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaac                                                      <  1:429553/139‑1 (MQ=255)                                                                        aTGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgcgc                                                                                                                            <  2:194318/66‑1 (MQ=255)                                                                        aTGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgcgc                                                                                                                            >  1:194318/1‑66 (MQ=255)                                                                         tGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAacaaca                                                  <  2:244934/139‑1 (MQ=255)                                                                           tACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACAtt                                                >  1:301649/1‑139 (MQ=255)                                                                            ccTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgcgc                                                                                                                            >  1:502648/2‑62 (MQ=255)                                                                            aCTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTCAAACTTTATTTTTgcgc                                                                                                                            <  2:502648/62‑1 (MQ=255)                                                                             cTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGAc                                                                 <  2:776083/120‑1 (MQ=255)                                                                             cTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGAc                                                                 >  1:776083/1‑120 (MQ=255)                                                                             cTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGa                                              >  2:692717/1‑139 (MQ=255)                                                                             cTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGa                                              >  1:178750/1‑139 (MQ=255)                                                                              tGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaa                                             <  2:35139/139‑1 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGt                                                                                      <  1:540371/97‑1 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGt                                                                                      >  2:540371/1‑97 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCa                                                                     >  1:248188/1‑114 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCa                                                                     <  2:248188/114‑1 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaaa                                            >  1:576351/1‑139 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaaa                                            >  2:117311/1‑139 (MQ=255)                                                                               gggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaaa                                            <  1:523537/139‑1 (MQ=255)                                                                                ggCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAAc                                           <  2:412032/139‑1 (MQ=255)                                                                                 gCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACg                                          >  1:83559/1‑139 (MQ=255)                                                                                     cATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAacaaca                                                  <  2:753657/127‑1 (MQ=255)                                                                                     cATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAacaaca                                                  >  1:753657/1‑127 (MQ=255)                                                                                     cATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCa                                      >  1:766951/1‑139 (MQ=255)                                                                                     cATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCa                                      <  1:345920/139‑1 (MQ=255)                                                                                       ttGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                          >  1:496159/1‑53 (MQ=255)                                                                                       ttGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCg                                                                                                                          <  2:496159/53‑1 (MQ=255)                                                                                               gccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgcgc                                                                                                                            >  1:502406/2‑43 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgcgc                                                                                                                            <  2:502406/43‑1 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgc                                                                                                                              <  2:103040/41‑1 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTgc                                                                                                                              >  1:103040/1‑41 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCg                                                                                 <  2:105668/86‑1 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCg                                                                                 >  1:105668/1‑86 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaaa                                            <  2:412049/123‑1 (MQ=255)                                                                                               cccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGaaa                                            >  1:412049/1‑123 (MQ=255)                                                                                                ccAATAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCAttatta                           <  1:117311/139‑1 (MQ=255)                                                                                                  aaTAAGTTTACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAAt                         <  1:565215/139‑1 (MQ=255)                                                                                                          tACTCTTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCAc                 >  2:407473/1‑139 (MQ=255)                                                                                                              ctTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       >  1:73276/1‑97 (MQ=255)                                                                                                              ctTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTcaa                                                       <  2:73276/97‑1 (MQ=255)                                                                                                              ctTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCa                  <  1:210258/134‑1 (MQ=255)                                                                                                              ctTCCGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCa                  >  2:210258/1‑134 (MQ=255)                                                                                                                  cGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCa                                      >  1:406459/1‑110 (MQ=255)                                                                                                                  cGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCa                                      <  2:406459/110‑1 (MQ=255)                                                                                                                  cGTCTTTAAACTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTaaa         >  2:459164/1‑139 (MQ=255)                                                                                                                         aaaCTTTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTAAAACTcgcg  >  1:280443/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                      TCCACAGCATCAGGAAGTGCATTCTAAACAAACAGCGTAACTCGTTATTGTTTGCGATCTACAATATCTAAATGTACTGGGCACCATTGCGGTGCCCAATAAGTTTACTCTTCCGTCTTTAAACTCTATTTTTGCGCCGATTCTTTACTGTCAGCTTTCGTTTCTAGTTCATCGTTATCGCTGAAAATACCATGACCTGAGTTCAACAACATTGAAACGAGCATATCATTATTAATCTCTTCACCATTTAAAACTCGCG  >  NC_000913/1308891‑1309149

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