Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,882,516 +A 5.9% coding (291/381 nt) atpI ← F0F1 ATP synthase subunit I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,882,5161.→A 5.9% 2.3 ‑2.5 ‑0.034coding (291/381 nt)atpIF0F1 ATP synthase subunit I
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (24/8):  total (24/10)
Fisher's exact test strand distribution p-value = 8.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.43e-01

CCCACGTAACGATCAGCGGCAAGAATACCGCCTTT.AAAACCGCCAACGCCACCACCAGTAACACCAACA  >  REL606/3882482‑3882550
                                   |                                  
cccACGTAACGATCAGCGGCAAGAATACCGCCTTT.a                                   >  1:2760742/1‑36 (MQ=255)
 ccACGTAACGATCAGCGGAAAAAATACCGCCTTT.aa                                  <  1:7747773/36‑1 (MQ=25)
 ccACGGAACGATCAGCGGCAAGAATACCGCTTTT.aa                                  >  1:7626149/1‑36 (MQ=37)
  cACGTAACGATCAGCGGCAAGAATACCGCCTTT.aaa                                 >  1:732151/1‑36 (MQ=255)
  cACGAAACGATCAGCGGCAAGAATACCGCCTTT.aaa                                 <  1:6044408/36‑1 (MQ=255)
  cACATAACGATCAGCGGAAAAAATACCGCCTTT.aaa                                 <  1:4704531/36‑1 (MQ=21)
   aCGTAACGATCAGCGGCAAGAATACCGCCTTT.aaaa                                >  1:905778/1‑36 (MQ=255)
   aCGTAACGATCAGCGGCAAGAATACCGCCTTT.aaaa                                >  1:7807538/1‑36 (MQ=255)
   aCGTAACGATCAGCGGAAAGAATACCGCCTTT.aaaa                                <  1:4041259/36‑1 (MQ=255)
     gTAACGATCAGCGGCAAGAATACCGCCTTT.AAAAcc                              >  1:3519421/1‑36 (MQ=255)
      tAACGATCAGCGGCAAGAATACCGCCTTT.TAAACCg                             >  1:410126/1‑36 (MQ=255)
      tAACGATCAGCGGCAAGAATACCGCCTTT.AAAACCg                             >  1:5936630/1‑36 (MQ=255)
        aCGATCAGCGGCAAGAATACCGCCTTT.AAAACCgac                           >  1:6055524/1‑34 (MQ=255)
        aCGATCAGCGGCAAGAATACCGCATTT.AAAACCGcc                           >  1:4745902/1‑36 (MQ=255)
          gATCAGCGGCAAGAATACCGCCTTT.AAAACCGCCaa                         >  1:679569/1‑36 (MQ=255)
          gATCAGCGGCAAGAATACCGCCTTT.AAAACCGCCaa                         >  1:7354081/1‑36 (MQ=255)
           aTCAGCGGCAAGAATACCGCCTTT.AAAACCGCCAAc                        >  1:1430161/1‑36 (MQ=255)
             cAGCGGCAAGAATACCGCCTTT.AAAACCGCCAACGc                      >  1:6601641/1‑36 (MQ=255)
             aaGCGGAAAAAAAACCGCCTTT.AAAACCGCCAACGc                      <  1:5934848/35‑1 (MQ=21)
              aGCGGCAAGAATACCGCCTTT.AAAACCGCCAACGcc                     >  1:4402213/1‑36 (MQ=255)
               gCGGCAAGAATACCGCCTTT.AAAACCGCCAACGcca                    >  1:7291016/1‑36 (MQ=255)
                cGGCAAGAATACCGCCTTT.AAAACCGCCAACGccac                   >  1:1345236/1‑36 (MQ=255)
                cGGCAAGAATACCGCCTTT.AAAACCGCCAACGccac                   >  1:7698994/1‑36 (MQ=255)
                   cAAGAATACCGCCTTT.AAAACCGCCAACGccaccac                >  1:2188488/1‑36 (MQ=255)
                    acgAATACACCCTTTAAAAACCGCCAACGccaccac                <  1:1464115/34‑1 (MQ=16)
                     aGAATACCGCCTTT.AAAACCGCCAACGccaccacca              >  1:251724/1‑36 (MQ=255)
                      gAATACCGCCTTT.AAAACCGCCAACGCCACCACCAg             >  1:5504635/1‑36 (MQ=255)
                      gAATACCGCCTTT.AAAACCGCCAAAGCCACCACCAg             >  1:688496/1‑36 (MQ=255)
                       aaTACCGCCTTT.AAAACCGCCAACGCCACCACCAGt            >  1:2243185/1‑36 (MQ=255)
                        aTACCGCCTTT.AAAACCGCCAACGCCACCACCAGTa           >  1:1822643/1‑36 (MQ=255)
                        aTACCGCCTTT.AAAACCGCCAACGCCACCACCAGTa           >  1:7247880/1‑36 (MQ=255)
                         aaCCCCCTTA.AAAACCGCCAACGCCACCACCAGTaa          <  1:4879329/35‑1 (MQ=25)
                          aCCGCCTTT.AAAACCGCCAACGCCACCACCAGTAac         >  1:7169075/1‑36 (MQ=255)
                          aCCGACTTT.AAAACCGCCAACGCCACCACCAGTAac         >  1:7486878/1‑36 (MQ=255)
                           cccccTTA.AAAACCGCCAACGCCACCACCAGTAaca        <  1:6987352/33‑1 (MQ=37)
                            ccccTTTAAAAACCGCCAACGCCACCACCAGTAaca        <  1:3795606/34‑1 (MQ=25)
                            ccccTTA.AAAACCGCCAACGCCACCACCAGTAacac       <  1:87935/34‑1 (MQ=37)
                            cGCCTTT.AAAACCGCCAACGCCACCACCAGTAacac       >  1:2446645/1‑36 (MQ=255)
                            cGCCTTA.AAAACCGCCAACGCCACCACCAGTAacac       <  1:4764601/36‑1 (MQ=255)
                             gCCTTT.AAAACCGCCAACGCCACCACCAGTAACAcc      >  1:6242087/1‑36 (MQ=255)
                             cccTTA.AAAACCGCCAACGCCACCACCAGTAACAcc      <  1:243683/35‑1 (MQ=255)
                              ccTTA.AAAACCGCCAACGCCACCACCAGTAACACCa     <  1:1711365/36‑1 (MQ=255)
                               cTTA.AAAACCGCAAACGCCACCACCAGTAACACCaa    <  1:6019910/36‑1 (MQ=25)
                                tat.AAAACCGCCAACTCCACCACCAGTAACACCAAc   <  1:5267941/34‑1 (MQ=38)
                                 tt.AAAACCGCCAACGCCACCACCAGTAACACCAACa  >  1:797834/1‑36 (MQ=255)
                                   |                                  
CCCACGTAACGATCAGCGGCAAGAATACCGCCTTT.AAAACCGCCAACGCCACCACCAGTAACACCAACA  >  REL606/3882482‑3882550

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG