Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,436,697	G→A	E60E (GAG→GAA)	trkA →	NAD‑binding component of Trk potassium transporters

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,436,697	0	G	A	100.0%	54.2 / NA	20	E60E (GAG→GAA)	trkA	NAD‑binding component of Trk potassium transporters
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/9);  total (11/9)

GGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTG  >  NC_000913/3436458‑3436883                                                                                                                                                                                                                                                |                                                                                                                                                                                           ngCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCcg                                                                                                                                                                                  >  1:150849/2‑250 (MQ=255)       ggCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGCTAATTATCTTTCTGGGTGCCGGCCAGGTTGGCGGCACTCTGGCGGAATACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGc                                                                                                                                                                           <  1:480818/251‑1 (MQ=255)                               aaaaaGTGATGAGAGAACGGGTCGCGACTGATTAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATTCCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGGAGGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAg                                                                                                                                                   <  1:351936/251‑1 (MQ=255)                                                      gcgACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAACACCTGGTTGACGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACATCACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGAGTATTGCGGGAAGCAGGTGCCGAAGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGttt                                                                                                                            >  1:350628/1‑250 (MQ=255)                                                                                                  ncacaCTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGt                                                                                                                                                  >  1:311728/2‑185 (MQ=255)                                                                                                     cacTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCt                                                                             >  1:115169/1‑251 (MQ=255)                                                                                                                          gCGAGNNCANNNATNNTACTGTTGTCGATACCAACGGTGNGCGTCTGCGGACCTTACAGGATCAATTTGNCNTNCCGGNNGTGNAGGGGCATGGNTCTCATCCACGCGTCTTGCGGGAAGCAGGTGCCGACGACGCCGATATGNTNGTTGCTGTAACCAGTTCAGATGAAACCAATANGGTtgcctgcc                                                                                                                      >  2:374193/1‑189 (MQ=255)                                                                                                                          gCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTGGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTAGTGCGGGAAGCAGGTTCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcctgcc                                                                                                                      <  1:374193/189‑1 (MQ=255)                                                                                                                                 nnnnnnnnnnaCNNNNNTCGNTACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGt                                                                                                                                                  <  2:311728/144‑1 (MQ=255)                                                                                                                                                 ntCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAata                                                                                                                                  >  1:302453/2‑154 (MQ=255)                                                                                                                                                 nnnnnnnnnnaCNNNNNGCGNCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAata                                                                                                                                  <  2:302453/144‑1 (MQ=255)                                                                                                                                                 gtCGATACCAACGGTGAGCGTATGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAATCAATATGGTTGCCTGCCAg                                                                                                                    >  1:488043/1‑168 (MQ=255)                                                                                                                                                       aCCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACttt                                                                                                    >  1:219736/1‑178 (MQ=255)                                                                                                                                                                                 gnaGGNTNNANNTGCNCGGCGGGTCGTGCAGGGGCATNNCTCTCATCCNCGCGTATGGCGGGAAGCAGGTGCCGACGACGCCGATAAGCTGGTTTCTGTAACNAGTTCAGATNNAANCAATANGNNNGNNNGCNAGGTAGCCTACTCACTTTTNAANNCCCCTANTCGCATCGCTCGTATCCGCTCACCAg                                                             <  2:115169/189‑1 (MQ=255)                                                                                                                                                                                                       gTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTan                                                                                                                <  1:309922/118‑2 (MQ=255)                                                                                                                                                                                                       gTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAg                                                                                                                >  2:309922/1‑118 (MQ=255)                                                                                                                                                                                                                   cATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGACGACGACGCCGATATGCTGGTTGCTGTCACAAGTTCAGATGAAAcc                                                                                                                                      >  1:341579/1‑84 (MQ=255)                                                                                                                                                                                                                    aTGGCTCTCATCCACGCGTATTGCGGNAAGCAGNNGcn                                                                                                                                                                                   >  2:117434/1‑37 (MQ=255)                                                                                                                                                                                                                    aTGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGcc                                                                                                                                                                                   <  1:117434/38‑1 (MQ=255)                                                                                                                                                                                                                                          gCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTg  <  1:68085/193‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                           GGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTG  >  NC_000913/3436458‑3436883

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