Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,574,203 A→T 32.2% 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→T32.2% 6.6 ‑2.1 ‑0.117intergenic (+45/+773)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (6/5):  total (12/5)
Fisher's exact test strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.97e-01

CAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTG  >  REL606/2574171‑2574238
                                |                                   
cAATAAAAAACCGCCGAATTTGGCGGTTTTTTTttt                                  >  1:3773867/1‑35 (MQ=255)
 aaTAAAAAACCGCCAAATTCGGCGGCTTTTTTtttt                                 >  1:1754050/1‑34 (MQ=17)
  aTAAAAAACCGCCGAATTTGGCGGGTTTTTGttttt                                >  1:4936253/1‑33 (MQ=14)
     aaaaaCCGCCGAATTTGGCGGCTTTTTTTTTTTAGt                             >  1:5105746/1‑36 (MQ=12)
            gccgCATTTGGCGGTTTTTTATTGCTAGTCTGGTTc                      <  1:2703443/36‑1 (MQ=255)
            gccgAATTTGGCGGTTTTTTTTTTTTTGTTTGGTTc                      >  1:1358190/1‑36 (MQ=9)
            gccgAATTTGGCGGTTTTTTTTTTCTATTTTTGTTc                      >  1:972590/1‑36 (MQ=12)
               gAATTTGGCGGTTTTTTATTTCTAGTCTGGTTcgcg                   >  1:216633/1‑36 (MQ=255)
                        ggTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca          <  1:1483135/36‑1 (MQ=255)
                          ttttttATTGCTAGTCTGGTTCGCGGCCTTTCcagc        >  1:529511/1‑36 (MQ=255)
                            ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag      <  1:3466720/36‑1 (MQ=255)
                            ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag      <  1:1421073/36‑1 (MQ=255)
                            ttttATTGCTAGTCTGGTTCGCGACCTTTCcagcag      >  1:2342438/1‑36 (MQ=255)
                             tttATTGCTAGTCTGGTTCGCGGCATTTCCAGCAgg     >  1:1969427/1‑36 (MQ=255)
                              ttATTGCTAGTCTGTTTCTCGGCCTTTCCAGCAGGt    <  1:3124950/36‑1 (MQ=18)
                              ttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt    >  1:1472300/1‑36 (MQ=255)
                               tATTGCTAGTCTGGTTCGCGGCCTTTCGAGCAGGtt   >  1:5892136/1‑36 (MQ=255)
                                aTTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg  >  1:3748747/1‑36 (MQ=255)
                                aTTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg  >  1:5010379/1‑36 (MQ=255)
                                |                                   
CAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTG  >  REL606/2574171‑2574238

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG