Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,856,086 G→T 27.4% intergenic (+117/+150) yeaP → / ← yeaQ predicted diguanylate cyclase/conserved inner membrane protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,856,0860G→T27.4% 2.4 ‑2.3 ‑0.118intergenic (+117/+150)yeaP/yeaQpredicted diguanylate cyclase/conserved inner membrane protein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (6/7):  total (11/7)
Fisher's exact test strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.43e-01

ATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTA  >  REL606/1856052‑1856121
                                  |                                   
aTTTTTAGCATGGTTTTTTATTGGCGGCTAtgctgc                                    <  1:7660821/36‑1 (MQ=255)
  ttttAGCATGGTTTTTTATTGGCGGCTATGCTGccc                                  <  1:7444582/36‑1 (MQ=255)
     tAGCATGGTTTTTTATTGGCGGCTATTCTTcccccg                               >  1:2666897/1‑34 (MQ=37)
        cATGGTTTTTTATTGGCGGGTATTCTTTCCCggggg                            >  1:6910312/1‑34 (MQ=16)
        cATGGTTTTTTATTGGCGGCTATGTTTCCCCcgggg                            >  1:3150088/1‑34 (MQ=18)
           ggTTTTTTATTGGCGGCTATTTCTCCCCGGGAgggt                         >  1:6894207/1‑33 (MQ=12)
           ggTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCAt                         <  1:7806357/36‑1 (MQ=255)
            gTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATa                        <  1:3136202/36‑1 (MQ=255)
               ttttATTGGCGGCTATTCTTCCCCCGGGGCATAAAg                     >  1:58007/1‑36 (MQ=16)
                 ttATTGGCGGCTATGCTGCCCCGGGAGCATAAAGAt                   <  1:895445/36‑1 (MQ=255)
                   aTTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGa                 <  1:1954186/36‑1 (MQ=255)
                    ttGGCGGCTATGCTGCCCCGGGAGCATAAAGATGaa                <  1:4412054/36‑1 (MQ=255)
                      ggcggcTATGCTGCCCCGGGAGCATAAAGATGaaaa              >  1:5815602/1‑36 (MQ=255)
                      ggcggcTATGCTGCCCCGGGAGCATAAAGATGaaaa              >  1:5404119/1‑36 (MQ=255)
                        cggcTATGCTGCCCCGGGAGCATAAAGATGaaaaaa            >  1:135163/1‑36 (MQ=255)
                        aggcTATGATGCCCCGGGAGCATAAAGATGaaaaaa            <  1:3833577/35‑1 (MQ=255)
                         ggcTATGCTGCCCCGGGAGCATAAAGATGaaaaaaa           >  1:2870882/1‑36 (MQ=255)
                          gcTATGCTGCCCCGGGAGCATAAAGATGAAAAAaac          >  1:2484928/1‑36 (MQ=255)
                             aTGCTGCCCCGGGAGCATAAAGATGAAAAAaacaac       >  1:227485/1‑36 (MQ=255)
                              tgctgcCCCGGGAGCATAAAGATGAAAAAAACAACg      >  1:5830513/1‑36 (MQ=255)
                                ctgcCCCGGGAGCATAAAGATGAAAAAAACAACGat    >  1:1633640/1‑36 (MQ=255)
                                 tgcCCCGGGAGCATAACGATGAAAAAAACAACGatt   >  1:4561573/1‑36 (MQ=255)
                                  gcCCCGGGAGCATAAAGATGAAAAAAACAACGatta  >  1:7526419/1‑36 (MQ=255)
                                  |                                   
ATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTA  >  REL606/1856052‑1856121

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG