Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,034,547	A→G	E468G (GAG→GGG)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,034,547	0	A	G	96.5%	87.7 / NA	29	E468G (GAG→GGG)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/12);  minor base T (0/1);  total (16/13)
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

ACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGAGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAG  >  NC_000913/4034409‑4034686                                                                                                                                           |                                                                                                                                            aCTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTc                                                                                                                                         >  1:325574/1‑143 (MQ=255)       cTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACa                                                                                                                                   >  1:561967/1‑143 (MQ=255)        tGCCTTTGCGTCTGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACAt                                                                                                                                  <  1:138531/143‑1 (MQ=255)         tccattGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACAtt                                                                                                                                 <  1:424449/139‑1 (MQ=255)         gCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACAtt                                                                                                                                 >  1:630031/1‑143 (MQ=255)           ctctGCGTCCGGAGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGc                                                                                                                               <  1:493584/140‑1 (MQ=255)                 gTCGTTTGTGGCGACATTGAATGACCTGGGGCCAGGGCTCGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGAGCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGc                                                                                                                         <  1:238237/143‑1 (MQ=255)                                     aTAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCCTTCCATTGCTGGTGCTCTTTACCCCCACTATCggg                                                                                                     >  2:323541/1‑143 (MQ=255)                                      tAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGGTCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGc                                                                                                    <  1:472397/143‑1 (MQ=255)                                                 cAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGATCTTTACCCCGACTTTCTGGCGTGAATGATgg                                                                                         >  1:173298/1‑143 (MQ=255)                                                  aGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGa                                                                                        <  1:257371/143‑1 (MQ=255)                                                       ttGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATGCTGATTGCCAACGTGGTGTTTGGTCGTCTCGGGGTCTTTAGGTTGCTGGTGCTGTTTGCCCCGGCTGTCTGGCGTGAAGGGTGGAGtgag                                                                                   >  1:437618/1‑140 (MQ=255)                                                                ttGCTGGTAACTGGACCAGTAGGAGCCCGGGGGCTAAAGGGATGCTGATTGCCAACATGCTGTTCGGTCGGGTCGGGGTCTGTACCTTGCTGGTGCGCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGaaaa                                                                          <  1:135370/143‑1 (MQ=255)                                                                 tGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAAc                                                                         <  1:324623/143‑1 (MQ=255)                                                                    tGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGACTTTACAGTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGGAGACAtt                                                                      >  1:622036/1‑143 (MQ=255)                                                                           tttACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGGGAATGATGGAGTAAACCGTCAAAGCAGGAATTGtt                                                               >  2:44658/1‑143 (MQ=255)                                                                             tACCAGTATGGACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCtttt                                                             <  1:593785/143‑1 (MQ=255)                                                                                 aGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGCATGATGGAGTAATACGTGAACACATTCATTCGTTCCTca                                                         >  1:306058/1‑143 (MQ=255)                                                                                    aTGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTCCATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTccaca                                                      >  1:736257/1‑143 (MQ=255)                                                                                    aTGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGCGTAATACATGAAAACATTAATTTTTTTCTcaaca                                                      >  1:401823/1‑143 (MQ=255)                                                                                                    aaTGGAGCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAAcgc                                      <  1:630572/143‑1 (MQ=255)                                                                                                     aTGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCTGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTGAACAAGGGACGGACAAAcgcg                                     >  2:706470/1‑143 (MQ=255)                                                                                                                ttGTCAACATCCTGTTTGGTCGTCTCGGGGTCTTTGAATTGCTTGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGcct                          <  2:561967/143‑1 (MQ=255)                                                                                                                  gCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGcctcc                        >  1:413444/1‑143 (MQ=255)                                                                                                                    cAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGcctcctc                      >  2:626383/1‑142 (MQ=255)                                                                                                                              gttGGTCGGCGCGTGGTCGTTACTTTGCTGGTGCTGTTTGCCCCGACTTTCTGGCGTGAATGAGGGAGCAATACGTGCGAACATTAATTCTTTTCTCCACAAAGGACGGCCATACGCGCGAGATTGCCTCCTGCCTGGCTTCg            <  2:27563/142‑1 (MQ=255)                                                                                                                                      tctcGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATAAGAGAAAAAATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCATCGGAACTGaa    >  2:758386/1‑143 (MQ=255)                                                                                                                                       ctcGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGaaa   <  1:518163/143‑1 (MQ=255)                                                                                                                                        tcGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTGCGGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAg  >  2:439734/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                            ACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGAGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAG  >  NC_000913/4034409‑4034686

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