Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,432,540 C→A 16.9% intergenic (‑150/‑124) 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,5400C→A16.9% 6.1 ‑3.9 ‑0.224intergenic (‑150/‑124)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 (14/6):  total (14/10)
Fisher's exact test strand distribution p-value = 1.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.59e-01

ATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATG  >  REL606/1432506‑1432569
                                  |                             
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:5483200/1‑36 (MQ=255)
               gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATaa               >  1:7205709/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)
                            agagGGCAAAAAAATAGCCATAAAAAAGTAAAAATg  >  1:4780022/1‑36 (MQ=255)
                                  |                             
ATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATG  >  REL606/1432506‑1432569

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