Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,135 G→C 14.9% G824G (GGC→GGG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,1350G→C14.9% 23.0 ‑4.0 ‑0.074G824G (GGC→GGG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (10/1):  ref base (31/32):  total (41/33)
Fisher's exact test strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.58e-01

AGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTCA  >  NC_014034/681052‑681235
                                                                                   |                                                                                                    
aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGc                                                                                     >  1:3916708/1‑101 (MQ=25)
aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGc                                                                                     >  1:1506451/1‑101 (MQ=25)
   tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                                  >  1:3539066/1‑98 (MQ=25)
   tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                                  >  1:2103216/1‑98 (MQ=25)
    gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                                                 >  2:3074367/1‑97 (MQ=25)
     cgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccac                                                                                >  1:2113462/1‑96 (MQ=25)
     cgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccaa                                                                                >  2:2617187/1‑96 (MQ=25)
      cccttGCTCAGCGCCGCCACCTGCTGCGAGGTTATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                                               <  1:1625932/97‑7 (MQ=21)
       ccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGctccactg                                                                              >  1:3933321/1‑94 (MQ=25)
       ccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactg                                                                              >  1:3101068/1‑94 (MQ=25)
         gTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggc                                                                            >  1:2777898/1‑92 (MQ=21)
                                                   tgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgt                                  <  1:503909/92‑1 (MQ=38)
                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc                                 <  2:1253766/93‑1 (MQ=255)
                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc                                 <  1:520833/93‑1 (MQ=255)
                                                     cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg                                >  1:2463190/8‑101 (MQ=255)
                                                      gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt                               >  1:2199396/7‑101 (MQ=255)
                                                       ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg                              <  1:2791572/96‑1 (MQ=255)
                                                       ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg                              <  2:1750863/96‑1 (MQ=255)
                                                        cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc                             >  2:570188/5‑101 (MQ=255)
                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt                            >  2:4012875/4‑101 (MQ=255)
                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt                            <  1:2684716/98‑1 (MQ=255)
                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGCCGTGCt                            >  1:3482095/4‑101 (MQ=255)
                                                         gtcgTCAGCGCCACGATGTCAGCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt                            >  1:3118963/4‑101 (MQ=255)
                                                          tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc                           <  1:3109555/99‑1 (MQ=255)
                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa                          <  2:3768971/100‑1 (MQ=255)
                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa                          <  1:1727738/100‑1 (MQ=255)
                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa                          <  2:3173425/100‑1 (MQ=255)
                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg                         >  2:386776/1‑101 (MQ=255)
                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg                         <  2:2268405/101‑1 (MQ=255)
                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg                         >  2:724654/1‑101 (MQ=255)
                                                             tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc                        <  1:188795/101‑1 (MQ=255)
                                                             gcAGCGCCACGATGTCATCGGTGCCGAGCCCCTTGGTCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCGGCGTCGTGCTCAgc                        <  2:2476925/100‑1 (MQ=38)
                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg                       >  1:2965703/1‑101 (MQ=255)
                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc                      >  2:3045110/1‑101 (MQ=255)
                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc                      <  1:116740/101‑1 (MQ=255)
                                                                gcgcGACGGTGTCATGGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGcgcg                     <  2:1124803/101‑1 (MQ=255)
                                                                 cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCgg                    <  2:1116459/101‑1 (MQ=255)
                                                                  gccacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGc                   >  1:4196210/4‑101 (MQ=255)
                                                                    cacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCa                 >  1:3998391/2‑101 (MQ=255)
                                                                    cacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCa                 <  2:2006800/100‑1 (MQ=255)
                                                                    cacgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCa                 >  2:621098/2‑101 (MQ=255)
                                                                     acgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCAt                >  2:3379531/1‑101 (MQ=255)
                                                                     acgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCAt                >  1:1402418/1‑101 (MQ=255)
                                                                      cgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATc               >  1:3538536/1‑101 (MQ=255)
                                                                      cgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATc               >  1:1690078/1‑101 (MQ=255)
                                                                      cgaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATc               <  1:2779193/101‑1 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              <  2:1090733/101‑1 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              >  1:1719755/1‑101 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              >  2:3726913/1‑101 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              >  2:49605/1‑101 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              <  2:2965481/101‑1 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              <  1:4028951/101‑1 (MQ=255)
                                                                       gaTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCt              >  1:3675456/1‑101 (MQ=255)
                                                                        aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg             <  2:2829093/101‑1 (MQ=255)
                                                                        aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg             <  2:3341033/101‑1 (MQ=255)
                                                                        aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg             >  2:3221008/1‑101 (MQ=255)
                                                                        aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg             <  1:48190/101‑1 (MQ=255)
                                                                        aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg             <  1:3953337/101‑1 (MQ=255)
                                                                        aTGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTg             <  2:677734/101‑1 (MQ=255)
                                                                         tGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgc            >  2:1060578/1‑101 (MQ=255)
                                                                         tGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgc            >  2:2621255/1‑101 (MQ=255)
                                                                         tGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgc            >  2:2774708/1‑101 (MQ=255)
                                                                         tGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgc            >  1:3131693/1‑101 (MQ=255)
                                                                         tGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACCCTGGCCGTATTCAGCGCCGCAAGCCGCGTCGCGCTCAGCGCGGcccgctcc            >  1:2240894/1‑95 (MQ=255)
                                                                          gTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcg           >  1:4051541/1‑101 (MQ=255)
                                                                           tCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcgc          <  2:931418/101‑1 (MQ=255)
                                                                           tCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcgc          <  1:74513/101‑1 (MQ=255)
                                                                           tCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTgcgc          <  1:2715068/101‑1 (MQ=255)
                                                                            cATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGcc         <  1:2788530/101‑1 (MQ=255)
                                                                               cGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAg      <  1:1618409/101‑1 (MQ=255)
                                                                                ggTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAgg     >  2:3039916/1‑101 (MQ=255)
                                                                                 gTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGt    >  2:2406040/1‑101 (MQ=255)
                                                                                 gTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGt    >  1:1052388/1‑101 (MQ=255)
                                                                                  tGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTc   <  1:2921293/101‑1 (MQ=255)
                                                                                  tGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTc   <  1:2774927/101‑1 (MQ=255)
                                                                                   gCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTCa  <  1:4077471/101‑1 (MQ=255)
                                                                                   gCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTCa  <  1:3727167/101‑1 (MQ=255)
                                                                                   gCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTCa  >  1:3204535/1‑101 (MQ=255)
                                                                                   |                                                                                                    
AGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGGCCATCTGCGCCGAGGTCA  >  NC_014034/681052‑681235

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG