Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,595	C→A	P151T (CCT→ACT)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,595	0	C	A	97.0%	67.6 / NA	34	P151T (CCT→ACT)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (21/11);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACGGGAGCCACTACGCTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCT  >  NC_000913/4033476‑4033719                                                                                                                        |                                                                                                                             cTACGGGAGCCACTACGCTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTACACTCg                                                                                                                     >  2:328787/1‑129 (MQ=255)          ccgctgcGCTGGTGGGGCTGGATGCTGTCCCGCACTCCATCGTCGGTGGTCGCCGGATGCTGGGATGGTTTGGCGGGATGGGGTGCATCGTGTTGGCGGTTGCGATACTGGCTATCCTAGGCgtgtgtg                                                                                                            <  2:564863/123‑1 (MQ=255)              tgcGCTGGTGGGGCTGGATGCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGAt                                                                                                        <  2:127402/127‑1 (MQ=255)                  cTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTTGCGGGTTTTGTATTTTTTTTTTTTCGGTTTCTATACTGACTATCCGCGGTGTGGGTGGGATGCAg                                                                                                    <  1:133216/129‑1 (MQ=255)                   tggtggGGCTGTATTCGCTCCCTCACGCCATCCTCTGTTATCGCCAGGTGCTGCAATGGTTTGGCGGGAGGGGGATCATTTTTTTGGCGGTTGTGATACTGACTATGCTCGGCGGGGGTGGGATGCAGc                                                                                                   <  1:394005/129‑1 (MQ=255)                     gtggGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGTTTCGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCGTCCGCGAGGGTGGGAGGCAActa                                                                                                 >  2:372999/1‑128 (MQ=255)                           cTGGATTCGCTCCCTCCCCCCCTCCTCTTTTATCGCCAGATACTAAAATGGATTGATGGGATAGAGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGAGTGGGTGGGATGCAGCTCTAtggg                                                                                           >  1:577704/1‑126 (MQ=255)                             gggttCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCgcgc                                                                                         <  2:36557/126‑1 (MQ=255)                               aTTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATACTAAAAAGGATTGGCGGGATGGAGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAg                                                                                       >  1:576955/1‑129 (MQ=255)                                 ccGCTACCTCTCGCGATCCTCTTTGCTCGCCAGATGCTGCAATGGGTTGGCGGGATGGGGATCATCGTGTGGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGaa                                                                                     <  2:532451/128‑1 (MQ=255)                                    cTCCCGCACGCCATCTTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATTGTTTTGTTTGTTTTTATACTGACTATCCTCGGCGTAGGTGGGATGCAGCTCTAGCGCGCAGAAATg                                                                                  <  1:66131/129‑1 (MQ=255)                                        cTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGAAATAAAGGAATGGGGATAATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATgcccg                                                                              >  1:146135/1‑129 (MQ=255)                                                           tagCCAGATGCTGCAATCGACATGCGTGATGGGGAGCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGAGTGGGTCGGATGCAGCACTATCGCGCAGAAAAGACCGGCCCGCGGAAAGATAACaa                                                           >  2:256623/3‑129 (MQ=255)                                                               cAGATGCTGCAATGGTTGGGCGGTATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATg                                                       >  2:365322/1‑129 (MQ=255)                                                                  aTGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTACCGCGCAGAAAAGCCCGGCCCGCTGAACGATAACAAAATgcgc                                                    >  2:240069/1‑129 (MQ=255)                                                                  aTGCTGCAATGGTTTGGCGGGATGGGGATCATCGTATTAAAAATTGAAAAACAAACAATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATgggc                                                    >  1:611888/1‑129 (MQ=255)                                                                          aTGGTTTGGCGGGATGGGGATCATCGTGGTGGCGGTGGCGATACTGACTAGCCTCCGCGTGGGTGGGATCCAGCTCTAACGCGCAGATATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAAt                                            <  2:87919/129‑1 (MQ=255)                                                                           tGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCACGCTGAAAGATAACAAAATGAGCCCGCGAAtt                                           >  2:481974/1‑129 (MQ=255)                                                                            ggTTTGGCGGGCTGGGGATCATCGTGTTGGCGGTAGAGAAACTGAAAAACCTCGACGTAGGTGGGATGCAGCTCTATAGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAAttc                                          >  1:344420/1‑128 (MQ=255)                                                                            ggTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTg                                          >  2:3110/1‑129 (MQ=255)                                                                            ggTTTGGCGGGATCGGGTTCATCGTGTTGGCGGTTGCGATACTGACTATCCTCGACGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGCATTg                                          >  1:127133/1‑129 (MQ=255)                                                                                     ggATGGGGATCATCGTGTTGGCGGTTGCGATACTGACTATCATAGACGAGAGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACgg                                 >  1:538960/1‑129 (MQ=255)                                                                                              tcatcGTGTTTGCGGTTGCGTTACTGACTATCCTCGGCGGGGGTGGGATGCAGCTCTATCGCGCAGAAAGGCCCGGCCCGCTGAAAGATAACAAAATGTGCCCGCGAATTGCGGAAACGGCGAAAAccc                        <  2:606220/129‑1 (MQ=255)                                                                                              tcatcGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAAccc                        <  2:156849/129‑1 (MQ=255)                                                                                               catcGAGGTGACGGCTGCGATACTGATTATCCTCGGCGTGGGTAGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCAGCGAAATCCGGAAACGGCGAAAACCCt                       >  2:118202/1‑129 (MQ=255)                                                                                                atcGTGTTGGCGGTTGCGATACTGACTATCCTCAACGTGGGAGAGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCtg                      >  1:491114/1‑129 (MQ=255)                                                                                                atcGTGTTGGCGGTTGCGATACTGACTATCCTAGGCGTGGAAAGGATAAAGCTCTATCGCGCAGAAATGCACGGCCAGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCtg                      >  1:265168/1‑129 (MQ=255)                                                                                                atcGTGTTGGCGGTTGCGATACTGACTATCATAGGAGAAGAAGGAATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTTCGGAAACGGCGAAAACCCtg                      >  1:145775/1‑129 (MQ=255)                                                                                                  cGTGTTGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATAGCACAGAAATGAACAGGCCGGTGAAGGATAACAAAGTGCGCGCGGGAGTTGCGGAGGCGGCGAAAGCACtgag                    >  1:427795/1‑127 (MQ=255)                                                                                                   gtgtTGGCGGTTGCGATACTGACTCTCCTCGACGTGGGTGAAATGCAGCTCTATAGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTgg                   >  1:449743/1‑129 (MQ=255)                                                                                                      ttGGCGGTTGCGATACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCTCGCTGAATGATATTTTTTTTTGTTCGTGAATTGAGGAAACGGCGAAAACCCGGTGGTTg                <  1:328787/129‑1 (MQ=255)                                                                                                         gCGGTTGCGATACTGACTATCCTCGGCGTGGGTGAGAAACAAATAAAAAAAAAAAAAAAACACAGCCAGATGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACGCTGTGGTTGAtt             >  1:36584/1‑129 (MQ=255)                                                                                                                  aTACTGACTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGACTTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTg    >  2:385631/1‑129 (MQ=255)                                                                                                                    aCTGACTATCCTCGGCGTGGGTGGGATGCACCTCTATCGAGAAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGTTTGCTAATTGCGGAAATGGCGAAAACCCTGTGGTTGATTTATGTCTTGCt  <  1:574290/129‑1 (MQ=255)                                                                                                                        |                                                                                                                             CTACGGGAGCCACTACGCTGGTGGGGCTGGATTCGCTCCCTCACGCCATCCTCTTTTATCGCCAGATGCTGCAATGGTTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCT  >  NC_000913/4033476‑4033719

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