Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,432,540 C→A 19.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→A19.9% 2.7 ‑4.1 ‑0.029intergenic (‑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/6):  ref base (14/9):  total (14/15)
Fisher's exact test strand distribution p-value = 1.69e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.82e-01

CATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGC  >  REL606/1432505‑1432570
                                   |                              
caTCAGCGTAATGAATGTTTTGATCAAACAGAGGGc                                >  1:6355107/1‑36 (MQ=255)
  tCAGCGTAATGAATGTCTTGATCAAACAGAGGGCaa                              <  1:2672019/36‑1 (MQ=255)
   cAGCGTAATGAATGTTTTGATCAAACAGAGGGCaaa                             <  1:39309/36‑1 (MQ=255)
   cAGCGTAATGAATGTTTTGATCAAACAGAGGGCaaa                             <  1:244574/36‑1 (MQ=255)
    aGCGTAATGAATGTTTTGATCAAACAGAGGGCaaaa                            >  1:1834965/1‑36 (MQ=255)
    aGCGTAATGAATGTTTTGATCAAACAGAGGGCaaaa                            <  1:590544/36‑1 (MQ=255)
    aGCGTAATGAATGTTTTGATCAAACAGAGGGCaaaa                            >  1:2712957/1‑36 (MQ=255)
      cGTAATGAATGTTTTGATCAAACAGAGGGCaaaaaa                          <  1:6642751/36‑1 (MQ=255)
      cGTAATGAATGTTTTGATCAAACAGAGGGCaaaaaa                          <  1:2323231/36‑1 (MQ=255)
        tAATGAATGTTTTGATCAAACAGAGGGCAAAAAAAt                        >  1:4821152/1‑36 (MQ=255)
        tAATGAATGTTTTGATCAAACAGAGGGCAAAAAAAt                        <  1:148892/36‑1 (MQ=255)
             aatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCa                   >  1:4520897/1‑36 (MQ=255)
             aatgTTTTGATAAAAAAGAGGGCAAAAAAATAGCCa                   <  1:2363041/36‑1 (MQ=37)
              atgTTTTGATCAAACAGAGGGCAAAAAAATAGCCAt                  <  1:2424275/36‑1 (MQ=255)
              atgTTTTGATCAAACAGAGGGAAAAAAAATAGCCAt                  <  1:6919774/36‑1 (MQ=255)
              atgTTTTGATAAAAAAGAGGGAAAAAAAATAGCCAt                  <  1:7443315/36‑1 (MQ=25)
               tgTTTTGATAAAAAAGAGGGAAAAAAAATAGCCATa                 <  1:3019820/36‑1 (MQ=25)
                gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATaa                >  1:2371099/1‑36 (MQ=255)
                gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATaa                >  1:5583793/1‑36 (MQ=255)
                gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATaa                >  1:825184/1‑36 (MQ=255)
                 ttttGATCAAACAGAGGGCAAAAAAATAGCCATata               >  1:3748907/1‑34 (MQ=255)
                  tttGATCAAACAGAGGGCAAAAAAATAGCCATaaaa              >  1:6107443/1‑36 (MQ=255)
                  tttGAAAAAAAAGAGGGAAAAAAAATAGCCATaaaa              <  1:660121/36‑1 (MQ=19)
                     gATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAg           >  1:3993996/1‑36 (MQ=255)
                       aaaaaAAGAGGAAAAAAAAATAGCCATAAAAAAGTa         <  1:4417474/34‑1 (MQ=21)
                          aaCAGAGGGCAAAAAAATAGCCATAAAAAAGTaaaa      >  1:5326389/1‑36 (MQ=255)
                          aaCAGAGGGCAAAAAAATAGCCATAAAAAAGTaaaa      >  1:1182222/1‑36 (MQ=255)
                           aaagagGGAAAAAAAATAGCCATAAAAAAGTaaaaa     <  1:445737/34‑1 (MQ=38)
                              gagGGCAAAAAAATAGCCATAAAAAAGTAAAAATGc  >  1:4676502/1‑36 (MQ=255)
                              gagGGCAAAAAAATAGCCATAAAAAAGTAAAAATGc  >  1:104031/1‑36 (MQ=255)
                                   |                              
CATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGC  >  REL606/1432505‑1432570

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