Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 74,984 A→C 6.8% I409S (ATT→AGT)  trkA ← TrK potassium transporter NAD‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294774,9840AC 6.8% 142.9 / 2.1 51I409S (ATT→AGT) trkATrK potassium transporter NAD‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/18);  new base C (2/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCACGACGGACGTGC  >  NC_002947/74835‑75133
                                                                                                                                                     |                                                                                                                                                     
ccTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCaa                                                                                                                                                       <  1:1043401/150‑1 (MQ=255)
 cTAGAAGAAACGCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTCTCGTCGTGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAAt                                                                                                                                                      >  1:174299/1‑150 (MQ=255)
       gaaACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAgtggt                                                                                                                                                >  1:1009674/1‑150 (MQ=255)
           cTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCACTATTGGTGcga                                                                                                                                            <  1:507877/150‑3 (MQ=255)
            tCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCAGCGCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCg                                                                                                                                           >  1:63247/1‑150 (MQ=255)
            tCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCACTATTGGTGcga                                                                                                                                            <  1:5556/149‑3 (MQ=255)
            tCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCACTATTGGTGcga                                                                                                                                            >  2:5556/1‑147 (MQ=255)
              aaGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGt                                                                                                                                         >  1:636751/1‑150 (MQ=255)
                              gCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCt                                                                                                                         <  2:576947/150‑1 (MQ=255)
                               cTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCATGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTtctc                                                                                                                        >  1:1200314/1‑150 (MQ=255)
                                           gcgAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTc                                                                                                                        >  1:921307/1‑138 (MQ=255)
                                           gcgAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTc                                                                                                                        <  2:921307/138‑1 (MQ=255)
                                                     tttttATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTg                                                                                                  >  2:883426/1‑150 (MQ=255)
                                                      ttttATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACACCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGc                                                                                                 <  1:927565/150‑1 (MQ=255)
                                                                aaCGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGa                                                                                       >  1:195305/1‑150 (MQ=255)
                                                                aaCGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCGGTCTCGACTTCGa                                                                                       >  1:548994/1‑150 (MQ=255)
                                                                   gAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATc                                                                                    >  2:168959/1‑150 (MQ=255)
                                                                    aaGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGt                                                                                                                                         >  1:558041/1‑96 (MQ=255)
                                                                    aaGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGt                                                                                                                                         <  2:558041/96‑1 (MQ=255)
                                                                         gATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATATCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCAtg                                                                              >  1:734966/1‑150 (MQ=255)
                                                                           tCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTc                                                                                               <  1:655916/131‑1 (MQ=255)
                                                                           tCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTc                                                                                               >  2:655916/1‑131 (MQ=255)
                                                                                  gTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGc                                                                     <  2:1009674/150‑1 (MQ=255)
                                                                                   tCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCAt                                                                               <  2:1008873/139‑1 (MQ=255)
                                                                                   tCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCAt                                                                               >  1:1008873/1‑139 (MQ=255)
                                                                                     cccTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTc                                                                  <  1:224781/150‑1 (MQ=255)
                                                                                       cTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGa                                                                >  1:588675/1‑150 (MQ=255)
                                                                                       cTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGa                                                                >  1:500012/1‑150 (MQ=255)
                                                                                       cTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGa                                                                <  2:63247/150‑1 (MQ=255)
                                                                                               aTCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCg                                                        >  2:417113/1‑150 (MQ=255)
                                                                                                tCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGCTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCCTGTGCCACCGCCTCGACCGCTTCgg                                                       >  2:28534/1‑150 (MQ=255)
                                                                                                 cATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGc                                                      <  1:883426/150‑1 (MQ=255)
                                                                                                      tgtCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCAc                                                 >  1:216477/1‑150 (MQ=255)
                                                                                                            cctGGGCCATCATCACTTCTTCCTCACGGATCATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACGTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCgcggc                                           <  2:216477/148‑1 (MQ=255)
                                                                                                              tGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGgcgc                                         >  2:392182/1‑150 (MQ=255)
                                                                                                                     tcatcaCTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGt                                  >  2:91190/1‑150 (MQ=255)
                                                                                                                        tcaCTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATc                                                              <  1:806685/119‑1 (MQ=255)
                                                                                                                        tcaCTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCAGGTGCCACCGCCTCGATc                                                              >  2:806685/1‑119 (MQ=255)
                                                                                                                            ttcttcGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACt                           >  1:592367/1‑150 (MQ=255)
                                                                                                                            ttcttcGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACt                           >  1:519779/1‑150 (MQ=255)
                                                                                                                               ttcGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGAc                        <  2:588675/150‑1 (MQ=255)
                                                                                                                                 cGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGa                      <  1:168959/150‑1 (MQ=255)
                                                                                                                                      cGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCg                 >  2:186560/1‑150 (MQ=255)
                                                                                                                                          tgatCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCac             <  1:392182/150‑1 (MQ=255)
                                                                                                                                           gatCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCacg            >  2:936111/1‑150 (MQ=255)
                                                                                                                                           ccgggCACCACTATTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCacg            >  2:1198262/5‑150 (MQ=255)
                                                                                                                                            atCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCacga           >  2:588230/1‑150 (MQ=255)
                                                                                                                                            atCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCacga           <  1:710730/150‑1 (MQ=255)
                                                                                                                                                cACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCACgacgga       >  2:605776/1‑150 (MQ=255)
                                                                                                                                                 aCCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGcaca                             >  2:789973/1‑127 (MQ=255)
                                                                                                                                                 aCCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGcaca                             <  1:789973/127‑1 (MQ=255)
                                                                                                                                                 aCCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCACgacggac      <  1:130664/150‑1 (MQ=255)
                                                                                                                                                    aaTAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCg                                                        >  1:1145147/1‑97 (MQ=255)
                                                                                                                                                    aaTAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCg                                                        <  2:1145147/97‑1 (MQ=255)
                                                                                                                                                     aTAGTGGTGCCAGGTGGACGCGTGATGTCATCAATTCTCTTTCCTTCCACCTTGCTCTACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCACGACGGACGTgc  <  2:734966/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTAGAAGAAACTCAAGCCGACGTGGAACAGCTTCTCCACATCGCGAATATGCTTTTTATCCACAACGAAGAGGATCACATGGTCCCCTGATTCGATCATGGTGTCGTCATGGGCAATCATCACTTCTTCGTCACGGATGATCGCACCAATAGTGGTGCCCGGTGGCAGCGTGATGTCCTCGATTGCCTTGCCGACCACCTTGCTCGACTTCGAATCGCCATGTGCCACCGCCTCGATCGCTTCGGCCGCACCGCGGCGCAGCGAGTGCACACTGACGATGTCGCCACGACGGACGTGC  >  NC_002947/74835‑75133

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