Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 88,137 A→G 8.9% intergenic (+21/‑297) leuO → / → ilvI leucine transcriptional activator/acetolactate synthase III large subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60688,1370A→G 8.9% 2.9 ‑1.9 ‑0.345intergenic (+21/‑297)leuO/ilvIleucine transcriptional activator/acetolactate synthase III large subunit
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (17/23):  total (17/28)
Fisher's exact test strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.64e-01

ATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTAT  >  REL606/88102‑88169
                                   |                                
aTTTGCAAACGCTAACTGATTGCAGAATAGGTCAGa                                  >  1:2378495/1‑36 (MQ=255)
 tttGCAAACGCTAACTGATTGCAGAATAGGTCAGAc                                 <  1:102885/36‑1 (MQ=255)
 tttGCAAACGCTAACTGATTGCAGAATAGGTCAGAc                                 >  1:7365567/1‑36 (MQ=255)
 tttGCAAACGCTAACTGATCGCAGAATAGGTCAGAc                                 <  1:4098159/36‑1 (MQ=255)
  ttGCAAACGCTAACTGATTGCAGAATAGGGCAGACa                                >  1:482840/1‑36 (MQ=255)
   tGCAAACGCTAACTGATTGCAGAATAGGTCAGCCAt                               <  1:745332/36‑1 (MQ=255)
    gCAAACGCTAACTGATTGCAGAATAGGTCAGACATg                              >  1:7731386/1‑36 (MQ=255)
    gCAAACGCTAACTGATTGCAGAATAGGTCAGACATg                              >  1:1774329/1‑36 (MQ=255)
    gCAAACGCTAACTGATTGCAGAATAGGTCAGACATg                              >  1:103061/1‑36 (MQ=255)
     cAAACGCTAACTGATTGCAGAATAGGTCAGACATGa                             <  1:4530520/36‑1 (MQ=255)
       aaCGCTAACTGATTGCAGCATAGGTCAGACATGAAt                           <  1:7413748/36‑1 (MQ=255)
       aaCGCTAACTGATTGCAGAATAGGGCAGACATGAAt                           >  1:5811040/1‑36 (MQ=255)
       aaCGCTAACTGATTGCAGAATAGGACAGACATGAAt                           >  1:5736267/1‑36 (MQ=255)
          gCTAACTGATTGCAGAATAGGTCAGGCATGAATGTc                        <  1:2008943/36‑1 (MQ=255)
           cTAACTGATTGCAGAATAGGTCAGGCATGAATGTCt                       <  1:7322113/36‑1 (MQ=255)
           cTAACTGATTGCAGAATAGGTCAGACATGAATGTCt                       <  1:361788/36‑1 (MQ=255)
           cTAACTGATTGCAGAATAGGTCAGACATGAATGTCt                       <  1:3571346/36‑1 (MQ=255)
             aaCTGATTGCAGAATAGGTCAGACATGAATGTCTgg                     >  1:1277266/1‑36 (MQ=255)
             aaCTGATTGCAGAATAGGTCAGACATGAATGTCTgg                     <  1:6200159/36‑1 (MQ=255)
              aCTGATTGCAGAATAGGTCAGACATGAATGTCTGGt                    <  1:5070592/36‑1 (MQ=255)
               cTGATTGCAGAATAGGTCAGACATGAATGTCTGGtt                   <  1:2137032/36‑1 (MQ=255)
                tGATTGCAGAATAGGTCAGACATGAATGTCTGGttt                  >  1:1627866/1‑36 (MQ=255)
                 gATTGCAGAATAGGTCAGACATGAATGTCTGGTTTa                 <  1:859177/36‑1 (MQ=255)
                 gATTGCAGAATAGGTCAGACATGAATGTCTGGTTTa                 <  1:4029636/36‑1 (MQ=255)
                  aTTGCAGAATAGGTCAGACATGAATGTCTGGTTTAt                >  1:7321903/1‑36 (MQ=255)
                   ttGCAGAATAGGTCAGGCATGAATGTCTGGTTTAtt               <  1:6213382/36‑1 (MQ=255)
                   ttGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt               <  1:6382649/36‑1 (MQ=255)
                   ttGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt               <  1:7858463/36‑1 (MQ=255)
                   ttGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt               <  1:1437963/36‑1 (MQ=255)
                   ttGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt               <  1:7638690/36‑1 (MQ=255)
                    tGCAGAATAGGTCAGACATGAATGTCTGGTTTATTc              <  1:5424614/36‑1 (MQ=255)
                    tGCAGAATAGGTCAGACATGAATGTCTGGTTTATTc              >  1:2341867/1‑36 (MQ=255)
                     gCAGAATAGGTCAGACATGAATGTCTGGTTTATTCt             >  1:7733704/1‑36 (MQ=255)
                     gCAGAATAGGTCAGACATGAATGTCTGGTTTATTCt             <  1:3654221/36‑1 (MQ=255)
                       aGAATAGGTCAGACATGAATGTCTGGTTTATTCTGc           >  1:6455468/1‑36 (MQ=255)
                        gAATAGGTCAGACATGAATGTCTGGTTTATTCTGCa          >  1:4674742/1‑36 (MQ=255)
                         aaTAGGTCAGACATGAATGTCTGGTTTATTCTGCAt         <  1:1503840/36‑1 (MQ=255)
                          aTAGGTCAGACATGAATGTCTGGTTTATTCTGCAtt        <  1:972473/36‑1 (MQ=255)
                             ggTCAGACATGAATGTCTGGTTTATTCTGCAttttt     <  1:7711401/36‑1 (MQ=255)
                             ggTCAGACATGAATGTCTGGTTTATTCTGCAttttt     >  1:1216450/1‑36 (MQ=255)
                              gTCAGACATGAATGTCTGGTTTATTCTGCAtttttt    <  1:7632689/36‑1 (MQ=255)
                              gTCAGACATGAATGTCTGGTTTATTCTGCAtttttt    <  1:5811985/36‑1 (MQ=255)
                              gTCAGACATGAATGTCTGGTTTATTCTGCAtttttt    >  1:5757570/1‑36 (MQ=255)
                               tCAGACATGAATGTCTGGTTTATTCTGCATTTTTTa   <  1:7360895/36‑1 (MQ=255)
                                cAGGCATGAATGTCTGGTTTATTCTGCATTTTTTAt  <  1:7859840/36‑1 (MQ=255)
                                cAGACATGAATGTCTGGTTTATTCTGCATTTTTTAt  >  1:5791474/1‑36 (MQ=255)
                                   |                                
ATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTAT  >  REL606/88102‑88169

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG