Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,574,191 T→C 15.1% intergenic (+33/+785) 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,1910T→C15.1% 2.2 ‑1.6 ‑0.820intergenic (+33/+785)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (7/10):  total (10/10)
Fisher's exact test strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.71e-01

TAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTT  >  REL606/2574156‑2574226
                                   |                                   
tAAGTATTTTACAGTCAATAAAAAACCGCCGAAttt                                     >  1:5498155/1‑36 (MQ=255)
   gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                  >  1:7888369/1‑36 (MQ=255)
   gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                  >  1:2554472/1‑36 (MQ=255)
    tATTTTAAAGGCAAAAAAAAACCGCCGAATTTGGCg                                 <  1:7864226/36‑1 (MQ=25)
         tACAGGCAATAAAAAACCGCCGAATTTGGCGGtttt                            >  1:200515/1‑36 (MQ=255)
          aaaGGCAATAAAAAACCGCCGAATTTGGCGGttttt                           >  1:3231099/3‑36 (MQ=255)
            aGGCAATAAAAAACCGCCGAATTTGGCGGttttttt                         >  1:365500/1‑35 (MQ=255)
             agCAATAAAAAACCGCCAAATTCGGCGGgtgttttt                        >  1:5210470/2‑34 (MQ=18)
                aaaaaaaaaCCGCCGAATTTGGCGGTTTTTTATTGc                     <  1:1390388/33‑1 (MQ=255)
                aaTAAAAAACCGCCGAATTTGGCGGTTTTTTTtttt                     >  1:5195756/1‑34 (MQ=17)
                aaTAAAAAACCGCCAAATTCGGCGGTTTTTTTtttt                     >  1:2333750/1‑34 (MQ=17)
                aaTAAAAAACAGCCAAATTCGGCGGTTTTTTTttat                     >  1:4369436/1‑34 (MQ=17)
                  aaaaaacCCGCCGAATTTGGCGGTTTTTTATTGCTa                   <  1:7322680/35‑1 (MQ=255)
                     aaatCCGCCGAATTTGGCGGTTTTTTATTGCTAGTc                <  1:6068504/36‑1 (MQ=255)
                     aaacCCGCCGAATTTGGCGGTTTTTTATTGCTAGTc                <  1:6477024/36‑1 (MQ=255)
                     aaaaCCGCCGAATTTGGCGGCGTTTTTTTTTTAGtg                >  1:6110003/1‑35 (MQ=12)
                        cccgccgAATTTGGCGGTTTTTTATTGCTAGTTTgg             <  1:1411216/35‑1 (MQ=255)
                           gccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTc          <  1:6634153/36‑1 (MQ=255)
                             cgAATATGGCGGTTTTTTATTGCTAGTCTGGTTcgc        <  1:7285214/36‑1 (MQ=255)
                               aaTTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCgg      <  1:5592711/36‑1 (MQ=255)
                                aTTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc     <  1:7124829/36‑1 (MQ=255)
                                  ttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt   <  1:6616335/36‑1 (MQ=255)
                                   tGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCtt  <  1:6368869/36‑1 (MQ=255)
                                   |                                   
TAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTT  >  REL606/2574156‑2574226

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG