Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,611,375 G→T 20.8% T116N (ACC→AAC)  ydfP ← conserved protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,611,3750G→T20.8% 4.4 ‑3.0 ‑0.024347ydfPconserved protein
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (5/14):  total (9/15)
Fisher's exact test strand distribution p-value = 4.74e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.16e-01

AGCATGTTGGTGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCCTCTGACGGCAGCGGCAAGATTTGCG  >  REL606/1611342‑1611410
                                 |                                   
ggCATGTTGGTGAGCATTCCTTCGGCGGTTTTGgtt                                   <  1:121464/35‑1 (MQ=255)
aGCATGTTGGTGAGCATTCCTTCGGCGGTTTTTgtt                                   >  1:469226/1‑36 (MQ=255)
aGCATGTTGGTGAGCATTCCTTCGGCGGTTTTGgtt                                   <  1:1338213/36‑1 (MQ=255)
aGCATGTTGGTGAGCATTCCTTCGGCGGTTTTGgtt                                   <  1:3886901/36‑1 (MQ=255)
aGCATGTTGGTGAGCATTCCTTCGGCGGGTTTTgtt                                   >  1:44641/1‑36 (MQ=37)
  cATGTTGGTGAGCATTCCTTCGGCGGTTTTGgttgt                                 <  1:744466/36‑1 (MQ=255)
  cATGTTGGTGAGCATTCCTTCGGCGGTTTTGgttgt                                 <  1:5921653/36‑1 (MQ=255)
   aTGTTGGTGAGCATTCCTTCGGCGGTTTTTgttttt                                >  1:4488777/1‑36 (MQ=37)
    tGTTGGTGAGCATTCCTTCGGCGGTTTTTGTTTttt                               >  1:5739153/1‑36 (MQ=37)
     gTTGGTGAGCATTCCTTCGGCGGTTTTTTTTTtttt                              >  1:4398162/1‑36 (MQ=25)
     gTTGGTGAGCATTCCTTCGGCGGTTTTTTTTTtttt                              >  1:3356461/1‑36 (MQ=25)
      ttGGTGAGCATTCCTTCGGCGGTTTTTTTTTttttt                             >  1:2606050/1‑35 (MQ=21)
        ggTGAGCATTCCTTCGGCGGTTTTTTTTTttttttt                           >  1:3673682/1‑33 (MQ=16)
        ggTGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGcc                           <  1:108248/36‑1 (MQ=255)
         gTGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCct                          <  1:5504005/36‑1 (MQ=255)
         gTGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCct                          <  1:2184168/36‑1 (MQ=255)
          tGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCctc                         <  1:3811153/36‑1 (MQ=255)
          tGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCctc                         <  1:5759432/36‑1 (MQ=255)
           gAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCctct                        <  1:5438365/36‑1 (MQ=255)
             gCATTCCTTCGGCGGTTTTGGTTGTTTTGCCTCTGa                      <  1:750653/36‑1 (MQ=255)
             gCATTCCTTCGGCGGTTTTGGTTGTTTTGCCTCTGa                      <  1:3331504/36‑1 (MQ=255)
               aTTCCTTCGGCGGTTTTGGTTGTTTTGCCTCTGACg                    <  1:2483966/36‑1 (MQ=255)
                  ccTTCGGCGGTTTTGGTTGTTTTGCCTCTGACGGCa                 <  1:2099170/36‑1 (MQ=255)
                  ccTTCGGCGGTTTTGGTTGTTTTGCCTCTGACGGCa                 <  1:4804783/36‑1 (MQ=255)
                    ttCGGCTGTTTTGGTTGTTTTGCCTCTGACGGCAGc               <  1:4005867/36‑1 (MQ=255)
                         cggTTTTGGTTGTTTTGCCTCTGACGGCAGCGGCaa          <  1:1286504/36‑1 (MQ=255)
                          ggTTTTGTTTGTTTTGCCTCTGACGGCAGCGGCAAg         <  1:4025949/36‑1 (MQ=255)
                           gTTTTGGTTGTTTTGCCTCTGACGGCAGCGGCAAGa        >  1:544313/1‑36 (MQ=255)
                           gTTTTGGTTGTTTTGCCTCTGACGGCAGCGGCAAGa        >  1:5039088/1‑36 (MQ=255)
                           gTTTTGGTTGTTTTGCCTCTGACGGCAGCCGCCAGa        >  1:2951669/1‑36 (MQ=37)
                           gTTTTGGTTGTTTTGCCTCTGACCGCAGGGGCAAgt        >  1:4159132/1‑35 (MQ=25)
                              ttGGTTGTTTTGCCTCTGACGGCAGCGGCAAGAttt     <  1:4583581/36‑1 (MQ=255)
                                ggTTGTTTTGCCTCTGACGGCAGCGGGAAGAttttc   >  1:2132212/1‑34 (MQ=37)
                                ggTTGTTTTGCCTCTGACGGCAGCGGCAAGATTTGc   >  1:3158258/1‑36 (MQ=255)
                                ggTTGTTTTGCCTCTGACGGCAGCGGCAAGATTTGc   <  1:1076336/36‑1 (MQ=255)
                                 gttgttTTGCCTCTGACGGCCGCGGCAAGATTTGCg  <  1:2748968/36‑1 (MQ=255)
                                 gttgttTTGCCTCTGAAGGCAGAGGCAAGATTTGCg  >  1:2633240/1‑36 (MQ=37)
                                 |                                   
AGCATGTTGGTGAGCATTCCTTCGGCGGTTTTGGTTGTTTTGCCTCTGACGGCAGCGGCAAGATTTGCG  >  REL606/1611342‑1611410

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG