Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,034,217	G→A	R358H (CGC→CAC)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,034,217	0	G	A	100.0%	72.8 / NA	23	R358H (CGC→CAC)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/21);  total (2/21)

CGCTGGTGGTTATTTGTACCCTCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGG  >  NC_000913/4033992‑4034454                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              cGCTGGTGGTTATTTGTACCCCGGTACTGGGTTTTCATAATGACTACGTTTGGGCGTTGATGACATTTAACCGGGCGTTGTTCCAGGGGGTGTCGATGGCTACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTGTATGTTCAGCCTTTATCGGCGGTTGCGCCGGTTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAgggt                                                                                                                                                                                                                      <  2:294168/251‑2 (MQ=255)     ggtgttTATTTGTCCCCTCGTACTTTGGTTTCATAATGTCTACAGTCCGGCTCGGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAAcc                                                                                                                                                                                                                  <  2:304577/251‑1 (MQ=255)                                                           gataaCAATTAACCCGGCGTTTTCCCGGGTGGTGTCGATGCCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTGTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTTAGCTGGGGa                                                                                                                                                            <  2:29277/247‑1 (MQ=255)                                                                     agccAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACt                                                                                                                                                  <  1:308277/249‑1 (MQ=255)                                                                        cAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGcc                                                                                                                                               <  1:33827/251‑1 (MQ=255)                                                                          ggtgTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCg                                                                                                                                              <  2:245307/247‑1 (MQ=255)                                                                                     ggttggtgTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTc                                                                                                                                  <  1:119141/248‑1 (MQ=255)                                                                                        gggggtCGAGGGCGGCAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGaa                                                                                                                               <  2:288587/247‑1 (MQ=255)                                                                                                        aaCCGCCCGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGAtttt                                                                                                               <  2:401348/251‑1 (MQ=255)                                                                                                         aCCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGAttttt                                                                                                              <  1:138236/251‑1 (MQ=255)                                                                                                                   ttcaaaCTGACAGCTTTGCCCGATGGCCGCTTTTTTGGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGGTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCTTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTTGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGATCGTATACTCGAAGCCGTTTGGGGTTTTTTCTCCGCCTa                                                                                                     <  2:385906/246‑1 (MQ=255)                                                                                                                            aatgCATTGCCCGCTGGCCGCTCTTTTTGCCTGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCACCGCCTATGCATTGgtg                                                                                           <  2:285938/248‑1 (MQ=255)                                                                                                                               gCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGttt                                                                                        <  1:338415/251‑1 (MQ=255)                                                                                                                                   tGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGCGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTAtt                                                                                     <  1:267778/250‑1 (MQ=255)                                                                                                                                                         ccGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCGGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCg                                                               <  2:301026/250‑1 (MQ=255)                                                                                                                                                            gTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGTAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCAc                                                           <  2:238718/251‑1 (MQ=255)                                                                                                                                                                 gCTTTGATGTTCAGCCTTTATCGGCGGTTGCGGCGGGTCAACGGGCGGGGGCCTGAAGGTGATCCACAGCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCGCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACggg                                                        <  1:25282/249‑1 (MQ=255)                                                                                                                                                                       aTGTTCAGACTTTATCGGCGGTTGCGCCAGGTCAACGGGCGGTGGCCTGAAGGTGATCTACATCCTGCTGCTGGTTAAGCAGGGGAACCGTGATCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGTGTAATCGTGCACTGCCGGAACGTATCCTCGAAGCCGTGTGGGGATGTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCAGCAGGGGGGTGGATg                                                <  2:162786/251‑1 (MQ=255)                                                                                                                                                                                      cggcggTTGCGGCGGGTCAACGGGCGGTGGCCTGAAGGTGAGCCACATCCTGCTGCTGTTTAAGCGGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATTCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCttt                                  <  2:149493/250‑1 (MQ=255)                                                                                                                                                                                           gTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGAGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGCACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCg                            >  2:49509/1‑251 (MQ=255)                                                                                                                                                                                              gcgcCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCTCTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGgt                          <  1:7450/250‑1 (MQ=255)                                                                                                                                                                                                                gTGGCCTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCGCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAAc       <  1:49509/251‑1 (MQ=255)                                                                                                                                                                                                                     cTGAAGGTGATCCACATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGACGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTggg  >  1:204850/1‑251 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              CGCTGGTGGTTATTTGTACCCTCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTCTTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCACTGCCGGAACGTATCCTCGAAGCCGTTTGGGGATTTTTCTCCGCCTATGCATTGGTGTTTATTGTCAGTATGCTGGCGATTATCGCCACGGGCGTGGATGACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGG  >  NC_000913/4033992‑4034454

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