Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,856,086 G→T 20.5% 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→T20.5% 2.2 ‑1.3 ‑0.120intergenic (+117/+150)yeaP/yeaQpredicted diguanylate cyclase/conserved inner membrane protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (10/6):  total (14/6)
Fisher's exact test strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.99e-01

ATTTTTAGCATGGTTTTTTATTGGCGGC.TATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATT  >  REL606/1856052‑1856120
                                   |                                  
aTTTTTAGCATGGTTTTTTATTGGCGGC.TAtgctgc                                   <  1:911565/36‑1 (MQ=255)
  ttttAGCATGGTTTTTTATTGGCGGG.TATTCTTccc                                 >  1:1722787/1‑36 (MQ=25)
       gCATGGTTTTTTATTGTCGGC.TATGCTGCCCCGGGa                            <  1:7613440/36‑1 (MQ=255)
        cATGGTTTTTTATTGGCGGC.TATTTTTCCCCggggg                           >  1:944223/1‑34 (MQ=21)
        cATGGTTTTTTATTGAGGGC.TATTCTGCCCCGGGAg                           <  1:44918/36‑1 (MQ=25)
         aTGGTTTTTTATTGGCGGG.TATTTTGCCCCgggggg                          >  1:6935900/1‑33 (MQ=16)
         aTGGTTTTTTATTGGCGGC.TTTTCTTCCCCcggggg                          >  1:5806240/1‑33 (MQ=16)
          tGGTTTTTTATTGGCGGCTTTTGTTGCCCCGGGGGc                          >  1:406998/1‑36 (MQ=16)
          tGGTTTTTTATTGGCGGC.TATTCTTCCCCGGGGGcc                         >  1:192826/1‑35 (MQ=25)
              tttttATTGGCGGC.TATGCTGCCCCGGGAGCATaaa                     <  1:2273180/36‑1 (MQ=255)
               ttttATTGGCGGC.TATGCTGCCCCGGGAGCATAAAg                    <  1:69512/36‑1 (MQ=255)
                     tGGCGGC.TATGCTGCCCCGGGAGCATAAAGATGaaa              >  1:1098277/1‑36 (MQ=255)
                     tGGCGGC.TATGCTGCCCCGGGAGCATAAAGATGaaa              >  1:7851585/1‑36 (MQ=255)
                      ggcggc.TATGCTGCCCCGGGAGCATAAAGATGaaaa             <  1:6916277/36‑1 (MQ=255)
                      ggcggc.TATGCTGCCCCGGGAGCATAAAGATGaaaa             >  1:398361/1‑36 (MQ=255)
                      ggcggc.TATGCTGCCCCGGGAGCATAAAGATGaaaa             >  1:2726620/1‑36 (MQ=255)
                         ggc.TATGCTGCCCCGGGAGCATAAAGATGaaaaaaa          >  1:314873/1‑36 (MQ=255)
                          gc.TATGCTGCCCCGGGAGCATAAAGATGAAAAAaac         <  1:4422011/36‑1 (MQ=255)
                          gc.TATGCTGCCCCGGGAGCATAAAGATGAAAAAaac         >  1:6733799/1‑36 (MQ=255)
                             tATGCTGCCCCGGGAGCATAAAGATGAAAAAaacaa       >  1:1942961/1‑36 (MQ=255)
                             tATGCTGCCCCGGGAGCATAAAGATGAAAAAaacaa       >  1:1494457/1‑36 (MQ=255)
                                gctgcCCCGGGAGCATAAAGATGAAAAAAACAACGa    >  1:2487621/1‑36 (MQ=255)
                                 ctgcCCCGGGAGCATAAAGATGAAAAAAACAACGat   >  1:1801445/1‑36 (MQ=255)
                                 ctgcCCCGGGAGAATAAAGATGAAAAAAACAACGat   <  1:3257572/36‑1 (MQ=255)
                                  tgcCCCGGGAGCATAAAGATGAAAAAAACAACGatt  >  1:7460498/1‑36 (MQ=255)
                                   |                                  
ATTTTTAGCATGGTTTTTTATTGGCGGC.TATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATT  >  REL606/1856052‑1856120

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG