Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,465,033	C→A	22.8%	F355L (TTC→TTA)	trkH‑I →	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,465,033	0	C	A	22.8%	56.1 / 2.2	25	F355L (TTC→TTA)	trkH‑I	potassium transporter
Reads supporting (aligned to +/- strand):  ref base C (8/11);  new base A (4/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTAACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTC  >  NC_002947/4464906‑4465162                                                                                                                                |                                                                                                                                  tagtaACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTACGTTTCCa                                                                                                                           >  1:20303/1‑136 (MQ=255) tagtaACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTACGTTTCCa                                                                                                                           >  1:388134/1‑136 (MQ=255)                   ttGCCGTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTc                                                                                                                             >  2:271241/1‑116 (MQ=255)                   ttGCCGTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTc                                                                                                                             <  1:271241/116‑1 (MQ=255)                        gTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTTAAAATATTACGTTTACAGGTGGCCGCGGCCTTGCTTGAca                                                                                                    >  1:159774/1‑136 (MQ=255)                              gATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGcc                                                                                              <  2:90864/136‑1 (MQ=255)                                         ctGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTg                                                                                   >  1:347575/1‑136 (MQ=255)                                           ggggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGat                                                                                 >  2:333099/1‑136 (MQ=255)                                           ggggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCCGCGGTTTGAAAATAATACGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGAACCCAGAAGCAGTTGat                                                                                 >  1:481566/1‑136 (MQ=255)                                             ggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGatcc                                                                               >  1:443424/1‑136 (MQ=255)                                             ggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGatcc                                                                               >  1:355801/1‑136 (MQ=255)                                                     cGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACgg                                                                       >  2:329680/1‑136 (MQ=255)                                                              gCTGTTCTTTTACCTGACTTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGAtc                                                              <  2:416337/136‑1 (MQ=255)                                                                     ttttACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCgaaga                                                       >  2:345349/1‑136 (MQ=255)                                                                                  tcgtcgGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCa                                          <  2:308004/136‑1 (MQ=255)                                                                                  tcgtcgGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCAACCACGGGCAACGATCTCGAAGAAATACAATGGGCa                                          <  1:64938/136‑1 (MQ=255)                                                                                    gtcgGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCAATCCACGGGCAACGATCTCGAAGAAATACAATGGGCAcc                                        <  2:10345/136‑1 (MQ=255)                                                                                             tGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATg                               <  2:349604/136‑1 (MQ=255)                                                                                                         aCCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTAc                   <  2:355801/136‑1 (MQ=255)                                                                                                         aCCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTAc                   <  2:390515/136‑1 (MQ=255)                                                                                                         aCCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTAc                   <  2:443424/136‑1 (MQ=255)                                                                                                         aCCGCCGGCGGTTTGAAAATATTACGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTAc                   <  2:20303/136‑1 (MQ=255)                                                                                                         aCCGCCGGCGGTTTGAAAATATTACGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTAc                   <  2:388134/136‑1 (MQ=255)                                                                                                                 cggTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCActg           >  2:458830/1‑136 (MQ=255)                                                                                                                          aaTATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttc  <  1:329680/136‑1 (MQ=255)                                                                                                                                |                                                                                                                                  TAGTAACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTC  >  NC_002947/4464906‑4465162

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