Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,626	A→G	Y161C (TAT→TGT)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,626	0	A	G	100.0%	184.7 / NA	53	Y161C (TAT→TGT)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/28);  total (25/28)

GGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTG  >  NC_000913/4033551‑4033699                                                                            |                                                                          ggTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTcgcg                                                                      >  2:743956/1‑81 (MQ=255)  gTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCgcgc                                                                     <  2:743885/81‑1 (MQ=255)     tGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAg                                                                   >  1:743961/1‑80 (MQ=255)      ggCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGa                                                                  <  1:743887/80‑1 (MQ=255)       gCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGaa                                                                 <  1:743889/80‑1 (MQ=255)       gCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGaa                                                                 <  2:743887/80‑1 (MQ=255)          ggATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATgc                                                             >  2:743996/1‑81 (MQ=255)             tGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATgcccg                                                          <  2:743889/81‑1 (MQ=255)               gggATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATgcccgg                                                         <  2:743890/80‑1 (MQ=255)               gggATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATgcccgg                                                         <  1:743892/80‑1 (MQ=255)                 gATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATgcccggccc                                                      <  1:743893/81‑1 (MQ=255)                  atcatcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATgcccggcccg                                                     >  2:743952/1‑81 (MQ=255)                     atcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTg                                                  <  1:743897/81‑1 (MQ=255)                     atcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTg                                                  <  1:743894/81‑1 (MQ=255)                     atcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTg                                                  <  1:743896/81‑1 (MQ=255)                     atcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTg                                                  <  1:743895/81‑1 (MQ=255)                        gtgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGaaa                                               <  1:743900/81‑1 (MQ=255)                        gtgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGaaa                                               <  1:743899/81‑1 (MQ=255)                        gtgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGaa                                                <  1:743898/80‑1 (MQ=255)                        gtgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGa                                                 >  1:743951/1‑79 (MQ=255)                          gtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGa                                                 >  2:743969/1‑77 (MQ=255)                           ttGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGa                                             <  2:743900/80‑1 (MQ=255)                              gCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAAc                                         <  1:743902/81‑1 (MQ=255)                                 gTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACaa                                       >  2:743971/1‑80 (MQ=255)                                  ttGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACaaaa                                     >  1:743958/1‑81 (MQ=255)                                    gCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAAt                                    <  2:743902/80‑1 (MQ=255)                                      gATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATgcg                                 >  1:743983/1‑81 (MQ=255)                                      gATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATgcg                                 >  2:743955/1‑81 (MQ=255)                                       aTACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATgcgc                                >  1:743969/1‑81 (MQ=255)                                        tACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGcc                               >  2:743959/1‑81 (MQ=255)                                          cTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCcg                             <  1:743904/81‑1 (MQ=255)                                           tGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCcgc                            >  1:743973/1‑81 (MQ=255)                                              cTAGCCTCGGCGTGGGTGAAATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGaa                         <  2:743904/81‑1 (MQ=255)                                                aTCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAAtt                       >  2:743988/1‑81 (MQ=255)                                                 tCCTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAAtag                      <  2:743905/81‑3 (MQ=255)                                                  ccTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGc                     >  1:743955/1‑81 (MQ=255)                                                   cTCGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGc                     <  2:743906/80‑1 (MQ=255)                                                     cGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCgg                   >  1:743979/1‑80 (MQ=255)                                                     cGGCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCgg                   >  2:743965/1‑80 (MQ=255)                                                      ggCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCgg                   <  1:743908/79‑1 (MQ=255)                                                      ggCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGa                  <  2:743908/80‑1 (MQ=255)                                                      ggCGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGa                  >  2:743961/1‑80 (MQ=255)                                                        cGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCGGGGATTTCCGGaacg               >  1:743948/1‑79 (MQ=255)                                                        cGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAAc               >  1:743971/1‑81 (MQ=255)                                                        cGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAAc               <  1:743911/81‑1 (MQ=255)                                                        cGTGGGTGGGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAAc               <  2:743909/81‑1 (MQ=255)                                                         gtgggtggGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACg              >  2:743963/1‑81 (MQ=255)                                                         gtgggtggGATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACg              >  2:743984/1‑81 (MQ=255)                                                           gggtgggATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGc            <  2:743911/81‑1 (MQ=255)                                                             gtgggATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGa          <  1:743913/81‑1 (MQ=255)                                                               gggATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGaa         >  1:743981/1‑80 (MQ=255)                                                                ggATGCAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAATATGCGCCCGCGAATTGCGGAAACGGCGaaaa       >  2:743958/1‑81 (MQ=255)                                                                     cAGCTCTGTCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCtg  >  2:743975/1‑81 (MQ=255)                                                                            |                                                                          GGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTG  >  NC_000913/4033551‑4033699

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