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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,034,5470AG88.0% 72.9 / ‑1.6 33E468G (GAG→GGG) trkHK(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base G (19/9);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGAGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAG  >  NC_000913/4034408‑4034686
                                                                                                                                           |                                                                                                                                           
gACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGt                                                                                                                                          >  1:219549/1‑143 (MQ=255)
   ttttCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCGGTTTGGCCGTCTCCGGGTcgg                                                                                                                                       >  2:340594/1‑141 (MQ=255)
   ttttCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGGTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCGAACATGCTGTTTGGTCGGCTCGGGGTGtt                                                                                                                                       >  2:459848/1‑143 (MQ=255)
       cTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGACTGCCAACATGCGCTTTGGTCGTCTCGGGGtgtgcata                                                                                                                                   >  2:509706/1‑136 (MQ=255)
                                                  aGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGACAACATGCGGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGa                                                                                        <  1:517157/143‑1 (MQ=255)
                                                     gCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGTGTGAATGATGGAGTa                                                                                     >  1:73667/1‑143 (MQ=255)
                                                          gCGTGGTTGCTGATAACCTTACCAGTATGATCCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACg                                                                                <  1:325336/143‑1 (MQ=255)
                                                             tGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGCCTTTCTGGCGTGAATGATGGACTAGTActaga                                                                             >  2:371440/1‑139 (MQ=255)
                                                              ggTTGCTGTTGCCTTTATCAGTGTGTGCCCGGTGGCTAGGTTGAGCCTGATTGCCATCATGCTGTGTGGTCGTTTCGGGGTCTTTACGTTGCTGGTGCTCTTTTCCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGaa                                                                            <  2:181157/143‑1 (MQ=255)
                                                                  gCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATCGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCGTTGCCCCGACTGTCTGGCGTCAATGCTGGAGCGATACGTgggaccg                                                                        >  2:426660/1‑137 (MQ=255)
                                                                  gCTGATAACTTTACCAGGATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACa                                                                        <  1:259166/143‑1 (MQ=255)
                                                                          cTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGGCTGTACATTGCTGGTGCTCTTAACACCGACTCTCTGGCGTGAATGATGGAGAAGTACGCAAAAACATAAAGTca                                                                >  2:355903/1‑142 (MQ=255)
                                                                           tttACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACCTTGCTGTTGCTCTTTACCCCGACTTTCTGACGTGAATGATGGAGTACTACGTGAAAACATTAATTCtt                                                               >  1:33450/1‑143 (MQ=255)
                                                                           tttACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGCCGCGAATGATGGTTGAGTACGTGCAAACATTAATTCtt                                                               >  1:18526/1‑143 (MQ=255)
                                                                             tACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAATATACTATTGGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGCCTTTCTGGCGTGAATGATGGAGTAATACGTGCAAGCATTAATTCtttt                                                             >  1:548987/1‑143 (MQ=255)
                                                                                     tGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTGGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTcaacaa                                                     <  1:428358/143‑1 (MQ=255)
                                                                                      gAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCGCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCCACGAg                                                    >  2:312332/1‑143 (MQ=255)
                                                                                          ccGGTTGCTAAATGCATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGAGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGggac                                                <  1:160429/143‑1 (MQ=255)
                                                                                           cGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGggacg                                               >  1:555633/1‑143 (MQ=255)
                                                                                           cGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGCAGTAATAGGAGACACAATGAATATTTTTCTCAAGAAGggacg                                               >  2:142775/1‑143 (MQ=255)
                                                                                                     ttGGATCCTGGTTGCCATCATGCTGTTTGGTCGTCTCGGGGTGTTTACATTGCTGGTGCGCTGTGCCCCGACTTTCTGGCGTGACTGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACCAGGGACGGACAAAcgcg                                     <  2:286277/142‑1 (MQ=255)
                                                                                                     aTGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACCAGGGACGGACAAAcgcc                                     >  2:130220/1‑142 (MQ=255)
                                                                                                          tCCTGATTGCCAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCgaga                                >  1:557517/1‑143 (MQ=255)
                                                                                                                    cAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTTGTGCTCTTTACCCCGACTTTTTCGAGTCCAATAAGCAGAAATTCGTGGTCCCATTAACTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTa                      <  1:576628/143‑1 (MQ=255)
                                                                                                                    cAACATGCTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATCCTCACCACCACAATGGACGGACAAACGCGCGAGATTGcctcctc                      >  2:406760/1‑142 (MQ=255)
                                                                                                                       cATGCTTTTTGGTCGTCTCGCCGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCt                   <  2:106182/143‑1 (MQ=255)
                                                                                                                        attcTGTTTGGGCGTCTCGGGCTGTGTCCATTGCTGGGGCTTTTTATACTAAGTATCTGGCGTGAATGATGGAGTGGTACGTGAGAACATGAATTCTTGCCTCAACTCGGGATGGACAAACGGGCGAGATTGCCTCCTACCTg                  <  2:229395/140‑1 (MQ=255)
                                                                                                                           cTGTTTGGTCGTCTCGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAAAACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCt               >  1:302826/1‑143 (MQ=255)
                                                                                                                              tttGGTCGGCTCGAGGTCTTTACATTGCTGGTGCGCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCg            <  2:508125/143‑1 (MQ=255)
                                                                                                                              ttcggTCGTCCCCGGGTCTTTACATTGCTGGTGCTCTTCACCACGACTGTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGTACAAACGCGCGAGATTGCCTCCTACCTGGCTTCg            <  1:157833/140‑1 (MQ=255)
                                                                                                                               tttgtcgttTGGAGGTCTTTACATTGCTGTTGCTCTTTCCCCCGATTTTCTGGCGTGAATGATGTAGTAATACGTGAAAACATTAATTCTTTCCTCAACAAGTGTCGGATAAACGCTCGAGATTGACTCCTTCCTGGCTTCgn           <  2:73667/140‑2 (MQ=255)
                                                                                                                                      tgtcGAGGTCGTTACCTTGCTGGTGCTCTTTACCCCGACTTGCTGGCGTGAATGATGGAGTAATACGTGGAAACATTAATTCTTTTCTCAACCAGGGTCGCACGAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGaa    <  2:316704/141‑1 (MQ=255)
                                                                                                                                        tcGGGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGCACTGaaat  >  1:2857/1‑142 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GACTTTTCTGCCTTTGCGTCGGTTGTTGCGACATTGAATAACCTGGGGCCAGGGCTTGGCGTGGTTGCTGATAACTTTACCAGTATGAACCCGGTGGCTAAATGGATCCTGATTGCCAACATGCTGTTTGGTCGTCTCGAGGTCTTTACATTGCTGGTGCTCTTTACCCCGACTTTCTGGCGTGAATGATGGAGTAATACGTGAAAACATTAATTCTTTTCTCAACAAGGGACGGACAAACGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAG  >  NC_000913/4034408‑4034686

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