Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,432,533 C→A 15.4% intergenic (‑143/‑131) ydbK ← / → ydbJ fused predicted pyruvate‑flavodoxin oxidoreductase: conserved protein/conserved protein/FeS binding protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,432,5330C→A15.4% 3.0 ‑3.8 ‑0.026intergenic (‑143/‑131)ydbK/ydbJfused predicted pyruvate‑flavodoxin oxidoreductase: conserved protein/conserved protein/FeS binding protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (15/7):  total (15/11)
Fisher's exact test strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.46e-01

TCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAA  >  REL606/1432498‑1432567
                                   |                                  
tCCCCCACATCAGCGTAATGAATGTTTTGATCAAAc                                    <  1:5441345/36‑1 (MQ=255)
     cacaTCAGCGTAATGAATGTTTTGATCAAACAGAgg                               <  1:1890922/36‑1 (MQ=255)
      acaTCAGCGTAATGAATGTTTTGATCAAACAGAggg                              >  1:3637147/1‑36 (MQ=255)
      acaTCAGCGTAATGAATGTTTTGATCAAACAGAggg                              >  1:3390384/1‑36 (MQ=255)
        aTCAGCGTAATGAATGTTTTGATCAAACAGAGGGCa                            >  1:892989/1‑36 (MQ=255)
         tCAGCGTAATGAATGTTTTGATCAAACAGAGGGCaa                           <  1:2810072/36‑1 (MQ=255)
           aGCGTAATGAATGTTTTGATCAAACAGAGGGCaaaa                         <  1:4324160/36‑1 (MQ=255)
            gCGTAATGAATGTTTTGATCAAACAGAGGGCaaaaa                        <  1:3749141/36‑1 (MQ=255)
              gTAATGAATGTTTTGATCAAACAGAGGGCaaaaaaa                      <  1:6403846/36‑1 (MQ=255)
                 atgaatgTTTTGATCAAACAGAGGGCAAAAAAATAg                   >  1:7421451/1‑36 (MQ=255)
                  tgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGc                  >  1:1082893/1‑36 (MQ=255)
                   gaatgTTTTGATCAAACAGAGGGCAAAAAAATAGcc                 <  1:6466513/36‑1 (MQ=255)
                   gaatgTTTCGATCAAACAGAGGGCAAAAAAATAGcc                 <  1:6228241/36‑1 (MQ=255)
                    aatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCa                >  1:508888/1‑36 (MQ=255)
                    aaggTTTTGATCAAACAGAGGGCAAAAAAATAGCCa                >  1:582921/4‑36 (MQ=255)
                       gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATaa             >  1:7205709/1‑36 (MQ=255)
                       gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATaa             >  1:5483200/1‑36 (MQ=255)
                       gTTTTGATAAAAAAGAGGGAAAAAAAATAGCCATaa             <  1:3669274/36‑1 (MQ=25)
                             aaaaaaAAGAGGCAAAAAAAATAGCCATAAAAAAGt       <  1:5205329/33‑1 (MQ=16)
                             aTCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGt       >  1:7727925/1‑36 (MQ=255)
                              tCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTa      >  1:2807549/1‑36 (MQ=255)
                              aaaaaAAGAGGAAAAAAAAATAGCCATAAAAAAGTa      <  1:3753493/34‑1 (MQ=21)
                               cAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTaa     >  1:1559925/1‑36 (MQ=255)
                               aaaaAAGAGGGAAAAAAAATAGCCATAAAAAAGTaa     <  1:5909126/35‑1 (MQ=37)
                                aaaCAGAGGGCAAAAAAATAGCCATAAAAAAGTaaa    >  1:2689463/1‑36 (MQ=255)
                                  aCAGAGGGCAAAAAAATAGCCATAAAAAAGTaaaaa  >  1:2544874/1‑36 (MQ=255)
                                  aCAGAGGGCAAAAAAATAGCCATAAAAAAGAaaaaa  >  1:2233168/1‑36 (MQ=255)
                                   |                                  
TCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAA  >  REL606/1432498‑1432567

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG