Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,192,978 G→C 100% intergenic (+146/+147) rrfE → / ← yjaB 5S ribosomal RNA/predicted acetyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,192,9780G→C100.0% 15.2 11intergenic (+146/+147)rrfE/yjaB5S ribosomal RNA/predicted acetyltransferase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (5/0):  total (5/6)
Fisher's exact test strand distribution p-value = 2.16e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATAAAGCG  >  REL606/4192964‑4193011
              |                                 
gTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCa              >  1:6304094/1‑36 (MQ=38)
 gggCGGACAGAACCCCCGCCATAAACTGCCAGGCAt             <  1:5450019/35‑1 (MQ=14)
  ggCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc            >  1:6484110/1‑36 (MQ=38)
    ggggAAGCACCCCCGCCATAAACTGCCAGGCATCaa          <  1:6414575/35‑1 (MQ=14)
     gggCAGGACGCCCGCCATAAACTGCCAGGCATCaaa         >  1:7055441/1‑36 (MQ=38)
     gggAAGAACCCCCGCCATAAACTGCCAGGCATCaaa         <  1:1574653/36‑1 (MQ=14)
      ggCAGGACGCCCGCCATAAACTGCCAGGCATCAAAt        >  1:2786987/1‑36 (MQ=38)
        cAGGACCACCGAGATAAACTGCCAGGCATCAAATaa      <  1:1708875/36‑1 (MQ=9)
         gagACCCCCGCCATAAACTGCCAGGCATCAAATaaa     <  1:572932/34‑1 (MQ=17)
           gACCCCCGCCATAAACTGCCAGGCATCAAATAAAGc   <  1:3661669/36‑1 (MQ=255)
            aCGCCCGCCATAAACTGCCAGGGATCAAATAAAGCg  >  1:5458434/1‑36 (MQ=255)
              |                                 
GTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATAAAGCG  >  REL606/4192964‑4193011

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 10 ≤  ATCG  < 25 ≤  ATCG  < 37 ≤  ATCG