Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,835,688 Δ1 bp 6.1% intergenic (‑167/‑111) dgoR ← / → yidX predicted DNA‑binding transcriptional regulator/predicted lipoproteinC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,835,6880A→. 6.1% 2.2 ‑0.7 ‑0.133intergenic (‑167/‑111)dgoR/yidXpredicted DNA‑binding transcriptional regulator/predicted lipoproteinC
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (14/17):  total (14/19)
Fisher's exact test strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.26e-01

TACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAG  >  REL606/3835654‑3835723
                                  |                                   
tactacAATTTAGATCACAAAAAGAACAATGCATaa                                    >  1:5133828/1‑36 (MQ=255)
tactacAATTTAGATCACAAAAAGAACAATGCATaa                                    >  1:2172828/1‑36 (MQ=255)
 actacAATTTAGATCACAAAAAGAACAATGCATaaa                                   >  1:7145455/1‑36 (MQ=255)
  ctacAATTTAGATCACAAAAAGAACAATGCATaaaa                                  >  1:5069134/1‑36 (MQ=255)
  atacAATTTAGATCACAAAAAGAACAATGCATaaaa                                  <  1:4561070/35‑1 (MQ=255)
   tacAATTTAGATCACAAAAAGAACAATGCATaaaaa                                 >  1:738476/1‑36 (MQ=255)
    acAATTTAGATCACAAAAAGAACAATGCATaaaaaa                                <  1:310210/36‑1 (MQ=255)
    acAATTTAGATCACAAAAAGAACAATGCATaaaaaa                                >  1:1642512/1‑36 (MQ=255)
    acAATTTAGATCAAAAAAAGAACAATGCATaaaaaa                                <  1:3168032/36‑1 (MQ=255)
          tAGATCACAAAAAGAACAATGCATAAAAAAATGACa                          >  1:385125/1‑36 (MQ=255)
           aGATCACAAAAAGAACAATGCATAAAAAAATGACAt                         <  1:1325312/36‑1 (MQ=255)
           aGATCAAAAAAAGAACAATGCATAAAAAAATGACAt                         <  1:2198412/36‑1 (MQ=255)
            gATCACAAAAAGAACAATGCATAAAAAAATGACATg                        >  1:776714/1‑36 (MQ=255)
            gATCACAAAAAGAACAATGCATAAAAAAATGACATg                        >  1:3479529/1‑36 (MQ=255)
              tCACAAAAAGAACAATGCATAAAAAAATGACATGCg                      <  1:3865494/36‑1 (MQ=255)
               cacaAAAAGAACAATGCATAAAAAAATGACATGCGt                     >  1:2186159/1‑36 (MQ=255)
               cacaAAAAGAACAATGCATAAAAAAATGACATGCGt                     <  1:3982072/36‑1 (MQ=255)
               aacaAAAAGAACAATGCATAAAAAAATGACATGCGt                     <  1:7138082/35‑1 (MQ=255)
                acaAAAAGAACAATGCATAAAAAAATGACATGCGTc                    <  1:6390490/36‑1 (MQ=255)
                 caAAAAGAACAATGCATAAAAAAATGACATGCGTCg                   >  1:5628198/1‑36 (MQ=255)
                  aaaaaGAACAATGCATAAAAAAATGACATGCGTCgg                  >  1:779620/1‑36 (MQ=255)
                   aaaaGAACAATGCATAAAAAAATGACATGCGTCggg                 >  1:5264672/1‑36 (MQ=255)
                   aaaaGAACAATGCATAAAAAAATGACATGCGTCggg                 <  1:7277717/36‑1 (MQ=255)
                    aaaGAACAATGCATAAAAAAATGACATGCGTCGGGc                <  1:3062586/36‑1 (MQ=255)
                    aaaGAACAATGCATAAAAAAATGACATGCGTCGGGc                <  1:3895041/36‑1 (MQ=255)
                       gAACAATGCATAAAAAAATGACATGCGTCGGGCAGa             >  1:4544209/1‑36 (MQ=255)
                        aaCAATGCATAAAAAAATGACATGCGGCGGGCAGaa            >  1:6125267/1‑36 (MQ=255)
                         aCAATGCATAAAAAAATGACATGCGTCGGGCAGaaa           >  1:1397927/1‑36 (MQ=255)
                         aCAATGCATAAAAAAATGACATGCGTCGGGCAGaaa           >  1:1258283/1‑36 (MQ=255)
                         aCAATGCATAAAAAAATGACATGCGTCGGGCAGaaa           >  1:1920761/1‑36 (MQ=255)
                         aCAATGCATAAAAAAATGACATGCGTCGGGCAGaaa           <  1:3093857/36‑1 (MQ=255)
                          cAATGCATAAAAAAATGACATGCGTCGGGCAGAAAt          <  1:4156886/36‑1 (MQ=255)
                           aaTGCATAAAAAAATGACATGCGTCGGGCAGAAATc         <  1:6224270/36‑1 (MQ=255)
                           aaTGCATAAAAAAATGACATGCGTCGGGCAGAAATc         >  1:1904112/1‑36 (MQ=255)
                            aTGCAT.AAAAAATGACATGCGTCGGGCAGAAATCTg       <  1:2954121/36‑1 (MQ=37)
                             tGCATAAAAAAATGACATGCGTCGGGCAGAAATCTg       <  1:4658945/36‑1 (MQ=255)
                             tGCAT.AAAAAATGACATGCGTCGGGCAGAAATCTGa      <  1:1727262/36‑1 (MQ=37)
                              gCATAAAAAAATGACATGCGTCGGGCAGAAATCTGa      <  1:2959430/36‑1 (MQ=255)
                               cATAAAAAAATGACATGCGTCGGGCAGAAATCTGaa     <  1:867251/36‑1 (MQ=255)
                                aTAAAAAAATGACATGCGTCGGGCAGAAATCTGaaa    <  1:4969671/36‑1 (MQ=255)
                                aTAAAAAAATGACATGCGTCGGGCAGAAATCTGaaa    >  1:1113344/1‑36 (MQ=255)
                                 tAAAAAAATGACATGCGTCGGGCAGAAATCTGaaaa   <  1:5359425/36‑1 (MQ=255)
                                  aaaaaaaTGACATGCGTCGGGCAGAAATCTGAAAAg  >  1:1045388/1‑36 (MQ=255)
                                  aaaaaaaTGACATGCATCGGGCACAAATCTGtaaaa  >  1:2918281/1‑35 (MQ=25)
                                  aaaaaaaAGACATGCGTCGGGCAGAAATCTGAAAAg  <  1:1518358/36‑1 (MQ=255)
                                  |                                   
TACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAG  >  REL606/3835654‑3835723

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG