Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,212,828 +G 5.7% coding (726/939 nt) fruK ← 1‑phosphofructokinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,212,8281.→G 5.7% 3.3 ‑0.7 ‑0.335coding (726/939 nt)fruK1‑phosphofructokinase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (18/15):  total (20/15)
Fisher's exact test strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.07e-01

CCTGCGCCAACGGTGCTTACGACATCGACTGACGGT.GGTTTGGCGATCCATTCGCCGGAGGCATTAACCCA  >  REL606/2212793‑2212863
                                    |                                   
ccTTCGCCAACGGTGCTTACGACATCGACTGACggt                                      <  1:1753990/36‑1 (MQ=255)
     gcCAACGGTGCTTACGACATCGACTGACGGG.GGctt                                >  1:1872964/1‑36 (MQ=37)
        aaCGGTGCTTACGACATTGACTGACGGTGGGTTTgg                              >  1:3318199/1‑36 (MQ=25)
        aaCGGTGCTTACGACATCGACTGACGGTGGGTTTgg                              >  1:3311388/1‑36 (MQ=37)
        aaCGGTGCTTACGACATCGACTGACGGT.GGTTTTGc                             >  1:3233811/1‑36 (MQ=255)
        aaCGGTGCTTACGACATCGACTGACGGT.GGTTTGGc                             <  1:334721/36‑1 (MQ=255)
        aaCGGTGCTTACGACATCGACTGACGGT.GGTTTGGc                             >  1:2078578/1‑36 (MQ=255)
         aCGGTGCTTACGACATCGACTGACGGT.GGGTTGGCg                            >  1:1453988/1‑36 (MQ=255)
           ggTGCTTACGACATCGACTGACGGT.GGTTTGGCGAt                          >  1:557512/1‑36 (MQ=255)
           ggTGCTTACGACATCGACTGACGGT.GGTTTGGCGAt                          >  1:2046596/1‑36 (MQ=255)
           ggTGCTTACGACATCGACTGACGGT.GGTTTGGCGAt                          >  1:4603572/1‑36 (MQ=255)
            gTGCTTACGACATCGACTGACGGT.GGTTTGGCGATc                         >  1:5126436/1‑36 (MQ=255)
             tGCTTACGACATCGACTGACGGT.GGTTTGGCGAttc                        >  1:4439783/1‑34 (MQ=255)
              gCTTACGACATCGACTGACGGT.GGTTTGGCGATCCa                       <  1:3897908/36‑1 (MQ=255)
               cTTACGACATCGACTGACGGT.GGTTTGGCGATCCAt                      <  1:4179603/36‑1 (MQ=255)
                ttACGACATCGACTGACGGT.GGTTTGGCGATCCAtt                     <  1:5182405/36‑1 (MQ=255)
                ttACGACATCGACTGACGGT.GGTTTGGCGATCCAtt                     <  1:4934865/36‑1 (MQ=255)
                ttACGAAATCGACTGACGGT.GGCTTTGCGATCCAtt                     >  1:2214955/1‑36 (MQ=25)
                 tACGAGATCGGCAGACGGT.GGTTTGGCGCTCCATTc                    <  1:4475381/36‑1 (MQ=19)
                 tACGACATCGACTGACGGT.GGTTTGGCGATCCATTc                    <  1:165000/36‑1 (MQ=255)
                 tACGACATCGACTGACGGG.GGGTTAGCGATACactt                    >  1:4721251/1‑33 (MQ=255)
                  aaaaCATCGCCTGACGGT.GGTTTGGCGATCCATTCg                   <  1:2230231/33‑1 (MQ=37)
                  aCGACATCGACTGACGGT.GGTTTGGCGATCCATTCg                   >  1:724716/1‑36 (MQ=255)
                  aCGACATCGACTGACGGC.GGGTTTGCGATTCATTcc                   >  1:925818/1‑35 (MQ=16)
                   cGACATCGACTGACGGT.GGTTTGGCGATCCATTCGc                  >  1:1723583/1‑36 (MQ=255)
                     acctCGACTGACGGT.GGTTTGGCGATCCATTCGCCg                <  1:2949826/33‑1 (MQ=255)
                     aCATCGACTGACGGT.GGTTTGGCGATCCATTAGCCg                >  1:4166535/1‑36 (MQ=255)
                        tCGACTGACGGT.GGTTTGGCGATCCATTCGCCGGAg             >  1:3518339/1‑36 (MQ=255)
                          gACTGGCGGT.GGTTTGGCGATCCATTCGCCGGAGGc           <  1:1205987/36‑1 (MQ=255)
                           aCTGACGGT.GGTTTGGCGAGCCATTCGCCGGAGGCa          >  1:4766302/1‑36 (MQ=255)
                            cTGACGGT.GGTTTGGCGATCCATTCGCCGGAGGCAt         <  1:4004247/36‑1 (MQ=255)
                            cTGACGGT.GGTTTGGCGATCCATTCGCCGGAGGCAt         <  1:3750712/36‑1 (MQ=255)
                             tGACGGT.GGTTTGGCGATCCATTCGCCGGAGGCAtt        <  1:239152/36‑1 (MQ=255)
                              gACGGT.GGTTTGGCGATCCATTCTCCCGTGGCGTTa       >  1:3025190/1‑36 (MQ=16)
                               gcGGT.GGTTTGGCGATCCATTCGCCGGAGGCATTaa      <  1:2095414/35‑1 (MQ=255)
                                cGGT.GGTATGGCGATCCATTCGCCGGAGGCATTAAc     <  1:4690217/36‑1 (MQ=255)
                                 ggt.ggtTTGGCGATCCATTCGCCGGAGGCATTAAcc    <  1:166637/36‑1 (MQ=255)
                                   t.ggtTTGGCGATCCATTCGCCGGAGGCATTAACCCa  >  1:5082355/1‑36 (MQ=255)
                                    |                                   
CCTGCGCCAACGGTGCTTACGACATCGACTGACGGT.GGTTTGGCGATCCATTCGCCGGAGGCATTAACCCA  >  REL606/2212793‑2212863

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG