Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,134,710 G→T 9.9% intergenic (‑2/‑152) argE ← / → argC acetylornithine deacetylase/N‑acetyl‑gamma‑glutamyl‑phosphate reductase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,134,7100G→T 9.9% 2.4 ‑3.0 0.047intergenic (‑2/‑152)argE/argCacetylornithine deacetylase/N‑acetyl‑gamma‑glutamyl‑phosphate reductase
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (19/23):  total (24/23)
Fisher's exact test strand distribution p-value = 4.96e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTC  >  REL606/4134678‑4134745
                                |                                   
tctcGATAAATGGCGGTAATTTTTTTTTTAtttttt                                  >  1:1725145/1‑35 (MQ=25)
 ctcGATAAATGGCGGTAATTTTTTTTTTAttttttt                                 >  1:6784037/1‑34 (MQ=21)
 ctcGATAAATGGCGGAAATTTTTTTTTTAtttttta                                 >  1:8467592/1‑34 (MQ=16)
  tcGATAAATGGCGGTAATTTGTTTTTCATTGTTGac                                <  1:767350/36‑1 (MQ=255)
   cGATAAATGGCGGTAAATTTTTTTTTATTTTTTact                               >  1:5328125/1‑35 (MQ=12)
    gATAAATGGCGGTTATTTTTTTTTTATTTTTGtcac                              >  1:5369576/1‑36 (MQ=12)
     aTAAATGGCGGTAATTTTTTTTTTATTTTTTacacc                             >  1:181735/1‑35 (MQ=16)
     aTAAATGGCGGTAATTTGTTTTTCATTGTTGacaca                             <  1:1698205/36‑1 (MQ=255)
     aTAAATGGCGGTAATTTGTTTTTCATTGTTGacaca                             <  1:3835993/36‑1 (MQ=255)
     aTAAATGGCGGTAATGTGTTTTTCATTGTTGacaca                             <  1:4365779/36‑1 (MQ=255)
       aaaTGGCGGTAATTTGTTTTTCATTTTTTACACAcc                           >  1:928654/1‑36 (MQ=37)
        aaTGGCGGTAATTTGTTTTTCATTGTTGACACACct                          <  1:7831439/36‑1 (MQ=255)
          tGGCGGTAATTTGTTTTTCATTTTTTACACACctct                        >  1:463945/1‑36 (MQ=37)
          tGGCGGTAATTTGTTTTTCATTGTTGACACACctct                        <  1:5527472/36‑1 (MQ=255)
            gCGGTAATTTGTTTTTCATTGTTGACACACCTCTgg                      >  1:1173207/1‑36 (MQ=255)
             cGGTAATTTGTTTTTCATTGTTGACACACCTCTGGt                     <  1:4660921/36‑1 (MQ=255)
             cGGTAATTTGTTTTTCATTGTTGACACACCTCTGGt                     <  1:3270255/36‑1 (MQ=255)
             cGGTAATTTGTTTTTCATTGTTGACACACCTCTGGt                     <  1:6550783/36‑1 (MQ=255)
              ggTAATTTGTTTTTCATTGTTGACACACCTCTGGTc                    >  1:7170098/1‑36 (MQ=255)
               gTAATTTGTTTTTCATTGTTGACACACCTCTGGTCa                   <  1:7025044/36‑1 (MQ=255)
               gTAATTTGTTTTTCATTGTTGACACACCTCTGGTCa                   >  1:7834653/1‑36 (MQ=255)
                tAATTTGTTTTTCATTGTTGACACACCTCTGGTCag                  <  1:1636187/36‑2 (MQ=255)
                tAATTTGTTTTTCATTGTTGACACACCTCTGGTCAt                  <  1:6893426/36‑1 (MQ=255)
                 aaTTTGTTTTTCATTGTTGACACACCTCTGGTCATg                 <  1:3847555/36‑1 (MQ=255)
                  aTTTGTTTTTCATTGTTGACACACCTCTGGTCATGa                >  1:1920124/1‑36 (MQ=255)
                  aTTTGTTTTTCATTGTTGACACACCTCTGGTCATGa                <  1:4277345/36‑1 (MQ=255)
                  aTTTGTTTTTCATTGTTGACACACCTCTGGTCATGa                <  1:3489747/36‑1 (MQ=255)
                  aTTTGTTTTTCATTGTTGACACACCTCTGGTCATGa                <  1:2498993/36‑1 (MQ=255)
                   tttGTTTTTCATTGTTGACACACCTCTGGTCATGAt               >  1:6310432/1‑36 (MQ=255)
                   tttGTTTTTCATTGTTGACACACCTCTGGTCATGAt               <  1:3631689/36‑1 (MQ=255)
                   tttGTTTTTCATTGTTGACACACCTCTGGTCATGAt               <  1:4928206/36‑1 (MQ=255)
                   tttGTTTTTCATTGTTGACACACCTCTGGTCATGAt               >  1:4973072/1‑36 (MQ=255)
                    ttGTTTTTCATTGTTGACACACCTCTGGTCATGATa              <  1:5748424/36‑1 (MQ=255)
                    ttGTTTTTCATTGTTGACACACCTCTGGTCATGATa              >  1:457234/1‑36 (MQ=255)
                       tttttCATTGTTGACACACCTCTGGTCATGATAGTa           >  1:7111922/1‑36 (MQ=255)
                       tttttCATTGTTGACACACCTCTGGTCATGATAGTa           >  1:1119300/1‑36 (MQ=255)
                       tttttCATTGTTGACACACCTCTGGTCATGATAGTa           <  1:2316729/36‑1 (MQ=255)
                       tttttCATTGTTGACACACCTCTGGTCATGATAGTa           >  1:3391724/1‑36 (MQ=255)
                       tttttCATTGTTGACACACCTCTGGTCATGATAGTa           >  1:464535/1‑36 (MQ=255)
                        ttttCATTGTTGACACACCTCTGGTCATGATAGTAt          >  1:1532469/1‑36 (MQ=255)
                        ttttCATTGTTGACACACCTCTGGTCATGATAGTAt          >  1:2910227/1‑36 (MQ=255)
                         tttCATTGTTGACACACCTCTGGTCATGATAGTATc         <  1:3233408/36‑1 (MQ=255)
                         tttCATTGTTGACACACCTCTGGTCATGATAGTATc         <  1:2190104/36‑1 (MQ=255)
                          ttCATTGTTGACACACCTCTGGTCATGATAGTATCa        >  1:213317/1‑36 (MQ=255)
                           tCATTGTTGACACACCTCTGGTCATGATAGTATCaa       >  1:2481805/1‑36 (MQ=255)
                            cATTGTTGACACACCTCTGGTCATGATAGTATCAat      >  1:3487535/1‑36 (MQ=255)
                            cATTGTTGACACACCTCTGGTCATGATAGTATCAat      >  1:7644261/1‑36 (MQ=255)
                            cATTGTTGACACACCTCTGGTCATGATAGTATCAat      <  1:2010053/36‑1 (MQ=255)
                             aTTGTTGACACACCTCTGGTCATGATAGTATCAata     >  1:8148545/1‑36 (MQ=255)
                             aTTGTTGACACACCTCTGGTCATGATAGTATCAata     >  1:409883/1‑36 (MQ=255)
                              ttgttgACACACCTCTGGTCATGATAGTATCAatat    >  1:6424190/1‑36 (MQ=255)
                               tgttgACACACCTCTGGTCATGATAGTATCAATAtt   <  1:6472105/36‑1 (MQ=255)
                               tgttgACACACCTCTGGTCATGATAGTATCAATAtt   <  1:6501416/36‑1 (MQ=255)
                               tgttgACACACCTCTGGTCATGATAGTATCAATAtt   <  1:4508468/36‑1 (MQ=255)
                               tgttgACACACCTCTGGTCATGATAGTATCAATAtt   <  1:576956/36‑1 (MQ=255)
                                gttgACACACCTCTGGTCATGATAGTATCAATATTc  >  1:6830766/1‑36 (MQ=255)
                                gttgACACACCTCTGGTCATGATAGTATCAATATTc  <  1:4590438/36‑1 (MQ=255)
                                |                                   
TCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTC  >  REL606/4134678‑4134745

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG