Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,509,366 +A 3.6% coding (814/996 nt) gntR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,509,3661.→A 3.6% 2.2 ‑1.9 ‑0.056coding (814/996 nt)gntRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (34/20):  total (34/22)
Fisher's exact test strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.19e-01

CATCACCTGACCAATGTCATGACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGAACTTTTAACCCCAGAC  >  REL606/3509331‑3509401
                                     |                                   
cATCACCTGACCAATGTCATGACCGTGGAAA.CCGGc                                      >  1:224710/1‑36 (MQ=255)
cATCACCTGACCAATGTCATGACCGTGGAAA.CCGGc                                      >  1:2736605/1‑36 (MQ=255)
 aTCACCTGACCAATGTCATGACCGTGGAAA.CCGGC.a                                    >  1:6382916/1‑36 (MQ=255)
 aTCACCTGACCAATGTCATGACCGTGGAAA.CCGGC.a                                    >  1:6381819/1‑36 (MQ=255)
 aTCACCTGACCAACGACATGACCCTGGAAA.CCGGC.a                                    >  1:5275683/1‑36 (MQ=21)
 aTCACCTGACAAATGTCATGACCGTGGAAA.CCGGC.a                                    <  1:5223238/36‑1 (MQ=255)
   cACCTGACCAATGTCATGACCGTGGAAA.CCGGC.AAt                                  <  1:3039098/36‑1 (MQ=255)
    aCCTGACCACTGCCATGACCGTGGAAA.CCGGC.AATc                                 <  1:7436392/36‑1 (MQ=37)
    aCCTGACCAATGTCATGACCGTGGAAA.CCGGC.AATc                                 <  1:1801677/36‑1 (MQ=255)
    aCCTGACCAATGTCATGACCGTGGAAA.CCGGC.AATc                                 <  1:4949268/36‑1 (MQ=255)
    aCCTGACCAATGTCATGACCGTGGAAA.CCGGC.AATc                                 >  1:4406842/1‑36 (MQ=255)
    aCCTGACAAATGTCATGACCGTGGAAA.CCGGC.AATc                                 <  1:3452599/36‑1 (MQ=255)
     ccTGACCAATGTCATGACCGTGGAAA.CCGGC.AATCg                                >  1:6178695/1‑36 (MQ=255)
       tGACCAATGTCATGACCGTGGAAA.CCGGC.AATCGcc                              <  1:933776/36‑1 (MQ=255)
       tGACCAATGTCATGACCGTGGAAA.CCGGC.AATCGcc                              <  1:8353783/36‑1 (MQ=255)
       tGACCAATGTCATGACCGTGGAAA.CCGGC.AATCGcc                              <  1:537133/36‑1 (MQ=255)
        gACCAATGTCATGACCGTGGAAA.CCGGC.AATCGCCa                             >  1:2384125/1‑36 (MQ=255)
         aCCAATGTCATGACCGTGGAAA.CCGGC.AATCGCCat                            >  1:3130735/1‑36 (MQ=255)
          ccAATGTCATGACCATGGAAA.CCGGC.AATCGCCata                           >  1:6506713/1‑36 (MQ=255)
           cAATGTCATGACCGTGGAAA.CCGGC.AATCGCCatat                          <  1:6301403/36‑1 (MQ=255)
            aaTGTCATGACCGTGGAAA.CCGGT.AATCGCCATATc                         >  1:4925380/1‑36 (MQ=255)
             aTGTCATGACCGTGGAAA.CCGGC.AATCGCCATATCg                        <  1:2753711/36‑1 (MQ=255)
             aTGTCAGGACCGGGGAAA.CCGGA.AATCGCCATATCg                        <  1:6273300/36‑1 (MQ=25)
              tGTCATGACCGTGGAAA.CCGGC.AATCGCCATATCGt                       >  1:3080866/1‑36 (MQ=255)
              tGTCATGACCGTGGAAA.CCGGC.AATCGCCATATCGt                       >  1:6454602/1‑36 (MQ=255)
               gTCATGACCGTGGAAA.CCGGC.AATCGCCATATCGTc                      >  1:8358278/1‑36 (MQ=255)
                 cATGACCGTGGAAA.CCGGC.AATCGCCATATCGTCAg                    >  1:4210450/1‑36 (MQ=255)
                 aaTGACCGTGAAAA.CCGGA.AATCGCCATATCGTCAg                    <  1:733731/35‑1 (MQ=37)
                  aTGACCGTGGAAA.CCGGC.AATCGCCATATCGCCAgg                   >  1:7066913/1‑36 (MQ=255)
                   tGACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGa                  >  1:6269007/1‑36 (MQ=255)
                   tGACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGa                  >  1:6787026/1‑36 (MQ=255)
                   tGACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGa                  >  1:4553279/1‑36 (MQ=255)
                    gACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGaa                 >  1:8007028/1‑36 (MQ=255)
                    gACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGaa                 >  1:7184465/1‑36 (MQ=255)
                    gACCGTGGAAA.CCGGC.AAGCGCCATATCGCCAGGaa                 >  1:4230010/1‑36 (MQ=37)
                    gACCGTGAAAA.CCGGA.AATCGCCATATCGTCAGGaa                 <  1:141405/36‑1 (MQ=37)
                    gACCGGGAAAA.CCGGCAAATCGCCATATCGTCAGGa                  <  1:6007151/36‑1 (MQ=16)
                    acccGGGGAAACCCGGA.AATCGCCATATCGTCAGGa                  <  1:4459204/34‑1 (MQ=16)
                      ccGTGGAAA.CCGGC.AATCGCCATATCGTCAGGAACt               >  1:5928307/1‑36 (MQ=255)
                       cggggAAA.ACGGCAAATCGCCATATCGTCAGGAACt               <  1:6043263/33‑1 (MQ=16)
                       cGTGGAAA.CCGGC.AATCGCCATATCGTCAGGAACtt              >  1:4450884/1‑36 (MQ=255)
                       cGTGGAAA.CCGGC.AATCGCCATATCGTCAGGAACtt              >  1:5592037/1‑36 (MQ=255)
                        ggggAAA.CCGGA.AATCGCCATATCGTCAGGAACttt             <  1:8320563/34‑1 (MQ=37)
                        gTGGAAA.CCGGC.AATCGCCATATCGTCAGGAACttt             >  1:776201/1‑36 (MQ=255)
                        gTGGAAA.CCGGA.AATCGCCATATCGTCAGGAACttt             <  1:1081526/36‑1 (MQ=255)
                         tGGAAA.CCGGC.AATCGCCATATCGTCAGGAACtttt            >  1:4003278/1‑36 (MQ=255)
                         tGGAAA.CCGGC.AATCGCCATATCGTCAGGAACtttt            >  1:1647565/1‑36 (MQ=255)
                         tGGAAA.CCGGC.AATCGCCATATCGTCAGGAACtttt            >  1:2825963/1‑36 (MQ=255)
                           gAAA.CCGGC.AATCGCCATATCGTCAGGAACTTTTaa          >  1:4008998/1‑36 (MQ=255)
                           aaaa.CCGGC.AATCGCCATATCGTCAGGAACTTTTaa          <  1:813990/35‑1 (MQ=255)
                           aaaa.CCGGA.AATCGCCATATCGTCAGGAACTTTTaa          <  1:6166201/35‑1 (MQ=255)
                            aaa.ACGGC.AATCGCCATATCGTCATGCATTTTTAAc         >  1:4501239/1‑36 (MQ=19)
                              a.CCGGC.AATCGCCATATCGTCAGGAACTTTTAAccc       >  1:5692045/1‑36 (MQ=255)
                              a.CCGGC.AATCGCCATATCGTCAGGAACTTTTAAccc       >  1:2356224/1‑36 (MQ=255)
                              a.CCGGC.AATCGCCATATCGTCAGGAACTTTTAAccc       >  1:1255517/1‑36 (MQ=255)
                                ccGGC.AATCGCCATATCGTCAGGAACTTTTAAcccc      >  1:6625160/1‑36 (MQ=255)
                                ccGAC.AATCGCCATATCGTCAGGAACTTTTAAcccc      <  1:8344489/36‑1 (MQ=255)
                                 gggC.AATCGCCATATCGTCAGGAACTTTTAACCCCa     <  1:6870861/35‑1 (MQ=255)
                                 cgac.AATCGCCATATCGTCAGGAACTTTTAACCCCa     <  1:1547546/33‑1 (MQ=255)
                                 cGGC.AATCGCCATATCGTCAGGAACTTTTAACCCCa     >  1:4546176/1‑36 (MQ=255)
                                  gac.AATCGCCATATCGTCAGGAACTTTTAACCCCAg    <  1:7721115/34‑1 (MQ=255)
                                   gC.AATCGCCATATCGTCAGGAACTTTTAACCCCAGa   >  1:448403/1‑36 (MQ=255)
                                   gC.AATCGCCATATCGTCAGGAACTTTTAACCCCAGa   >  1:1351975/1‑36 (MQ=255)
                                   gC.AATCGCCATATCGTCAGGAACTTTTAACCCCAGa   >  1:1355974/1‑36 (MQ=255)
                                    c.AATCGCCATATCGTCAGGAACTTTTAACCCCAGAc  >  1:3868399/1‑36 (MQ=255)
                                     |                                   
CATCACCTGACCAATGTCATGACCGTGGAAA.CCGGC.AATCGCCATATCGTCAGGAACTTTTAACCCCAGAC  >  REL606/3509331‑3509401

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG