Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,606,814 G→A 3.7% A86A (GCG→GCA deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,606,8140G→A 3.7% 2.7 ‑0.4 ‑0.1134A86A (GCG→GCAdeoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (73/56):  total (75/59)
Fisher's exact test strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.62e-01

Only 100 of 158 total aligned reads displayed.

ACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTG  >  REL606/4606780‑4606848
                                  |                                  
aCTTCCCACACGGTAACGATGACATCGAAATcgcg                                    <  1:3880321/35‑1 (MQ=255)
aCTTCCCACACGGTAACGATGACATCGAAATcgcg                                    <  1:4835492/35‑1 (MQ=255)
 cTTCACACACGGTAACGATGACATCGAAATCgcgc                                   <  1:2796320/35‑1 (MQ=255)
  ttCCCACACGGTAACGATGACATCGAAATCGCGCt                                  >  1:17642076/1‑35 (MQ=255)
  ttCCCACACGGTAACGATGACATCGAAATCGCACt                                  <  1:16366218/35‑1 (MQ=255)
   tCCCACACGGTAACGATGACATCGAAATCGCGCTg                                 >  1:19881750/1‑35 (MQ=255)
   tCCCACACGGTAACGATGACATCGAAATCGCGCTg                                 >  1:3850645/1‑35 (MQ=255)
   tCCCACACGGTAACGATGACATCGAAATCGCGCTg                                 <  1:15466028/35‑1 (MQ=255)
    cccACACTGTAACGATGACATCGAAATCGCGCTgg                                <  1:16106837/35‑1 (MQ=255)
    cccACACGGTAACGATGACATCGAAATCGCGCTgg                                >  1:10212501/1‑35 (MQ=255)
    cccACACGGTAACGATGACATCGAAATCGCGCTgg                                >  1:8935200/1‑35 (MQ=255)
    cccACACGGTAACGATGACATCGAAATCGCGCTgg                                <  1:13907151/35‑1 (MQ=255)
    cccACACGGTAACGATGACATCGAAATCGCGCTgg                                <  1:621432/35‑1 (MQ=255)
    cccACACGGTAACGATGACATCGAAATCGCGCCgg                                >  1:1706282/1‑35 (MQ=255)
     cgacacGGTAACGATGACATCGAAATCGCGCTGGc                               <  1:5528533/33‑1 (MQ=255)
     ccACACGGTAACGATGACATCGAAATCGCGCTGGc                               >  1:15400068/1‑35 (MQ=255)
      cacgCGGTAACGATGACATCGAAATCGCGCTGGCa                              <  1:14267794/35‑1 (MQ=255)
      cacaCGGTAACGCTGACATCGAAATCGCGCTGGCa                              >  1:19942245/1‑35 (MQ=255)
      cacaCGGTAACGATGACATCGAAATCGCGCTGGCa                              <  1:13293619/35‑1 (MQ=255)
       acacGGTAACGATGACATCGAAATCGCGCTGGCAg                             >  1:1353848/1‑35 (MQ=255)
       acacGGTAACGATGACATCGAAATCGCGCTGGCAg                             >  1:4607626/1‑35 (MQ=255)
        cacGGTAACGATGACATCGAAATCGCGCTGGCAGa                            <  1:7878517/35‑1 (MQ=255)
        cacGGTAACGATGACATCGAAATCGCGCTGGCAGa                            <  1:8868122/35‑1 (MQ=255)
         acGGTAACGATGACATCGAAATCGCGCTGGCAGaa                           >  1:6564452/1‑35 (MQ=255)
          cGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                          <  1:16890826/35‑1 (MQ=255)
          cGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                          >  1:2067101/1‑35 (MQ=255)
          cGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                          >  1:6600411/1‑35 (MQ=255)
          cGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                          >  1:9143115/1‑35 (MQ=255)
          cGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                          >  1:10072549/1‑35 (MQ=255)
          cGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                          <  1:15053291/35‑1 (MQ=255)
           ngTAACGATGACATCGAAATCGCGCTGGCAGAAAc                         >  1:15952800/2‑35 (MQ=255)
           ggTAACGATGACATCGAAATCGCGCTGGCAGAAAc                         <  1:12415520/35‑1 (MQ=255)
            gTAACGATGACATCGAAATCGCGCTGGCAGAAAcc                        >  1:3995646/1‑35 (MQ=255)
             tAACGATGACATCGAAATCGCGCTGGCAGAAAccc                       >  1:14050382/1‑35 (MQ=255)
              aaCGATGACATCGAAATCGCGCTGGCAGAAACCCg                      <  1:19420494/35‑1 (MQ=255)
              aaCGATGACATCGAAATCGCGCTGGCAGAAACCCg                      <  1:17560166/35‑1 (MQ=255)
              aaCGATGACATCGAAATCGCGCTGGCAGAAACCCg                      >  1:6258311/1‑35 (MQ=255)
              aaCGATGACATCGAAATCGCGCTGGCAGAAACCCg                      >  1:6352990/1‑35 (MQ=255)
               aCGATGACATCGAAATCGCGCTGGCAGAAACCCGt                     <  1:12054596/35‑1 (MQ=255)
               aCGATGACATCGAAATCGCGCTGGCAGAAACCCGt                     >  1:11477586/1‑35 (MQ=255)
               aCGATGACATCGAAATCGCGCTGGCAGAAACCCGt                     >  1:12249940/1‑35 (MQ=255)
               aCGATGACATCGAAATCGCGCTGGCAGAAACCCGt                     >  1:10540212/1‑35 (MQ=255)
               aCGATGACATCGAAATCGCGCTGGCAGAAACCCGt                     <  1:2811125/35‑1 (MQ=255)
                cGATGACATCGAAATCGCGCTGGCAGAACCCCGTg                    <  1:9560211/35‑1 (MQ=255)
                cGATGACATCGAAATCGCGCTGGCAGAAACCCGTg                    <  1:14328853/35‑1 (MQ=255)
                cGATGACATCGAAATCGCGCTGGCAGAAACCCGTg                    <  1:20611069/35‑1 (MQ=255)
                cGATGACATCGAAATAGCGCTGGCAGAAACCCGTg                    <  1:9135962/35‑1 (MQ=255)
                 gATGACATCGAAATCGCGCTGGCAGAAACCCGTGc                   >  1:7349709/1‑35 (MQ=255)
                 gATGACATCGAAATCGCGCTGGCAGAAACCCGTGc                   <  1:9755501/35‑1 (MQ=255)
                 gATGACATCGAAATCGCGCTGGCAGAAACCCGTGc                   >  1:15562459/1‑35 (MQ=255)
                  aTGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                  <  1:16581706/35‑1 (MQ=255)
                  aTGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                  >  1:16883901/1‑35 (MQ=255)
                  aTGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                  >  1:11491253/1‑35 (MQ=255)
                  aTGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                  >  1:4549715/1‑35 (MQ=255)
                  aTGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                  <  1:341919/35‑1 (MQ=255)
                   tGACATCGAAATCGCGCTGGCAGAAACCCGTGCgg                 <  1:5899533/35‑1 (MQ=255)
                   tGACATCGAAATCGCGCTGGCAGAAACCCGTGCgg                 <  1:5759694/35‑1 (MQ=255)
                   tGACATCGAAATCGCGCTGGCAGAAACCCGTGCgg                 <  1:17558572/35‑1 (MQ=255)
                   tGACATCGAAATCGCACTGGCAGAAACCCGTGCgg                 <  1:927692/35‑1 (MQ=255)
                    gACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                >  1:8549546/1‑35 (MQ=255)
                     aCATCGAAATCGCGCTGGCAGAAACCCGTGCGGCa               >  1:20220367/1‑35 (MQ=255)
                     aCATCGAAATCGCGCTGGCAGAAACCCGTGCGGCa               <  1:7947605/35‑1 (MQ=255)
                      cATCGAAATCGCGCTGGCAGAAACCCGTGCGGCaa              <  1:6095898/35‑1 (MQ=255)
                      cATCGAAATCGCGCTGGCAGAAACCCGTGCGGCaa              <  1:19651388/35‑1 (MQ=255)
                      cATCGAAATCGCGCTGGCAGAAACCCGTGCGGCaa              >  1:14173708/1‑35 (MQ=255)
                      cATCGAAATCGCGCTGGCAGAAACCCGTGCGGCaa              >  1:2485696/1‑35 (MQ=255)
                       aTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAt             <  1:17359452/35‑1 (MQ=255)
                       aTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAt             <  1:8061986/35‑1 (MQ=255)
                       aTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAt             <  1:8388088/35‑1 (MQ=255)
                       aTCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAt             <  1:15657374/35‑1 (MQ=255)
                       aTCGAAATCGCACTGGCAGAAACCCGTGCGGCAAt             >  1:5978048/1‑35 (MQ=255)
                        tCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc            >  1:5169901/1‑35 (MQ=255)
                        tCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc            >  1:5298147/1‑35 (MQ=255)
                        tCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc            <  1:6743457/35‑1 (MQ=255)
                        tCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc            <  1:16748245/35‑1 (MQ=255)
                        gcGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc            <  1:5412673/34‑1 (MQ=255)
                         cGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCg           >  1:6841613/1‑35 (MQ=255)
                         cGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCg           >  1:19754497/1‑35 (MQ=255)
                         cGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCg           >  1:18229060/1‑35 (MQ=255)
                          gAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGc          <  1:1934197/35‑1 (MQ=255)
                          gAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGc          >  1:3582863/1‑35 (MQ=255)
                          gAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGc          >  1:13509282/1‑35 (MQ=255)
                          gAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGc          >  1:2271345/1‑35 (MQ=255)
                          gAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGc          >  1:17435382/1‑35 (MQ=255)
                          gAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGc          >  1:11618036/1‑35 (MQ=255)
                           aaaTCGCGCTGGCAGAAACCCGTGCGGCAATCGcc         >  1:1291359/1‑35 (MQ=255)
                           aaaTCGCGCTGGCAGAAACCCGTGCGGCAATCGcc         >  1:15118582/1‑35 (MQ=255)
                            aaTCGCGCTGGCAGAAACCCGTGCGGCAATCGCCt        >  1:8867813/1‑35 (MQ=255)
                             aTCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa       >  1:4252065/1‑35 (MQ=255)
                             aTCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTa       <  1:8407977/35‑1 (MQ=255)
                               cgcgCTGGCAGAAACCCGTGCGGCAATCGCCTACg     >  1:14999565/1‑35 (MQ=255)
                               cgcgCTGGCAGAAACCCGTGCGGCAATCGCCTACg     <  1:1189274/35‑1 (MQ=255)
                                gcgcTGGCAGAAACCCGTGCGGCAATCGCCTACgg    <  1:1068563/35‑1 (MQ=255)
                                gcgcTGGCAGAAACCCGTGCGGCAATCGCCTACgg    >  1:8747124/1‑35 (MQ=255)
                                gcgcTGGCAGAAACCCGTGCGGCAATCGCCTACgg    >  1:13239599/1‑35 (MQ=255)
                                 cgcTGGCAGAAACCCGTGCGGCAATCGCCTACgg  >  1:432973/1‑34 (MQ=255)
                                 cgcTGGCAGAAACCCGTGCGGCAATCGCATACGGt   <  1:17664439/35‑1 (MQ=255)
                                  gcTGGCAGAAACCCGTGCGGCAATCGCCTACGGTg  <  1:12793523/35‑1 (MQ=255)
                                  gcTGGCAGAAACCCGTGCGGCAATCGCCTACGGTg  >  1:12638577/1‑35 (MQ=255)
                                  gcTGGCAGAAACCCGTGCGGCAATCGCCTACGGTg  >  1:20551116/1‑35 (MQ=255)
                                  |                                  
ACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTG  >  REL606/4606780‑4606848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.