Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,867,898 +G 4.4% coding (427/1686 nt) fadD ← acyl‑CoA synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,867,8981.→G 4.4% 3.7 ‑2.7 ‑0.645coding (427/1686 nt)fadDacyl‑CoA synthase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (10/33):  total (12/33)
Fisher's exact test strand distribution p-value = 6.67e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.75e-01

ATCGCCCATACGGGTCAGAATTACGTGCTGAACGGC.GGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAA  >  REL606/1867863‑1867933
                                    |                                   
aTCGCCCATACGGGTCAGAATTACGTGCTGAAcggc                                      >  1:3622443/1‑36 (MQ=255)
aTCGCCCATACGGGTCAGAATTACGTGCTGAAcggc                                      <  1:5158344/36‑1 (MQ=255)
 tCGCCCATACGGGTCAGAATTACGTGCTGAAcggc.g                                    <  1:1807138/36‑1 (MQ=255)
   gCCCATACGGGTCAGAATTACGTGCTGAACGGC.GGt                                  <  1:6209227/36‑1 (MQ=255)
   gCCCATACGGGTCAGAATTACGTGCTGAACGGC.GGt                                  <  1:6643137/36‑1 (MQ=255)
    cccATACGGGTCAGAATTACGTGCTGAAcggc.gggg                                 >  1:866572/1‑34 (MQ=255)
    ccaaTACGGGTCAGAATTACGTGCTGAACGGC.GGtt                                 <  1:6459216/33‑1 (MQ=255)
     ccATACGGGTCAGAATTACGTGCTGAACGGC.GGttt                                <  1:5276575/36‑1 (MQ=255)
     ccATACGGGTCAGAATTACGTGCTGAACGGC.GGttt                                <  1:2374125/36‑1 (MQ=255)
     ccATACGGGTCAGAATTACGTGCTGAACGGC.GGgtt                                >  1:2133741/1‑36 (MQ=255)
      cATACGGGTCAGAATTACGTGCTGAACGGC.GGtttt                               <  1:3636955/36‑1 (MQ=255)
      cATACGGGTCAGAATTACGTGCTGAACGGC.GGtttt                               <  1:654235/36‑1 (MQ=255)
       aaaCGGGTCAGCATGACGTGCTGAACGGC.TGttttt                              >  1:4767240/3‑35 (MQ=16)
       aTACGGGTCAGAATTACGTGCTGAACGGC.GGgtttt                              >  1:6769958/1‑35 (MQ=255)
       aTACGGGTCAGAATTACGTGCTGAACGGC.GGgtttt                              >  1:2223666/1‑35 (MQ=255)
        tACGGGTCAGAATTACGTGCTGAACGGC.GGtttttt                             >  1:3757393/1‑34 (MQ=255)
        tACGGGTCAGAATTACGTGCTGAACGGC.GGgttttt                             >  1:4154494/1‑34 (MQ=37)
        tACGGGTCAGAATTACGTGCTGAACGGC.GGTTTTAt                             <  1:4123256/36‑1 (MQ=255)
         tcGGGTCAGAATTACGTGCTGAACGGC.GGTTTTATc                            <  1:912304/35‑1 (MQ=255)
         aCGGGTTAGAATTACGTGCTGAACGGC.GGTTTTATc                            <  1:7366256/36‑1 (MQ=255)
         aCGGGTCAGAATTACGTGCTGAACGGC.GGttttttt                            >  1:146122/1‑33 (MQ=255)
         aCGGGTCAGAATTACGTGCTGAACGGC.GGTTTTATc                            <  1:976156/36‑1 (MQ=255)
         aCGGGTCAGAATTACGTGCTGAACGGC.GGTTTTATc                            <  1:3569959/36‑1 (MQ=255)
          cGGGTCAGAATTACGTGCTGAACGGC.GGTTTTATCa                           <  1:935790/36‑1 (MQ=255)
          cGGGTCAGAATAACGTGCTGAACGGC.GGTTTTATCa                           <  1:5750463/36‑1 (MQ=255)
           tggTCATACTTACATGCTGAGCGTC.GGTTTTATCaa                          <  1:2739250/35‑1 (MQ=12)
            ggTCAGAATTACGTGCTGAACGGC.GGTTTTATCAAc                         <  1:607744/36‑1 (MQ=255)
             gTCAGAATTACGTGCTGAACGGC.GGTTTTATCAAcc                        <  1:4138087/36‑1 (MQ=255)
              tCAGAATTACGTGCTGAACGGC.GGTTTTATCAACCa                       <  1:1716017/36‑1 (MQ=255)
               cAGAATTACTTGCTGAACGGC.GGTTTTATCAACCAc                      <  1:5569835/36‑1 (MQ=255)
               cAGAATTACGTGCTGAACGGC.GGTTTTATCAACCAc                      <  1:2561202/36‑1 (MQ=255)
                 gAATTACGTGCTGAACGGC.GGTTTTATCAACCACtt                    <  1:6853369/36‑1 (MQ=255)
                    ttACGTGCTGAACGGC.GGGTTTTTTTAACACttttt                 >  1:3028206/1‑36 (MQ=12)
                     tACGTGCTGAACGGC.GGTTTTATCAACCACTTTTTc                <  1:93594/36‑1 (MQ=255)
                     tACGTGCTGAACGGC.GGGTTTTTTAAACACtttttt                >  1:7240920/1‑35 (MQ=16)
                      aCGTGCTGAACGTC.GGTTTTATCAACCACTTTTTcc               <  1:2718793/36‑1 (MQ=255)
                      aCGTGCTGAACGGCGGGTTTTTTTATCCCCtttttt                >  1:5037541/1‑35 (MQ=0)
                      aCGTGCTGAACGGCGGGTTTTTTTAACACCTTTTTc                >  1:4955492/1‑36 (MQ=2)
                      aCGTGCTGAACGGC.GGTTTTATCAACCACTTTTTcc               <  1:383209/36‑1 (MQ=255)
                       cGTGCTGAACGGC.GGTTTTATCAACCACTTTTTCCa              <  1:5177821/36‑1 (MQ=255)
                       cGTGCTGAACGGC.GGTTTTATCAACCACTTTTTCCa              <  1:1389430/36‑1 (MQ=255)
                         tGCTGAACGGC.GGTTTTATCAACCACTTTTTCCAgt            <  1:5529897/36‑1 (MQ=255)
                         tGCTGAACGGC.GGTTTTATCAACCACTTTTTCCAgt            <  1:6609253/36‑1 (MQ=255)
                           cTGAACGGC.GGTTTTATCACCCACTTTTTCCAgtgt          <  1:6200708/36‑1 (MQ=255)
                           cTGAACGGC.GGTTTTATCAACCACTTTTTCCAgtgt          <  1:7889751/36‑1 (MQ=255)
                             gAACGGC.GGTTTTTTTAACCACCTTTTTCCctgtgt        >  1:165626/1‑36 (MQ=2)
                              aaCGGC.GGTTTTATCAACCACTTTTTCCAGtgtgtg       <  1:2285020/36‑1 (MQ=255)
                               aCGGC.GGTTTTATCAACCACTTTTTCCGGTGTGTGa      <  1:6982619/36‑1 (MQ=255)
                               aCGGC.GGTTTTATCAACCACTTTTTCCAGTGTGTGa      <  1:59562/36‑1 (MQ=255)
                                 ggc.ggTTTTATCAACCACTTTTTCCAGTGTGTGAGc    <  1:2719501/36‑1 (MQ=255)
                                  gc.ggTTTTATCAACCACTTTTTCCAGTGTGTGAGCa   >  1:7115578/1‑36 (MQ=255)
                                   c.ggTTTTATCAACCACTTTTTCCAGTGTGTGAGCaa  >  1:5772280/1‑36 (MQ=255)
                                   c.ggTTTTATCAACCACTTTTTCCAGTGTGTGAGCaa  <  1:1325593/36‑1 (MQ=255)
                                    |                                   
ATCGCCCATACGGGTCAGAATTACGTGCTGAACGGC.GGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAA  >  REL606/1867863‑1867933

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG