Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,574,203 A→T 21.3% intergenic (+45/+773) sseA → / ← sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,574,2030A→T21.3% 2.4 ‑2.9 0.025intergenic (+45/+773)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (10/9):  total (16/9)
Fisher's exact test strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTG  >  REL606/2574171‑2574238
                                |                                   
cAATAAAAAACCGCCGAATTTGGCGGGTTTTTTttt                                  >  1:7469975/1‑35 (MQ=17)
cAATAAAAAACCGCCAAATTCGGCGGGTTTTTTttt                                  >  1:2541716/1‑35 (MQ=17)
 aaTAAAAAACCGCCAAATTCGGCGGTGTTTTTtttc                                 >  1:4236733/1‑34 (MQ=17)
  aaaaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCt                                <  1:7911174/34‑1 (MQ=255)
  aTAAAAAACCGCCGAATTTGGCGGTTTTTTTttttt                                >  1:7237136/1‑33 (MQ=17)
  aTAAAAAACCGCCGAAGTTGGCGGTTTTTTTtttgt                                >  1:3714005/1‑33 (MQ=17)
         cccgccgAATTTGGCGGTTTTTTATTGCTAGTCTgg                         <  1:8139377/35‑1 (MQ=255)
             ccgAATTTGGCGGTTTTTTATTTCTATTCTTGTTct                     >  1:5564722/1‑35 (MQ=21)
               gAATTTGGCGGTTTTTTTTTTCTATTCTTGTTtgcg                   >  1:1087243/1‑36 (MQ=12)
               gAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                   <  1:2668843/36‑1 (MQ=255)
                 aTTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                 <  1:4913446/36‑1 (MQ=255)
                 aTTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                 >  1:4931202/1‑36 (MQ=255)
                      gCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc            >  1:3319547/1‑36 (MQ=255)
                      gCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc            >  1:4602062/1‑36 (MQ=255)
                      gCGGTTTTTTATTGCTAGTCTGGTTAGCGGCCTTTc            >  1:6405121/1‑36 (MQ=255)
                       cGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc           <  1:8193421/36‑1 (MQ=255)
                        ggTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca          >  1:3707694/1‑36 (MQ=255)
                        ggTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca          >  1:3871584/1‑36 (MQ=255)
                        ggTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca          <  1:1216169/36‑1 (MQ=255)
                         gTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag         >  1:183084/1‑36 (MQ=255)
                          ttttttATTGCTAGTCTGGTTCGCGGCCTTTCcagc        <  1:7262536/36‑1 (MQ=255)
                            ttttATTGCTAGTCTGTTTCGCGGCCTTTCcagcag      >  1:2490225/1‑36 (MQ=255)
                            ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag      <  1:4121660/36‑1 (MQ=255)
                             tttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAgg     <  1:2062786/36‑1 (MQ=255)
                               tATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGtt   >  1:3525523/1‑36 (MQ=255)
                                aTTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg  <  1:56675/36‑1 (MQ=255)
                                |                                   
CAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTG  >  REL606/2574171‑2574238

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG