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	96.5%	94.9 / NA	29	E60E (GAG→GAA)	trkA	NAD‑binding component of Trk potassium transporters
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/16);  minor base C (1/0);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGC  >  NC_000913/3436451‑3436942                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAg                                                                                                                                                                                                                                                   <  2:177585/251‑1 (MQ=255)                cttcttTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCg                                                                                                                                                                                                                                    >  1:135319/1‑251 (MQ=255)                     tttACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAGACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATg                                                                                                                                                                                                                               >  2:13950/1‑251 (MQ=255)                               gctaaTCCAAAAGTGATGCGCTAACGGGTCGCGACTGATGAAAATTCTCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGACAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGCCCTTACAGGATAAATTTGACCTGAGGGTCGTGCAGGGGCATGGCTCTAATCCACGCGTATTGCGGGACGCAGGTGCCGACGACGCCGCTATGCTGGTTGCtt                                                                                                                                                                                                                     >  2:21825/1‑250 (MQ=255)                                             gatgaGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAGACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGAt                                                                                                                                                                                                       <  1:13950/251‑1 (MQ=255)                                                gagaTAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGaa                                                                                                                                                                                                    <  1:196185/251‑1 (MQ=255)                                                                                  tCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCa                                                                                                                                                                   <  1:231865/250‑1 (MQ=255)                                                                                                        ggcACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcgcat                                                                                                                                            <  1:123523/251‑1 (MQ=255)                                                                                                        ggcACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATGACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcgcat                                                                                                                                            <  1:271929/251‑1 (MQ=255)                                                                                                          cacaCTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACAAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGCCGTGCAGGGGCATGGATCTCCTCCACGCGTATTGCGGGAAGCAGGTGCCGACGAAGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcgcatcg                                                                                                                                          >  1:267829/1‑251 (MQ=255)                                                                                                            cacTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCGCTCATCCACGCGGATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCCCTTTTCAACACCCCTAatcgcatagcc                                                                                                                                        >  2:119467/1‑247 (MQ=255)                                                                                                             acTGGCGGAAAACCTGGTTGGCGAGAACAACGATAATACTGTTGGCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTc                                                                                                                                       <  2:14869/251‑1 (MQ=255)                                                                                                                 gCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGTTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTAt                                                                                                                                   <  1:92857/251‑1 (MQ=255)                                                                                                                         cctggattgccggaACAACGATAATACTGTTGTCGATACCAACGGTGAGCGTCGGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTTCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGGTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACGCACGTTTCAACTCCCCTAATCGCATCGCTCGTATCCGCTCAc                                                                                                                           <  2:102504/239‑1 (MQ=255)                                                                                                                            ggTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCAcc                                                                                                                          <  1:236291/249‑1 (MQ=255)                                                                                                                                                     gttgtCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGc                                                                                                <  1:187843/250‑1 (MQ=255)                                                                                                                                                        gtCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGcc                                                                                                                                                                                                                                     >  2:96459/1‑114 (MQ=255)                                                                                                                                                              aCCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGGTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTAGCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATttcatt                                                                                       <  2:65544/250‑1 (MQ=255)                                                                                                                                                                    ggTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCCTGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTCCTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCCTTCTGATGc                                                                                >  2:35136/1‑251 (MQ=255)                                                                                                                                                                     gTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCt                                                                               <  1:139674/251‑1 (MQ=255)                                                                                                                                                                                 ggACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCCAATCCACGCGTATTGCGGGAAGCAGGTGCCGCCGACGCAGAAATGCTGGTTGCTGTAACAAGTTCAGCTGAAACCAATATGGGTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACAAGACTACGGGCGCGATGCCGCTACGCTATTTCATTCAGATGCTGGGCCGATTggt                                                                   >  2:21103/1‑249 (MQ=255)                                                                                                                                                                                    ccTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGatcat                                                                <  1:14686/251‑1 (MQ=255)                                                                                                                                                                                               aaaTTTGCCCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCAGCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCTGATGCTGTGCCGATTGATCATCTGATCGCAc                                                      <  1:21825/250‑1 (MQ=255)                                                                                                                                                                                                 aTTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCAcc                                                     >  1:64560/1‑249 (MQ=255)                                                                                                                                                                                                                   gCAGGGGCATGGCTCTCATCCACGTGTATTGCGGTAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCATTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATACGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGGTGGTTATCGata                                 <  2:267829/251‑1 (MQ=255)                                                                                                                                                                                                                     aGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGataat                               >  1:248056/1‑251 (MQ=255)                                                                                                                                                                                                                        ggCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATAtt                            >  1:3019/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATcc        >  1:108433/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTCCTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGCCTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGACGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGgcgc  >  1:153178/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGC  >  NC_000913/3436451‑3436942

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