Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,754,270 A→G 16.4% V413A (GTG→GCG)  trkA ← NAD‑binding component of TrK potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,754,2700AG16.4% 164.2 / 23.5 73V413A (GTG→GCG) trkANAD‑binding component of TrK potassium transporter
Reads supporting (aligned to +/- strand):  ref base A (36/25);  new base G (8/4);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CATAATTACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTGCT  >  minE/2754200‑2754340
                                                                      |                                                                      
cATAATTACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTacca                                                                        >  1:1803417/1‑71 (MQ=255)
    aTTACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTaccacc                                                                      >  1:1869572/1‑69 (MQ=255)
    aTTACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCg                                                                     >  1:2100340/1‑70 (MQ=255)
     ttACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGc                                                                    >  1:2165132/1‑70 (MQ=255)
     ttACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCt                                                                   >  1:3147507/1‑71 (MQ=255)
      tACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTc                                                                  >  1:2868800/1‑71 (MQ=255)
       aCGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCt                                                                   >  1:125415/1‑69 (MQ=255)
       aCGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTc                                                                  >  1:1858950/1‑70 (MQ=255)
       aCGTGATAGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTc                                                                  >  1:3650954/1‑70 (MQ=255)
        cGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTcc                                                                 <  1:3819172/70‑1 (MQ=255)
           gATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTCCaat                                                              >  1:3284685/1‑70 (MQ=255)
             tCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTCCaataat                                                           >  1:1688976/1‑71 (MQ=255)
             tCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTCCaataa                                                            >  1:2802051/1‑70 (MQ=255)
             tCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCaataat                                                           >  1:3430814/1‑71 (MQ=255)
               gCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATAACGTCGTTTCCACGTACCACCGCTCCAATAATCg                                                         <  1:2184011/71‑1 (MQ=255)
                ccTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGt                                                        <  1:3503663/71‑1 (MQ=255)
                ccTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGt                                                        <  1:3916971/71‑1 (MQ=255)
                 cTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGtt                                                       <  1:3000130/71‑1 (MQ=255)
                 cTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGtt                                                       <  1:3000129/71‑1 (MQ=255)
                    gCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTc                                                      <  1:3659400/69‑1 (MQ=255)
                     cTCAATTCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTg                                                   <  1:3398333/71‑1 (MQ=255)
                     cTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTg                                                   <  1:1212550/71‑1 (MQ=255)
                       cAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTCCAATAATCGTTCCTGGc                                                 >  1:2961885/1‑71 (MQ=255)
                       cAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTc                                                      >  1:1426014/1‑66 (MQ=255)
                        aaTGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCg                                                >  1:3117864/1‑71 (MQ=255)
                        aaTGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCg                                                >  1:2849705/1‑71 (MQ=255)
                        aaTGCGCAGATTGTCATTGGCAATCATCACGTCGTTTACACGTACCACCGCTc                                                                  >  1:1600340/1‑53 (MQ=255)
                         aTGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTCCAATAATCGTTCCTGGCgg                                               <  1:3759239/71‑1 (MQ=255)
                            cgcAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAg                                            <  1:1709793/71‑1 (MQ=255)
                             gcAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGc                                           <  1:2221390/71‑1 (MQ=255)
                              atgATTGCCAATGACAATCATCACGTCGTTTCCACGTACCACCGCTCCaat                                                              <  1:2948345/49‑1 (MQ=25)
                               aGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCtt                                         >  1:2003543/1‑71 (MQ=255)
                               aGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCtt                                         >  1:2853621/1‑71 (MQ=255)
                                 aTTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGa                                       >  1:1189050/1‑71 (MQ=255)
                                 aTTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGa                                       >  1:859621/1‑71 (MQ=255)
                                  ttGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGAt                                      >  1:2958792/1‑71 (MQ=255)
                                  ttGTCATTGGCAATCATCACGTCGTTTCCACGTACAACCGCTCCAATAATCGTTCCTGGCGGTa                                             >  1:3012056/1‑64 (MQ=255)
                                  ttGTCATTGGCAATCATCAAGTCGTTTCCAGGTACCGCCGCTCCAATAATCGTTACTGGc                                                 >  1:890594/1‑60 (MQ=255)
                                   tGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGAtt                                     >  1:1826937/1‑71 (MQ=255)
                                    gTCATTGGCAATCATCACGTCGTTTCCACGTaccac                                                                       >  1:1360584/1‑36 (MQ=255)
                                    gTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGAtt                                     <  1:2406667/70‑1 (MQ=255)
                                    gTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGAtt                                     <  1:3198318/70‑1 (MQ=255)
                                     tCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTc                                   <  1:641584/71‑1 (MQ=255)
                                      cATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCg                                  >  1:753087/1‑71 (MQ=255)
                                       aTTGGCAATCATCACGTCGTTTCCACGTACCGCCGCTCCAATAATCGtt                                                       >  1:3712220/1‑49 (MQ=255)
                                       aTTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCg                                  <  1:892676/70‑1 (MQ=255)
                                       aTTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGt                                 >  1:2983317/1‑71 (MQ=255)
                                          ggCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCa                               <  1:3650378/70‑1 (MQ=255)
                                           gCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAAt                             <  1:2507579/71‑1 (MQ=255)
                                              atcatcACGTCGTTTCCACGTACCACCGCTCCAATAATCTTTCCTGGCGGTAGCTTGATTTCGTCAATGAc                          >  1:2761396/1‑71 (MQ=255)
                                               tcatcaCGTCGTTTCCACGTACCGCCGCTCCAATAATCGTTc                                                      >  1:373914/1‑42 (MQ=255)
                                                 atcaCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGActc                        <  1:2107216/70‑1 (MQ=255)
                                                  tcaCGTCGTTTCCACGTACCGCCGCTCCAATAATCGt                                                        <  1:570922/37‑1 (MQ=255)
                                                  tcaCGTCGTTTCCACGTACCGCCGCTCCAATAATCGt                                                        <  1:1551679/37‑1 (MQ=255)
                                                  tcaCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTg                      >  1:3796793/1‑71 (MQ=255)
                                                  tcaCGGCGGTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTg                      <  1:2570243/71‑1 (MQ=255)
                                                   caCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCt                                                    <  1:2032883/40‑1 (MQ=255)
                                                   caCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGc                     >  1:3702069/1‑71 (MQ=255)
                                                    aCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCaa                              >  1:1296338/1‑61 (MQ=255)
                                                    aCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGcc                    >  1:154556/1‑71 (MQ=255)
                                                    aCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGcc                    <  1:250831/71‑1 (MQ=255)
                                                    aCGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGcc                    <  1:3370326/71‑1 (MQ=255)
                                                     cgtcgtTTCCACGTACCGCCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGcc                    >  1:2126391/1‑70 (MQ=255)
                                                      gtcgtTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGa                  >  1:2901862/1‑71 (MQ=255)
                                                      gtcgtTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGa                  >  1:3834140/1‑71 (MQ=255)
                                                       tcgtTTCCACGTACCACCGCTCCAATAATCGTTCCTg                                                   >  1:822684/1‑37 (MQ=255)
                                                       tcgtTTCCACGTACCACCGCTCCAATAATCGTTCCTg                                                   >  1:2420591/1‑37 (MQ=255)
                                                       tcgtTTCCACGTACCACCGCTCCAATAATCGTTCCTg                                                   >  1:3364987/1‑37 (MQ=255)
                                                            tCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAAcgc            >  1:2617102/1‑71 (MQ=255)
                                                             ccACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAAcgcg           >  1:193923/1‑71 (MQ=255)
                                                             ccACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAAcgcg           >  1:1988682/1‑71 (MQ=255)
                                                              cACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGt          <  1:616661/71‑1 (MQ=255)
                                                              cACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGt          <  1:602488/71‑1 (MQ=255)
                                                                 gTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAg       >  1:3796840/1‑71 (MQ=255)
                                                                 gTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAg       >  1:3338160/1‑71 (MQ=255)
                                                                  tACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAAc              <  1:2485478/63‑1 (MQ=255)
                                                                   accaccGCTCCAATAATCGTTCCTGGCGGTAGCTTGAtt                                     <  1:2166471/39‑1 (MQ=255)
                                                                   accaccGCTCCAATAATCGTTCCTGGCGGTAGCTTGAtt                                     <  1:3368389/39‑1 (MQ=255)
                                                                    ccaccGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTg    <  1:2532657/71‑1 (MQ=255)
                                                                    ccaccGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTg    >  1:3448273/1‑71 (MQ=255)
                                                                     caccGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTGc   >  1:3439283/1‑71 (MQ=255)
                                                                      accGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGa                           <  1:809715/46‑1 (MQ=255)
                                                                      accGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTGCt  <  1:1307833/71‑1 (MQ=255)
                                                                      accGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTGCt  >  1:1972947/1‑71 (MQ=255)
                                                                      |                                                                      
CATAATTACGTGATCGCCTTGCTCAATGCGCAGATTGTCATTGGCAATCATCACGTCGTTTCCACGTACCACCGCTCCAATAATCGTTCCTGGCGGTAGCTTGATTTCGTCAATGACTCTGCCGACAACGCGTGAGGTGCT  >  minE/2754200‑2754340

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: