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

AATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTT  >  REL606/2574172‑2574237
                               |                                  
aaaaaaaacCCGCCGAATTTGGCGGTTTTTTATTGc                                <  1:5114470/33‑1 (MQ=17)
 aaaaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCt                               <  1:2842673/34‑1 (MQ=255)
 aTAAAAAACCGCCAAATTCGGCGGTTTTTTTttttt                               >  1:5484820/1‑33 (MQ=17)
 aTAAAAAACCGCAGAATTTGGCGGGTTTTTTttttt                               >  1:1829118/1‑33 (MQ=12)
  tAAAAAACCGCCGAATTTGGCGGTTTTTTTTTTCtt                              >  1:6955449/1‑35 (MQ=17)
       acccgccgAATTTGGCGGTTTTTTATTGCTAGTCTg                         <  1:751042/34‑1 (MQ=255)
          cgccgAATTTGGCGGTTTTTTTTTTTTTGTCTTGtt                      >  1:4590479/1‑36 (MQ=9)
                aTTTGGCGGTTTTTTATTTCTAGTCTTGTTCTCGGc                >  1:5905266/1‑36 (MQ=25)
                 tttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc               <  1:2382910/36‑1 (MQ=255)
                  ttGGCGGTTTTTTTTTTCTAGTCTGGTTCGCGGCCt              >  1:2812338/1‑36 (MQ=18)
                  ttGGCGGTTTTTTATTGCTAGTCTGTTTCGCGGCCt              <  1:3390390/36‑1 (MQ=255)
                     gCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc           <  1:6155738/36‑1 (MQ=255)
                      cGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc          >  1:3109828/1‑36 (MQ=255)
                       ggTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca         >  1:3635773/1‑36 (MQ=255)
                           ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag     >  1:3001264/1‑36 (MQ=255)
                           ttttATTGCTAGTCTGGTTCGCGGCCTTTCcagcag     >  1:6958167/1‑36 (MQ=255)
                             ttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt   >  1:3602673/1‑36 (MQ=255)
                             ttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt   >  1:4831840/1‑36 (MQ=255)
                             ttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt   <  1:6048503/36‑1 (MQ=255)
                             ttATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt   <  1:6530063/36‑1 (MQ=255)
                              tATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGCtt  >  1:187616/1‑36 (MQ=255)
                               |                                  
AATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTT  >  REL606/2574172‑2574237

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