Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,908,398 +G 7.1% coding (1329/2871 nt) xdhD → fused predicted xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit/Fe‑S binding subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,908,3981.→G 7.1% 8.7 ‑3.0 ‑0.342coding (1329/2871 nt)xdhDfused predicted xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit/Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (13/26):  total (16/26)
Fisher's exact test strand distribution p-value = 4.88e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.38e-01

CATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGACATCCTGCTGGAAG  >  REL606/2908363‑2908432
                                    |                                  
cATAAAGTTCATGTCATTAAAGAGCGAGTTggcggc                                     >  1:6605155/1‑36 (MQ=255)
 aTAAAGTTCATGTCATTAAAGAGCGAGTTGgcggc.g                                   >  1:3206923/1‑36 (MQ=255)
  tAAAGTTCATGTCATTAAAGAGCGAGTTGGcggc.gg                                  <  1:7949151/36‑1 (MQ=255)
  tAAAGTTCATGTCATTAAAGAGCGAGTTGGcggc.gg                                  <  1:7457968/36‑1 (MQ=255)
   aaaGTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGt                                 <  1:5778449/36‑1 (MQ=255)
   aaaGTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGt                                 <  1:2256622/36‑1 (MQ=255)
    aaGTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGtt                                <  1:6161127/36‑1 (MQ=255)
      gTTCATGTTATTAAAGAGCGAGTTGGCGGC.GGtttt                              <  1:4843691/36‑1 (MQ=255)
      gTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGtttt                              >  1:4687850/1‑36 (MQ=255)
      gTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGtttt                              >  1:6420532/1‑36 (MQ=255)
       ttCATGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTg                             <  1:6623403/36‑1 (MQ=255)
        tCATGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTgg                            <  1:5082726/36‑1 (MQ=255)
        tCATGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTgg                            <  1:407172/36‑1 (MQ=255)
         cATGTCATTAAAGAGCGAGTTGGCGGC.GGttttttt                           >  1:5916982/1‑33 (MQ=255)
         cATGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTTGt                           >  1:2069614/1‑36 (MQ=255)
           tGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTTGttt                         >  1:7309398/1‑35 (MQ=255)
           tGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTGGTTc                         <  1:2293315/36‑1 (MQ=255)
             tCATTAAAGAGCGAGTTGGCGGC.GGTTTTTGtttcc                       >  1:2494954/1‑33 (MQ=37)
             tCATTAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCa                       <  1:134118/36‑1 (MQ=255)
             tCATTAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCa                       <  1:6184096/36‑1 (MQ=255)
             tCATTAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCa                       <  1:5914718/36‑1 (MQ=255)
               aTTAAAGAGCGAGTTGGCGGCGGGTTTTGTTTCCaa                      >  1:5877128/1‑36 (MQ=25)
                ttAAAGAGCGAGTTGGCGGCGGGTTTTTTTTCCaac                     >  1:4431382/1‑35 (MQ=16)
                ttAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCAAAc                    <  1:5171150/36‑1 (MQ=255)
                 tAAAGAGCGAGTTGTCGGC.GGTTTTGGTTCCAAACa                   <  1:6422251/36‑1 (MQ=255)
                 tAAAGAGCGAGTTGGCGGCGGGTTTTGTTTCCaata                    >  1:5470028/1‑34 (MQ=21)
                 tAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCAAACa                   <  1:2497429/36‑1 (MQ=255)
                 tAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCAAACa                   <  1:8385949/36‑1 (MQ=255)
                   aaGAGCGAGTTGGCGGC.GGTTTTTGTTCCCAACAgg                 >  1:6145573/1‑36 (MQ=37)
                   aaGAGCGAGTTGGCGGC.GGGTTTTGTTTCCAACAgg                 >  1:3809638/1‑36 (MQ=19)
                    agagCGAGTTGGCGGC.GGTTTTTGTTTCCAACAggg                >  1:7603791/1‑35 (MQ=21)
                     gagCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGAc               <  1:4829763/36‑1 (MQ=255)
                     gagCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGAc               <  1:546910/36‑1 (MQ=255)
                     gagCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGAc               <  1:2600487/36‑1 (MQ=255)
                     gagCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGAc               <  1:2542104/36‑1 (MQ=255)
                     gagCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGAc               <  1:7848212/36‑1 (MQ=255)
                      agCGAGTTGGCGGC.GGTTTTTGTTTCAAACAgggga              >  1:6392214/1‑33 (MQ=25)
                       gCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGACAt             <  1:7153016/36‑1 (MQ=255)
                        cGAGTTGGCGGC.GGTTTTTGTTTCCAACAGGACATc            >  1:1225169/1‑36 (MQ=25)
                        cGAGTTGGCGGC.GGTTTTTGTTTACAACACGACAtt            >  1:2589909/1‑35 (MQ=12)
                        cGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGACATc            <  1:281281/36‑1 (MQ=255)
                         gAGTTGGCGGC.GGTTTTTGTTTCCAACAGGACATcc           >  1:6083223/1‑36 (MQ=25)
                         gAGTTGGCGGC.GGTTTTGGTTCCAAACAGGACATcc           <  1:3051068/36‑1 (MQ=255)
                          aGTTGGCGGC.GTTTTTGGTTCCAAACAGGACATCct          <  1:1880547/36‑1 (MQ=255)
                             tGGCGGC.GGTTTTGGTTCCAAACAGGACATCctgct       <  1:3122517/36‑1 (MQ=255)
                             tGGCGGC.GGTTTTGGGTCCAAACAGGACATCctgct       <  1:4327295/36‑1 (MQ=255)
                              ggcggc.GGTTTTTGTTTCCCACAGGACATCctgctt      >  1:1622256/1‑35 (MQ=16)
                                cggc.ggTTTTGGTTCCAAACAGGACATCCTGCTGga    <  1:3260657/36‑1 (MQ=255)
                                 ggc.ggTGTTGGTTCCAAACAGGACATCCTGCTGgaa   <  1:6060697/36‑1 (MQ=255)
                                  gc.ggTTTTGGTTCCAAACAGGACATCCTGCTGgaag  >  1:600441/1‑36 (MQ=255)
                                  gc.ggTTTTGGTTCCAAACAGGACATCCTGCTGgaag  <  1:5240753/36‑1 (MQ=255)
                                  gc.ggTTTTGGTTCCAAACAGGACATCCTGCTGgaag  <  1:3867560/36‑1 (MQ=255)
                                  gc.ggTTTTGGTTCCAAACAGGACATCCTGCTGgaag  <  1:2153666/36‑1 (MQ=255)
                                  gc.ggTTTTGGTTCCAAACAGGACATCCTGCTGgaag  >  1:944448/1‑36 (MQ=255)
                                    |                                  
CATAAAGTTCATGTCATTAAAGAGCGAGTTGGCGGC.GGTTTTGGTTCCAAACAGGACATCCTGCTGGAAG  >  REL606/2908363‑2908432

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 21 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG