Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,222,410 T→C 20.4% intergenic (+27/+56) ycgH → / ← ymgD predicted protein; b1170(b4491)/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,222,4100T→C20.4% 3.0 ‑4.5 ‑0.024intergenic (+27/+56)ycgH/ymgDpredicted protein; b1170(b4491)/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (13/6):  total (13/11)
Fisher's exact test strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.16e-01

CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCA  >  REL606/1222376‑1222443
                                  |                                 
cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATc                                  <  1:4037421/36‑1 (MQ=255)
 gCAGTGATGGTGGGCCTGAAATATAAGTTCTGATcc                                 <  1:4461768/36‑1 (MQ=255)
  cAGTGATGGTGGGCCTGAAATATAAGTTCTGATccc                                <  1:4107294/36‑1 (MQ=255)
     tgctggtggGCCTCAAATATAAGTTCTGATCccgcc                             <  1:324689/33‑1 (MQ=38)
      gATGGTGGGCCTGAAATATAAGTTCTGATCccgccg                            >  1:4981804/1‑36 (MQ=255)
         ggtggGCCTGAAATATAAGTTCTGATCCCGCCGTTa                         >  1:4766554/1‑36 (MQ=255)
             ggCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTa                     >  1:949322/1‑36 (MQ=255)
              gCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                    >  1:877822/1‑36 (MQ=255)
              gCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaa                    >  1:4289663/1‑36 (MQ=255)
               aatGACATAGAATTGCTGTTGCCGCCGTTAGCTaaa                   <  1:524148/34‑1 (MQ=0)
                cTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaa                  >  1:94318/1‑36 (MQ=255)
                cTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaa                  >  1:5742050/1‑36 (MQ=255)
                cTGAAATATAAGTTCTGATCCCGCCGCTAGCTaaaa                  >  1:4204477/1‑36 (MQ=255)
                  aaaaAATAATTCCTGACCCCGCCGTTAGCTaaaaaa                <  1:424789/35‑1 (MQ=16)
                     atatAAGTTCTGATCCCGCCGTTAGCTAAAAAACcg             >  1:1739487/1‑36 (MQ=255)
                       ataAGTTCTGACCCCGCCGTTAGCTAAAAAACcgcg           <  1:5690916/36‑1 (MQ=255)
                        taAGTTTTGACCCCGCCGTTAGCTAAAAAACCGcgt          <  1:5855342/36‑1 (MQ=38)
                         aatttCTGTCCCCGCCGTTAGCTAAAAAACCGcgtc         <  1:4573575/33‑1 (MQ=37)
                         aatttCAGACCCCGCCGTTAGCTAAAAAACCGcgtc         <  1:1940237/33‑1 (MQ=37)
                         aaGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc         >  1:4112821/1‑36 (MQ=255)
                            ttCTGATCCCGCCGTTGGCTAAAAAACCCCGTCGTa      >  1:2594522/1‑36 (MQ=38)
                             tCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt     >  1:997190/1‑36 (MQ=255)
                               tGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc   <  1:1086919/36‑1 (MQ=255)
                                gATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCa  >  1:5808772/1‑36 (MQ=255)
                                  |                                 
CGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCA  >  REL606/1222376‑1222443

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 16 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG