Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,047,519 +G 7.7% coding (241/330 nt) yccK ← predicted sulfite reductase subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,047,5191.→G 7.7% 3.2 ‑3.5 ‑0.026coding (241/330 nt)yccKpredicted sulfite reductase subunit
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (3/21):  total (5/21)
Fisher's exact test strand distribution p-value = 3.08e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

CGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCT.GTTACCTTTTTCTTCGCCAAATTTATTCGCCATCG  >  REL606/1047484‑1047554
                                    |                                   
cGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCt                                      >  1:1016782/1‑36 (MQ=255)
 ggACCCTTCGGAACCAGTCGGTATAAATAGCTGCtc                                     <  1:1737639/36‑2 (MQ=18)
  gACCTTTCGGAAACAGTCGGTATAAATAGCGGGT.gg                                   >  1:6080513/1‑35 (MQ=255)
  gACCTTTCGGAAACAGTCGGTATAAATAGCGGCT.Gt                                   <  1:1355338/36‑1 (MQ=255)
      tttCGGAAACAGTCGGTATAAATAGCGGCT.GTTAcc                               <  1:4023971/36‑1 (MQ=255)
         cGGCAACAGTCGGTATAAATAGCGGCT.GTTACCttt                            <  1:5328943/36‑1 (MQ=255)
           gAAACAGTCGGTATAAATAGCGGCT.GTTACCttttt                          <  1:2946434/36‑1 (MQ=255)
           gAAACAGTCGGTATAAATAGCGGCT.GTTACCttttt                          <  1:5482438/36‑1 (MQ=255)
            aaccAGTCGGTATAAATAGCGGCT.GTTACCTTTttc                         <  1:3006932/33‑1 (MQ=255)
            aaaCAGTCGGTATAAATAGCGGCT.GTTACCTTTttc                         <  1:7815002/36‑1 (MQ=255)
            aaaCAGCCGTTATACATCTCGGCT.GTTACCTTTttc                         <  1:4603794/36‑1 (MQ=12)
             aaCAGTCGGTATAAATAGCGGCT.GTTACCTTTttct                        <  1:7638859/36‑1 (MQ=255)
              aCAGTCGGTATAAATAGCGGCTGGTTACTttttttt                        >  1:1314000/1‑34 (MQ=25)
               cAGTCGGTATAAATAGCGGCT.GTTACCTTTttcttc                      <  1:3273297/36‑1 (MQ=255)
               cAGTCGGTATAAATAGCGGCT.GTTACCTTTttcttc                      <  1:3813398/36‑1 (MQ=255)
               cAGTCGGTATAAATAGCGGCT.GTTACCTTTttcttc                      <  1:1376960/36‑1 (MQ=255)
                aGTCGGTATAAATAGCGGCT.GTTACCTTTTTCTTCg                     <  1:5708409/36‑1 (MQ=255)
                      tataAATAGCGGCT.GTTACCTTTTTCTTCGCCAAAt               <  1:886047/36‑1 (MQ=255)
                      tataAATAGCGGCT.GTTACCTTTTTCTTCGCCAAAt               <  1:5565231/36‑1 (MQ=255)
                      tataAATAGCGGCT.GTTACCTTTTTCTTCGCCAAAt               <  1:5402930/36‑1 (MQ=255)
                       ataAATAGCGGCTGGTTACTTTTTTTTTTTCCaatt               >  1:7184717/1‑34 (MQ=2)
                       ataAATAGCGGCT.GTTACCTTTTTCTTCGCCAAAtt              <  1:6382403/36‑1 (MQ=255)
                       ataAATAGCGGCT.GGTTCCTTTTTTTTTTCCAAAtt              >  1:3573312/1‑36 (MQ=16)
                         aaaTAGCGGCT.GTTACCTTTTTCTTCGCCAAATTTa            <  1:7532742/36‑1 (MQ=255)
                           aTAGCGGCT.GGTACCTTTTTTTTTTCCAAATTTTtt          >  1:1131262/1‑36 (MQ=16)
                            tAGCGGCT.GTTACCTTTTTCTTCGCCAAATTTATTc         <  1:6719130/36‑1 (MQ=255)
                              gCGGCT.GTTACCTTTTTCTTCGCCAAATTTATTCGc       <  1:2079218/36‑1 (MQ=255)
                               cGGCT.GTTACCTTTTTCTTCGCCAAATTTATTCGcc      >  1:7842956/1‑36 (MQ=255)
                                  cT.GTTACCTTTTTCTTCGCCAAATTTATTCGCCATc   <  1:3437423/36‑1 (MQ=255)
                                   t.GTTACCTTTTTCTTCGCCAAATTTATTCGCCATCg  <  1:7152441/36‑1 (MQ=255)
                                    |                                   
CGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCT.GTTACCTTTTTCTTCGCCAAATTTATTCGCCATCG  >  REL606/1047484‑1047554

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG