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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,856,1450G→T25.5% 4.2 ‑4.1 ‑0.043intergenic (+176/+91)yeaP/yeaQpredicted diguanylate cyclase/conserved inner membrane protein
Reads supporting (aligned to +/- strand):  new base (11/1):  ref base (16/15):  total (27/16)
Fisher's exact test strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.65e-01

AAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGT  >  REL606/1856110‑1856180
                                   |                                   
aaaCAACGATTATTATGATGGGTGTGGCGATTATTg                                     <  1:4799292/36‑1 (MQ=255)
aaaCAACGATTATTATGATGGGTGTGGCGATTATTg                                     >  1:4556708/1‑36 (MQ=255)
 aacaacGATTATTATGATGGGTGTGGCGATTATTGt                                    <  1:1959031/36‑1 (MQ=255)
    aacGATTATTATGATGGGTGTGGCGATTATTGTCGt                                 <  1:856280/36‑1 (MQ=255)
    aacGATTATTATGATGGGTGTGGCGATTATTGTCGt                                 >  1:348705/1‑36 (MQ=255)
    aacGATTATTAAGATGCGTGTCGCGATTATTGTCGt                                 <  1:5915469/36‑1 (MQ=21)
     acGATTATTATGATGGGTGTGGCTATTATTTTCGTa                                <  1:4008794/36‑1 (MQ=38)
     acGATTATTATGATGGGTGTGGCGATTATTGTCGTa                                <  1:1292183/36‑1 (MQ=255)
       aattattATGATGGGTGTGGCGATTATTTTTGTagt                              >  1:1493955/2‑34 (MQ=25)
        attattATGATGGGTGTGGCGATTATTTTTGTattt                             >  1:1860735/1‑33 (MQ=25)
          tattatGATGGGTGTGGCGATTTTTTTTGTACTAgg                           >  1:4432617/1‑36 (MQ=16)
          tattatGATGGGTGTGGCGATTATTTTTGTACTctg                           >  1:2082779/1‑34 (MQ=25)
          tattatGATGGGTGTGGCGATTATTGTCGGACTCgg                           <  1:2409399/36‑1 (MQ=255)
           attatGATGGGTGTGGCGATTATTTTCGTACTTGGc                          >  1:5132368/1‑36 (MQ=37)
           attatGATGGGTGTGGCGATTATTGTCGTACTCGGc                          <  1:91481/36‑1 (MQ=255)
            ttatGATGGGTGTGGCGATTATTTTTGTACTGGGCa                         >  1:5622395/1‑36 (MQ=25)
            ttatGATGGGTGCGGCGATTATTGTCGTACTCGGCa                         >  1:4314467/1‑36 (MQ=255)
             tatGATGGGTGTGGCGATTATTTTTGTACTCGGccc                        >  1:2223801/1‑34 (MQ=25)
             tatGATGGGTGTGGCGATTATTGTCGTACTCGGccc                        >  1:252799/1‑34 (MQ=255)
             gatgatgGGAGTGGCTATTATTGTCGTACTCGGCAc                        <  1:617577/35‑1 (MQ=37)
               tgatgGGTGTGGCGATTATTTTCGTAATCGGCACTg                      >  1:5564351/1‑36 (MQ=37)
               tgatgGGTGTGGCGATTATTGTCTTACTCTGCattt                      >  1:1408114/1‑33 (MQ=25)
                gatgGGTGTGGCGATTCTTTTCGTAATCGGCACttt                     >  1:108216/1‑34 (MQ=21)
                 atgGGTGTGGCGATTATTTTTGTACTCGGCACTTAg                    >  1:1269028/1‑36 (MQ=25)
                 atgGGTGTGGCGATTATTTTAGTACTCGGCACttgg                    >  1:5548037/1‑33 (MQ=25)
                 atgGGTGTGGCGATTATTGTCGTACTCGGCACTGAg                    >  1:4324094/1‑36 (MQ=255)
                  tgGGTGTGGCGATTATTGTCGTACTCGGCACTGAGc                   <  1:1683064/36‑1 (MQ=255)
                  tgGGTGTGGCGATTATTGTCGTACTCGGCACTGAGc                   >  1:3958246/1‑36 (MQ=255)
                  tgGGTGTGGCAATTATTGTCGTACTCGGCACTGAGc                   <  1:863409/36‑1 (MQ=255)
                   gggTGTTGCGATTATTGTCGTACTCGGCGCTGAGCt                  <  1:1086489/36‑1 (MQ=37)
                   gggTGTGGCGATTATTGTCGTACTCGGCACTGGGCt                  >  1:5380582/1‑36 (MQ=255)
                   gggTGTGGCGATTATTGTCGTACTCGGCACTGAggt                  >  1:2817142/1‑34 (MQ=255)
                   gggTGTGGCGATTATTGTCGTACTCGGCACTGAGCt                  <  1:1328958/36‑1 (MQ=255)
                   gggTGTGGCGATTATTGTCGTAATCGGCATTGAGCt                  >  1:2641599/1‑36 (MQ=37)
                    ggTGCGGCGATTATTGCCGTAATCGGAACTCAGccg                 >  1:551418/1‑34 (MQ=9)
                     gtgtGGCGATTATTGTCGTAATAGGCACTGAGCTgg                >  1:4036078/1‑36 (MQ=37)
                       gtgGCGATTATTGTCGTACTCGGCACTGAGCTGGGa              >  1:5423904/1‑36 (MQ=255)
                        tgGCGATTATTGTCGTACTCGGCACTGAGCTGGGAt             >  1:5452976/1‑36 (MQ=255)
                         ggCGATTATTGTCGTACTCGGCACTGAGCTGGGAtg            >  1:3975326/1‑36 (MQ=255)
                            gATTATTGTCGTACTCGGCACTGAGCTGGGAtggtg         <  1:4631436/36‑1 (MQ=255)
                             attattGTCGTACTCGGCACTGAGCTGGGAtggtgg        >  1:4662088/1‑36 (MQ=255)
                              ttattGTCGTACTCGGCACTGAGCTGGGATggtggt       <  1:5225587/36‑1 (MQ=255)
                               tattGTCGTACTCGGCACTGAGCTGGGATGGTGGTa      >  1:1408428/1‑36 (MQ=255)
                                attGTCGTACTCGGCACTGAGCTGGGATGGTGGTaa     <  1:3342176/36‑1 (MQ=255)
                                 ttGTCGTACTCGGCACTGAGCTGGGATGGTGGTAAc    <  1:1374542/36‑1 (MQ=255)
                                   gTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGt  <  1:3691998/36‑1 (MQ=255)
                                   gTCGTACTCGGCAATGAGCAGGGATGGTGGTAACGt  >  1:5287406/1‑36 (MQ=38)
                                   |                                   
AAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGT  >  REL606/1856110‑1856180

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 37 ≤  ATCG