Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,464,862	C→A	21.1%	A298A (GCC→GCA)	trkH‑I →	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,464,862	0	C	A	21.1%	46.4 / 4.8	19	A298A (GCC→GCA)	trkH‑I	potassium transporter
Reads supporting (aligned to +/- strand):  ref base C (6/9);  new base A (2/2);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GTATGTCGCCACGCTGCGTGGCCATCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCG  >  NC_002947/4464733‑4464959                                                                                                                                  |                                                                                                  gTATGTCGCCACGCTGCGTGGCCATCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCtggtgg                                                                                             >  2:38562/1‑136 (MQ=255)         ccACGCTGCGTGGCCATCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGtt                                                                                     <  1:148043/136‑1 (MQ=255)                    ggCCATCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTg                                                                          >  1:4069/1‑136 (MQ=255)                    ggCCATCGAGGCATCCTGCTCAAGGACCATCAACTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTg                                                                          >  1:202895/1‑136 (MQ=255)                     gCCATCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGc                                                                         <  1:114191/136‑1 (MQ=255)                         tCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTg                                                                     <  1:147361/136‑1 (MQ=255)                                   cTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCt                                                           >  2:261064/1‑136 (MQ=255)                                     gCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCGCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCg                                                         <  1:326491/136‑1 (MQ=255)                                                 tCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGgtagt                                                    <  1:272035/129‑1 (MQ=255)                                                 tCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGgtagt                                                    >  2:272035/1‑129 (MQ=255)                                                             tttCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCtggt                                                                                               <  1:194634/74‑1 (MQ=255)                                                             tttCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCtggt                                                                                               >  2:194634/1‑74 (MQ=255)                                                                       tttCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGAtt                       <  1:64955/136‑1 (MQ=255)                                                                       tttCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGAtt                       <  2:195233/136‑1 (MQ=255)                                                                          ctcGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCATGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGAtt                       <  1:372258/133‑1 (MQ=255)                                                                          ctcGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCATGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGAtt                       <  1:381582/133‑1 (MQ=255)                                                                          ctcGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCATGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGAtt                       >  2:372258/1‑133 (MQ=255)                                                                          ctcGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCATGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGAtt                       >  2:381582/1‑133 (MQ=255)                                                                                            aCCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCg  <  1:178700/136‑1 (MQ=255)                                                                                                                                  |                                                                                                  GTATGTCGCCACGCTGCGTGGCCATCGAGGCATCCTGCTCAAGGACCATCAAGTTCGCGGTTTCATTGGTTTTCTCGTCTTCGTCTGGCTCACCGTCGGAACCTGGCTCTGCCTGCACAGCAATCTGGCCTGGTGGGATGCGTTTCGTATCGTTGCGGTGAATGTGACCTCGGTAGTAACCACAACAGGCATTGCCGTGGGTGATTACACACTCTGGGGCAGCTTCG  >  NC_002947/4464733‑4464959

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